genebank biopython get /chromosome=”22″ or /map=”22q13.33″

it is located in the qualifiers field of the feature

from Bio import SeqIO

recs = SeqIO.parse("input.gb", format="genbank")

# Go over each record
for rec in recs:

    # Go over each feature
    for feat in rec.features:
        chr = feat.qualifiers.get("chromosome")
        if chr:
            print (chr)

prints:

 ['22']


Login
before adding your answer.

Traffic: 2021 users visited in the last hour

Read more here: Source link