[Bioperl-l] Problems using Bio::Tools::Blast
gert thijs
gert.thijs@esat.kuleuven.ac.be
Fri, 11 Aug 2000 15:56:41 +0200
Hello,
I'm trying to write a script to blast several genes with using webblast.pm.
The script should read a fasta file from stdin and blast each sequence
remotely at ncbi. Reading in the sequences is not a problem, but creating the
Blast object is much more trouble.
I assume I do something wrong in my script but I cannot find what or where.
Could anyone give me some advice. Code and error message are attached below.
This is my code sofar:
---
$instream = Bio::SeqIO->new(-fh=>\*STDIN, -format=>'Fasta');
%runParam = (
-prog => 'blastn',
-method => 'remote',
-version => 2,
-database => 'nr',
-html => 'off',
-filter => 'default',
-email => undef,
);
while ( $seq = $instream->next_seq() ) {
$runParam{-seqs} = [$seq];
# create a BlastObject
eval {
$blast = Bio::Tools::Blast->new(-run=>\%runParam, -parse=>1);
};
if($@) {
croak "\n*** Trouble creating Blast object:\n$@\n\n";
}
print $blast,"\n"; # just for checking
#$results = $blast->run(\%runParam);
# here should come some interpretation
}
---
This script returns the following exception:
-------------------- EXCEPTION --------------------
MSG: Blast file gi_609339_gb_U18675.blastn2.nr.temp.html is incomplete: can't
find program name.
Possibly an unrecognized report format.
CONTEXT: Error in uNKNOWN CONTEXT
SCRIPT: ./ParseBlastReport2.pl
STACK:
Bio::Tools::Blast::Run::Webblast::_blast(1175)
Bio::Tools::Blast::Run::Webblast::blast_remote(580)
Bio::Tools::Blast::_run_remote(1244)
Bio::Tools::Blast::run(1178)
Bio::Tools::SeqAnal::_initialize(279)
Bio::Root::Object::new(462)
main::./ParseBlastReport2.pl(33)
---------------------------------------------------
--
==========================================================
+ 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 +
==========================================================