[Bioperl-l] Warning while running bioperl StandAloneBlast.
shalabh sharma
shalabh.sharma7 at gmail.com
Tue May 19 14:28:11 UTC 2009
Hi Praveen,
Your script is bit confusing, i don't understand what you
actually trying to do.
But what i see is that you are parsing blast report without using SearchIO
module.
You have to include Bio::SearchIO in your script.
Hope this helps.
-Shalabh Sharma
On Tue, May 19, 2009 at 8:03 AM, Praveen Surendran <
praveen_spillai at yahoo.co.in> wrote:
> Hi,
>
> I am getting some warning along with results which says --> "[NULL_CAPTION]
> Warning: Failed to initialize search. ISAM Error code is -5"
> Can someone please let me know what this is all about.
>
> Please find my script below.
>
>
> ----------------------------------------------------------------------------------
> use Bio::SeqIO;
> use Bio::Tools::Run::StandAloneBlast;
> open FILE, ">./blast.txt";
> my $Seq_in = Bio::SeqIO->new (-file => $ARGV[0],
> -format => 'fasta');
> my $query = $Seq_in->next_seq();
> my $factory = Bio::Tools::Run::StandAloneBlast->new( 'program' =>
> 'blastn',
> 'database' => 'nr',
> _READMETHOD => "Blast" );
> my $blast_report = $factory->blastall($query);
> my $result = $blast_report->next_result;
> while( my $hit = $result->next_hit()) {
> print FILE "\thit name: ", $hit->name(), " significance: ",
> $hit->significance(), "\n";
> # while( my $hsp = $hit->next_hsp()) {
> # print "E: ", $hsp->evalue(), "frac_identical: ",
> $hsp->frac_identical(), "\n"; }}
> }
> close FILE
>
> -----------------------------------------------------------------------------------
>
> Kind Regards,
>
> Praveen Surendran.
>
>
> Explore your hobbies and interests. Go to
> http://in.promos.yahoo.com/groups/
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
More information about the Bioperl-l
mailing list