[Bioperl-l] Running Blast
Lluis Armengol Dulcet
llarmengol@iro.es
Tue, 21 Nov 2000 15:44:08 +0100
Dear all,
I am experiencing some problems while trying to blast a Bio::Seq object
against databases using Bio::Tools::Blast . I do not actually have any
clue, I've tried all imaginable combinations of passing sequences to the
Bio::Seq object but no matter what I try I always get the following
error report:
> -------------------- EXCEPTION --------------------
> MSG: Blast failed: No file created or a error occured upon submission.
> CONTEXT: Error in uNKNOWN CONTEXT
> SCRIPT: bioperl6.pl
> STACK:
> Bio::Tools::Blast::Run::Webblast::_blast(1165)
> Bio::Tools::Blast::Run::Webblast::blast_remote(572)
> Bio::Tools::Blast::_run_remote(1207)
> Bio::Tools::Blast::run(1138)
> Bio::Tools::SeqAnal::_initialize(264)
> Bio::Root::Object::new(448)
> main::bioperl6.pl(15)
> ---------------------------------------------------
please find also enclosed the original script so you can take a look and
figure out what's going on ...
> #!/usr/local/bin/perl
>
> use Bio::Seq;
> use Bio::Tools::Blast;
>
> $seqObject = Bio::Seq->new( -file => '/root/Desktop/prova.tfa',
> -format => 'GenBank');
> %runParam = (
> -method => 'remote',
> -prog => 'blastn',
> -database => 'nr',
> -seqs => [$seqObject], # Bio::Seq.pm objects.
> );
>
> $blastObj = Bio::Tools::Blast->new( -run => \%runParam,
> -parse => 1,
> -signif => '1e-10',
> -strict => 1,
> );
Thanks a lot for your help!
Lluís