[Bioperl-l] skip "exception"

Qunfeng qfdong at iastate.edu
Sun Sep 26 11:13:48 EDT 2004


Hi,

When I try to use Bioperl to parse genbank sequence file downloaded from 
NCBI, I am getting the following "exception" with sequence (acc# AJ635582). 
It seems to be the format problem in this genbank file.

My question is: is there any way in bioperl for users to test if there is 
any format problem in the sequence object and then decide by the users if 
they want to throw an exception or simply skips to the next object? Thanks!

Qunfeng

-------------------- WARNING ---------------------
MSG: exception while parsing location line [1163^1163] in reading 
EMBL/GenBank/SwissProt, ignoring feature variation (seqid=AJ6355
82):
------------- EXCEPTION  -------------
MSG: Only adjacent residues when location type is IN-BETWEEN. Not [1163] 
and [1163]
STACK Bio::Location::Simple::location_type 
/usr/lib/perl5/site_perl/5.8.0/Bio/Location/Simple.pm:280
STACK Bio::Location::Simple::new 
/usr/lib/perl5/site_perl/5.8.0/Bio/Location/Simple.pm:95
STACK Bio::Factory::FTLocationFactory::_parse_location 
/usr/lib/perl5/site_perl/5.8.0/Bio/Factory/FTLocationFactory.pm:218
STACK Bio::Factory::FTLocationFactory::from_string 
/usr/lib/perl5/site_perl/5.8.0/Bio/Factory/FTLocationFactory.pm:157
STACK (eval) /usr/lib/perl5/site_perl/5.8.0/Bio/SeqIO/FTHelper.pm:124
STACK Bio::SeqIO::FTHelper::_generic_seqfeature 
/usr/lib/perl5/site_perl/5.8.0/Bio/SeqIO/FTHelper.pm:123
STACK Bio::SeqIO::genbank::next_seq 
/usr/lib/perl5/site_perl/5.8.0/Bio/SeqIO/genbank.pm:501



More information about the Bioperl-l mailing list