[Bioperl-l] One more odd little parsing problem for your list

Michael Muratet mam at torchconcepts.com
Sat May 31 15:00:45 EDT 2003


Greetings

I was parsing CDS features in Refseq human (hs.gbff.gz) when it died on
PSMAL/GCP III (NM_153696). The CDS was 527..1855. The error was from
Bio::PrimarySeq::subseq 'You have to start positive....sequence
[527:1855] Total III'. The parser in SeqIO is picking up the length from
the LOCUS line (as I recall) and for this record it sees 'III' and not
'1992' bp. It seems a lot to ask Bioperl to figure out every possible
configuration, maybe Genbank needs to have rules about whitespace in
gene names?

Cheers

Mike


More information about the Bioperl-l mailing list