|||
第一步:建库 使用命令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
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-27 12:26
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社