[Bioperl-l] RemoteBlast

Suhail M suhail.mn12 at yahoo.co.uk
Fri Mar 18 12:18:49 EDT 2011


Dear Sir or Madam,
I am using a remote blast and encounter either a too many hsp or CPU overload 
problem with my specific gene. The gene itself is not that big so I was 
wondering is there a parameter I can pass when using the submit option to 
restrict number of hsp and to avoid the cpu and too many hsp problem? (I know I 
can use repeat masking and change the e-value but I need something that does not 
require those changes)

my $prog = 'blastn';
#my $db   = 'swissprot';
my $e_val= '1e-10';

my @params = ( '-prog'       => $prog,
#       '-data'       => $db,
               '-expect'     => $e_val,
               '-readmethod' => 'SearchIO'
            );

my $factory = Bio::Tools::Run::RemoteBlast->new(@params);

$Bio::Tools::Run::RemoteBlast::HEADER{'ENTREZ_QUERY'} = 'Homo sapiens [ORGN]';

delete $Bio::Tools::Run::RemoteBlast::HEADER{'FILTER'};

my $str = Bio::SeqIO->new(-file            =>'temporary.fa' , 
 '-format'         => 'fasta');

my $r = $factory->submit_blast('temporary.fa');

Thank you for your help.
Suhail


      



More information about the Bioperl-l mailing list