cd /path/to/CASE gzip -d swissprot.gz mv swissprot swissprot.fa
mkdir data mv swissprot.fa data/ mv protein_query.fasta data/
export PATH=$BLAST_TOP/c++/ReleaseMT/bin:$PATH makeblastdb -in data/swissprot.fa -dbtype prot
执行完成后会在“data”目录下生成以下4个文件,可使用ls命令查看。
{ time blastp -query data/protein_query.fasta -db data/swissprot.fa ; } 2>&1 |tee output_protein_alignments.txt
需要查看日志中最后的time命令输出中的“real”数值,单位是时间“x分钟y秒”,数值越低性能越优。