[Biopython] read and write full ID line of EMBL SeqRecord?

Biocyberman biocyberman at gmail.com
Fri Feb 22 15:18:58 UTC 2013


Hi there,
I am using Biopython version 1.6.1 (latest).
My original ID line is:
ID   ACCESSION1; SV 1; linear; genomic DNA; HTG; PRO; 26402 BP.

But after reading and writing out, I got this:

ID   ACCESSION1; SV 1; ; DNA; ; PRO; 26402 BP.

How do I get the same ID line ?

Attached is the python script and input file.

Thanks for taking a look.
Biocyberman
-------------- next part --------------
A non-text attachment was scrubbed...
Name: input.embl
Type: application/octet-stream
Size: 1063 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/biopython/attachments/20130222/e7dfedbf/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: checkconvert.py
Type: application/octet-stream
Size: 249 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/biopython/attachments/20130222/e7dfedbf/attachment-0005.obj>


More information about the Biopython mailing list