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

Peter Cock p.j.a.cock at googlemail.com
Fri Feb 22 16:08:13 UTC 2013


On Fri, Feb 22, 2013 at 3:18 PM, Biocyberman <biocyberman at gmail.com> wrote:
> 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

This is probably part of https://redmine.open-bio.org/issues/2578
(the GenBank and EMBL code overlaps a lot).

Peter



More information about the Biopython mailing list