blast

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
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • 杜鹃花开 文/卓文 那年春天,我跟随一群山民到山上开垦荒地。山岗上开满了红艳艳的...
    卓女阅读 575评论 5 8
  • 有那么一阵风 有那么一场雨 有那么一座城 金色的阳光 深沉的拥抱 胡杨挺立着坚强 猫咪惬意的微笑 老黑眺望着远方 ...
    水榭歌沨阅读 715评论 48 27
  • 1.完成采访任务 2. 少吃少吃
    谁要陪我吃火锅阅读 166评论 0 1
  • 我真正开始思考有关婚姻和家庭是在看了《小猪佩奇》动画片之后。 猪妈妈永远都很温柔,对待两个孩子既有耐心又有策略,还...
    EvaHuang阅读 413评论 0 0
  • 前天花大价钱买了口铁锅,因为之前的锅用久了锅底涂层被我刷掉了些,炒菜油烟比较大,炒菜还特耗油!叫人清洗油烟机,换锅...
    梦扬云子阅读 496评论 5 1

友情链接更多精彩内容