[Bioperl-l] Doubt - Regarding Stand Alone Blast

[Dave Messina]

>
> formatdb -i sequences.fasta -p T -o T
>

The above command makes a database of the sequences.fasta file.



> blastall -p blastp -d homologene_result.fasta -i sequences.fasta -o stdout
> -a 1 -e 10
>

But then you ask for a database made from the homologene_result.fasta file.

For the above blastall command to work, you will have to use the following
formatdb command:

 formatdb -i homologene_result.fasta -p T -o T



Dave