blast

作者: 琼脂糖 | 来源:发表于2017-11-06 09:31 被阅读16次

1.

-dbtype <String, nucl',prot'>
Molecule type of target db
-parse_seqids
Option to parse seqid for FASTA input if set, for all other input types
seqids are parsed automatically

2.

# Build the blast index.
makeblastdb -in index/all-proteins.fa -dbtype prot -out index/all -parse_seqids

# Run the blast queries.
blastp -db index/all -query L.fa

# Customize the output even more.
blastp -db index/all -query L.fa -outfmt "7 qseqid sseqid pident"

# Three best alignments
blastp -db index/all -query L.fa -outfmt "6 qseqid sseqid pident" | sort -k3 -rn | head -5

# The worst alignments.
blastp -db index/all -query L.fa -outfmt "6 qseqid sseqid pident" | sort -k3 -rn | tail -5

outfmt 6为m8格式。

3.

fasta之外是makedb的输出。


db

相关文章

网友评论

      本文标题:blast

      本文链接:https://www.haomeiwen.com/subject/hzvkmxtx.html