[BioPython] Loading dbxrefs from a gbk file

Peter biopython at maubp.freeserve.co.uk
Sat Oct 25 17:22:27 UTC 2008


On Sat, Oct 25, 2008 at 2:03 AM, Sebastian Bassi <sbassi at gmail.com> wrote:
> I have a genbank file like this one: http://www.pastecode.com.ar/f231664eb
> ...
> If you look at lines 130 to 133 (I highlighted in yellow) of the
> genbank sequence, there is cross database information (db_xref), but
> it is not associated with the SeqRecord, it is an empty list.

What you have highlighted is part of a gene feature, and would not be
part of the SeqRecord's db_xref list.  It should however be present in
the relevant SeqRecord feature.  Try:

print my_record.features[1]

(seeing as this is the second feature in the file, i.e. feature 1
using zero-based counting).

Peter



More information about the Biopython mailing list