[Biopython] Problem with parsing strand in Homo_sapiens.GRCh37.68 genbank files

Peter Cock p.j.a.cock at googlemail.com
Tue Aug 14 14:46:54 UTC 2012


On Tue, Aug 14, 2012 at 3:10 PM, Susan Wilson <smwilson at hpc.unm.edu> wrote:
> Hi,
>
> I am parsing the gb files with biopython. My problem is that none of the
> seqfeature.strand values are returning the plus strand (value == 1).

That should happen with a protein sequence.

> The commands below are a bit fabricated. (For instance, I have left out the
> opening and closing of fout.) I have read in
> Homo_sapiens.GRCh37.68.chromosome.1.dat using SeqIO.read.

What URL are you getting that file from?

Which version of Biopython are you using? There were some strand
related changes recently (internally moving it from the SeqFeature to
the SeqFeature's location object).

Thanks,

Peter



More information about the Biopython mailing list