linux blast

建库

减压后,改名为blast,并在blas目录在建立db文件
1,建立数据库
makeblastdb -in db.fasta -dbtype nucl(prot) -parse_seqids -hash_index -out dbname
2、比对

blastn为例:
blastn -query tese.fa -db dbname -out test.name -evalue 1e-5 -outfmt 6 -num_threads 8 -max_target_seqs 5 #6是为一行一行的格式.7为带注释行的tab格式的输出 -num_threads 8 是用多少cup,-max_target_seqs 5 设置最多的的目标序列匹配数


blastall -i query.fasta -d database -o balast.out -p blastn -e 1e-10 -m 6

猜你喜欢

转载自www.cnblogs.com/zhanmaomao/p/9758256.html
今日推荐