[Bioperl-l] Doubt - Regarding Stand Alone Blast

Dave Messina David.Messina at sbc.su.se
Tue May 3 09:16:33 UTC 2011


>
> 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



More information about the Bioperl-l mailing list