blastn跟blastp用法（linux）
2022-05-06 本文已影响0人
geneonto
blastn
makeblastdb -in a.fasta -dbtype nucl -out a_db #建立参考核苷酸序列索引
blastn -query b.fasta -db a_db -evalue 1e-10 -outfmt 6 -max_target_seqs 6 -out result #-evalue为阈值，-outfmt为输出格式，-max_target_seqs为最大匹配基因对个数
blastp
makeblastdb -in a.fasta -dbtype prot -out a_db #建参考蛋白序列索引
blastp -query b.fasta -db a_db -evalue 1e-10 -outfmt 6 -out result -max_target_seqs 6 #比对分析
