#!/usr/bin/perl #sub search { use Bio::Tools::Run::StandAloneBlast; use Bio::SearchIO::blast; #my $query = Bio::Seq -> new ( -seq=>"$_[0]", # -id=>"query"); my $query=Bio::Seq->new(-seq=>"ctgtattctgggatgca"); my $len=$query->length(); #@db=('nt.nal'); #foreach my $db (@db) { my $factory = Bio::Tools::Run::StandAloneBlast->new( -program =>"blastn", -database =>'D:/blast/bin/ecoli.nt', -e =>1, -o=>'lclblastResult.txt'); my $rc = $factory->blastall($query);