linux blast

时间:2023-03-09 02:27:33
linux blast

建库

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

还有formatdb    http://boyun.sh.cn/bio/?p=1483

建立核算库:  formatdb -i  fa -p F -o T

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