[Bioperl-l] StandAlone BLAST error

Arvind Jayasundar Arvind.Jayasundar at portalplayer.com
Sun Feb 2 01:08:25 EST 2003


Hi
 
I am trying to blast sequences from a COG against a single sequence in a loop. After successfully blasting about 10 sequences, I get the following error. Can anybody throw some light on this ? Thanks in advance for any comments.
 
Arvind
 
------------- EXCEPTION  -------------
MSG:  No hit object found for bl2seq report
STACK Bio::Tools::BPbl2seq::next_feature C:/Perl/site/lib/Bio/Tools/BPbl2seq.pm:
239
STACK toplevel hw3.pl:83
--------------------------------------
Error removing C:\windows\TEMP\Jiu9H1ueKt at C:/Perl/lib/File/Temp.pm line 839,
<GEN26> line 36.
The offending code looks something like
...
      
my @params = ('program' => 'blastp');
my $factory = Bio::Tools::Run::StandAloneBlast->new(@params);
my $evalue = 0.0001;
$factory->e($evalue);

...
      $cogIO = Bio::SeqIO->new(-file => $cogfile, -format => 'fasta');
      # for every sequence in the cog file, blast against $myseq
      while ($seq = $cogIO->next_seq()) {
        if ($report = $factory->bl2seq($seq, $myseq)) {
         if ($hsp = $report->next_feature) {
           print "\n Hit:", $seq->display_id(), " ",
"e: ", $hsp->P(), " score: ". $hsp->score, "bits: ", $hsp->bits;
          }
        }
      }
...



More information about the Bioperl-l mailing list