[Bioperl-l] standalone blast on 2 databases

Damien Mattei Damien.Mattei@unice.fr
Thu, 30 May 2002 14:45:44 +0200


I want to do a standalone blast on 2 databases,
i know it is possible, as i can do it on the
command line specifying the name of databse with -d "nt htg" option
  for example,but when i do it in bio-perl i get an error:

@params = ('program'  => 'blastn',
             'database' => 'nt htg',
	    'outputfile' => 'blast_gmap.out' );

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

result in:
[blastall] ERROR: Arguments must start with '-' (the offending argument 
#5 was: 'htg')
------------- EXCEPTION  -------------
MSG: blastall call crashed: 256 /home/dmattei/ncbi/build/blastall -p 
blastn  -d  /home/dmattei/ncbi/db/nt htg  -i 
/home/dmattei/tmp/gi5622463  -o  blast1.out

is there a way to blast in both databases?

--------------------------------------
Damien Mattei
C.N.R.S / U.N.S.A - UMR 6549
Instabilite & Alterations des Genomes
Universite de Nice - Sophia Antipolis
Faculte de Medecine, Av. de Valombrose
06107 NICE Cedex 2 FRANCE
Phone: +33 (0)4 93 37 76 66
Fax:   +33 (0)4 93 53 30 71
mailto:mattei@unice.fr
http://www-iag.unice.fr/
--------------------------------------