[Bioperl-l] Standalone blast
florent.angly at gmail.com
Fri Apr 29 05:24:45 UTC 2011
To run BLAST the way Dave described, you need to have BLAST installed on
your computer, and you need to download BLAST databases to your computer
(or make them yourself with the formatdb command). There are plenty of
databases available on the NCBI FTP website: ftp://ftp.ncbi.nih.gov/.
And yes, some of these databases are very large and will take a long
time to download. By the way, the BLAST may also take a very long time
to execute if you use large databases, so, you'd better run the analysis
on a powerful computer or a server.
Also read this documentation:
It stipulates that you can BLAST an entire FASTA file (not just a
$inputfilename = 't/testquery.fa';
$blast_report = $factory->blastall($inputfilename);
On 29/04/11 14:53, khush ........ wrote:
> Dear Dave,
> Thank you for your support.
> If need to change the following lines like
> $blast_obj = Bio::Tools::Run::StandAloneBlast->new(-program => 'blastx',
> -database => 'nr.fa'));
> $seq_obj = Bio::Seq->new(-id =>"test query", -seq =>"file.fa");
> I have a simple and basic query for you, as I am beginners in bioperl, that
> if I need to download the whole nr database from NCBI to run the code or It
> will directly fetch information from the NCBI website. I do not understand
> it, because downloading the whole nr d/b itself takes long time for me.
> How could I read whole file instead of simple string "TTTATAGATAGAGACAG" in
> -seq (a fasta file). Is there a simple way to do the exercise according to
> my conditions.
> Thank you
> On Thu, Apr 28, 2011 at 12:59 PM, Dave Messina<David.Messina at sbc.su.se>wrote:
>> Hi Kamal,
>> This is covered in the beginners' HOWTO:
>> On Thu, Apr 28, 2011 at 07:22, khush ........<bioinfo.khush at gmail.com>wrote:
>>> I have some sequences ~250 and wanted to use BLASTX to blast against nr
>>> database of NCBI, as this is time consuming using web based search. Can
>>> one please tell me how to start BIOPERL with scuh problems. I know that
>>> is possible with bioperl, but do not know how.
>>> Any suggestion will be appreciable.
>>> Thanks in advance
>>> Bioperl-l mailing list
>>> Bioperl-l at lists.open-bio.org
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
More information about the Bioperl-l