« Return to Thread: blastall problem

blastall problem

by DeeGee () :: Rate this Message:

Reply to Author | View in Thread

hi all,
can anyone plz help me out with this problem that i've been dealing with for quite a while now. following is a part of my script that's not working for some reason. it is suppose to get the sequence from 'result/fasta.faa' and do the blast.

###my script ###########
......
my $Seq_in = Bio::SeqIO->new (-file => 'result/fasta.faa', '-format' => 'Fasta');
my $queryin = $Seq_in->next_seq();
my $factory = Bio::Tools::Run::StandAloneBlast->new('program'  => 'blastp',
                                                 'database' => '/export/home/database/nr',
                                                 _READMETHOD => 'Blast'
                                                  );
$factory->outfile("result/out.blast");
my $blastreport = $factory->blastall($queryin);
.....

when i paste the protein sequence into the textarea of my html page and save the same as 'result/fasta.faa', so that the above script would do the blast, i get the following error:

Software error:
------------- EXCEPTION  -------------
MSG:    not Bio::Seq object or array of Bio::Seq objects or file name!
STACK Bio::Tools::Run::StandAloneBlast::blastpgp /usr/perl5/5.6.1/lib/Bio/Tools/Run/StandAloneBlast.pm:611
STACK toplevel /usr/local/apache2/htdocs/remote_ncbi.pl:50
--------------------------------------
i would appreciate your help.
i would also like to add that the 'result/fasta.faa' has the sequence saved in it. i don't think that there's anything wrong with the open(OUTPUT,">result/fasta.faa"); line as i could get the 'fasta.faa' file with the sequence in it. it looks like the blast is not being able to read from the result/fasta.faa. ^ ^*

 « Return to Thread: blastall problem