[Bioperl-l] Trouble parsing PSI-BLAST
Rachita Sharma
rachitasharma at gmail.com
Tue Nov 3 22:12:11 UTC 2009
I am having trouble parsing PSI-BLAST results. Please help.
The code is:
my $in = new Bio::SearchIO( -format => 'blast',
-file =>
"BS_XFpsiRblastoutputs/e${ev}/bloutput${i}.txt");
while( my $result = $in->next_result ) {
while( my $hit = $result->next_hit ) {
$sth->execute($result->query_name, $hit->name, $hit->significance);
print "Query executed!\n";
}
}
The error is:
------------- EXCEPTION: Bio::Root::Exception -------------
MSG: no data for midline ***** No hits found ******
STACK: Error::throw
STACK: Bio::Root::Root::throw
/usr/lib/perl5/site_perl/5.8.8/Bio/Root/Root.pm:359
STACK: Bio::SearchIO::blast::next_result
/usr/lib/perl5/site_perl/5.8.8/Bio/SearchIO/blast.pm:1813
STACK: BSubVCpsiRblast.pl:92
-----------------------------------------------------------
More information about the Bioperl-l
mailing list