suviu的个人博客分享 http://blog.sciencenet.cn/u/suviu

博文

windows下diamond快速比对蛋白质序列

已有 4246 次阅读 2020-6-9 19:09 |个人分类:软件|系统分类:科研笔记| diamond, 蛋白质比对, 超级快, windows

第一步:建库
使用命令diamond makedb --in 建库文件.fa --db 输出文件名称
E:\Util_file\diamond-windows>diamond makedb --in Brassica_napus.annotation_v5.pep_rename.fa --db Bra
ssica  #建库示例命令
diamond v0.9.32.133 (C) Max Planck Society for the Advancement of Science
Documentation, support and updates available at http://www.diamondsearch.org

#CPU threads: 8
Scoring parameters: (Matrix=BLOSUM62 Lambda=0.267 K=0.041 Penalties=11/1)
Database input file: Brassica_napus.annotation_v5.pep_rename.fa
Opening the database file...  [0s]
Loading sequences...  [0.382s]
Masking sequences...  [0.624s]
Writing sequences...  [0.116s]
Hashing sequences...  [0.015s]
Loading sequences...  [0s]
Writing trailer...  [0.007s]
Closing the input file...  [0s]
Closing the database file...  [0s]
Database hash = 349da09b950d300f557b96e24122a623
Processed 101040 sequences, 33719134 letters.
Total time = 1.156s

E:\Util_file\diamond-windows>diamond blastp --db Brassica -q Gmax_508_Wm82.a4.v1.protein_primaryTran
scriptOnly.fa -o BrassicaWithGmax.txt  #比对实例命令
diamond v0.9.32.133 (C) Max Planck Society for the Advancement of Science
Documentation, support and updates available at http://www.diamondsearch.org

#CPU threads: 8
Scoring parameters: (Matrix=BLOSUM62 Lambda=0.267 K=0.041 Penalties=11/1)
Temporary directory:
Opening the database...  [0s]
#Target sequences to report alignments for: 25
Reference = Brassica.dmnd
Sequences = 101040
Letters = 33719134




https://blog.sciencenet.cn/blog-2557039-1237161.html

上一篇:关于基因家族的期刊统计及内容浅析
下一篇:VBox虚拟机中,deepin系统下共享文件夹
收藏 IP: 221.232.142.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-19 22:14

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部