[Bioperl-l] (no subject)
Parminder Kang
ppkang@hotmail.com
Tue, 12 Jun 2001 13:12:07
hi,
i am new to bioperl, trying to run local blast, but running into some
problems.
so far i have this
use Bio::Tools::Run::StandAloneBlast();
@params = ('program'=> 'blastn',
'database' => 'gpcr_nt.fasta');
my $factory = Bio::Tools::Run::StandAloneBlast->new(@params);
my $input = Bio::Seq->new('-id'=>"test query",
'-seq' =>"atggagcccaagcaaggaattctaaggg");
my $blast_report = $factory->blastall($input);
while(my $sbjct = $blast_report->nextSbjct) {
# print "$sbjct"
My problem is that when now I print the $sbjct, I only get Hash references.
May be I have to add something to the script I did not see anything else in
the documentation.
If anyone sees the what the problem is please let me know. If you have
examples scripts that would be helpful too. Also is this the best way to run
the local blast.
Thanks
PK
_________________________________________________________________
Get your FREE download of MSN Explorer at http://explorer.msn.com