[Biopython] User-defined SeqRecord annotations are trashed in INSDC formats?

[Uri Laserson]

If I load a GenBank-formatted record:

    a = SeqIO.parse('myfile.gb','gb').next()

then set an annotation:

    a.annotations['myannotation'] = 'saveme'

and then format the SeqRecord object as GenBank:

    a.format('gb')

then 'myannotation' is lost.

Is this expected behavior?  If so, that's a huge bummer...what is the
suggested method to store my own annotations in INSDC formats?

Uri



...................................................................................
Uri Laserson
Graduate Student, Biomedical Engineering
Harvard-MIT Division of Health Sciences and Technology
M +1 917 742 8019
laserson at mit.edu