[Bioperl-l] Finding significant hits with BPlite

gert thijs gert.thijs@esat.kuleuven.ac.be
Tue, 22 Aug 2000 14:02:38 +0200


Hello,

I am using BPlite to parse a blast report, which works fine, but sometimes I
get the following exception using 
while ( my $hsp = $subject->nextHSP ){
  ...	
}

-------------------- EXCEPTION --------------------
MSG: asking for tag value that does not exist signif
CONTEXT: Error in uNKNOWN CONTEXT
SCRIPT: PrintSelectedSequences.pl
STACK: 
Bio::SeqFeature::Generic::each_tag_value(470)
Bio::SeqFeature::Similarity::_tag_value(116)
Bio::SeqFeature::Similarity::significance(134)
Bio::SeqFeature::SimilarityPair::significance(339)
Bio::Tools::BPlite::HSP::_initialize(81)
Bio::Root::Object::new(462)
Bio::Tools::BPlite::Sbjct::nextHSP(125)
BlastSequence::BlastSequence(80)
main::PrintSelectedSequences.pl(114)
---------------------------------------------------

If I look at the report the HSP looks alright to me. I have no idea how I can
work around this error and grab the HSP from the report. 

Can anyone give me some advice here?

Thanks,
Gert Thijs


-- 
==========================================================
+ Gert Thijs              gert.thijs@esat.kuleuven.ac.be +
+                                                        +
+ Dept. Elektrotechniek ESAT-SISTA                       +
+ Kardinaal Mercierlaan, 94                              +
+ B-3001 HEVERLEE  Belgium                               +
+ Tel :  +32-16-32 18 84 ---- Fax : +32-16-32 19 70      +
==========================================================