生物信息_Blast比对_blast

生物信息_Blast比对_blast

 

先对ref建库
/bin/formatdb -i /dir/ref.fa -p F -o
进行blast比对，输出是m8格式，具体很多参数，参看《常用生物数据分析软件》
/bin/blastall -i input.fa -d /dir/ref.fa -o blast.out -p blastn -m 8