Biotech > FAQ > BioPerl FAQ (Frequently Asked Questions)

How do I tell BLAST to search multiple database ...

To see other biotech frequently asked questions, please visit http://biotech.fyicenter.com/faq/

(Continued from previous question...)

How do I tell BLAST to search multiple database using Bio::Tools::Run::StandAloneBlast?

Put the names of the databases in a variable. like so:

my $dbs = '"/dba/BMC.fsa /dba/ALC.fsa /dba/HCC.fsa"';
my @params = ( d	      => "$dbs",
program     => "BLASTN",
_READMETHOD => "Blast",
outfile     => "$dir/est.bls" );

my $factory =
Bio::Tools::Run::StandAloneBlast->new(@params);
my $seqio = Bio::SeqIO->new(-file=>'t/amino.fa',-format => 'Fasta' );
my $seqobj = $seqio->next_seq();
$factory->blastall($seqobj);

(Continued on next question...)

Other Frequently Asked Questions