介绍
使用hmmscan软件，将氨基酸序列比对pfam数据库
输入
序列文件：氨基酸序列，例如
>1_1 38720 921 tig00000462_pilon_pilon 38720 921
MKGSDGTSPNTTDAPIAVVGLSCRLPGAPDPATFRHLLRDGADAITEAPEGRWGGDTGAP
RRGGFLDRDRIDHFDAAFFGISPREAAAMDPQQRLTLELTWEALEDAGIVPDRLRDSRTG
VYIGVIADDYATLTRRG
E ：report models <span><</span>= this E-value threshold in output [10.0] (x>0)
domE ：report domains <span><</span>= this E-value threshold in output [10.0] (x>0)
结果
pfam.domtblout
例如，
# target name accession tlen query name accession qlen E-value score bias # of c-Evalue i-Evalue score bias from to from to from to acc description of target
#------------------- ---------- ----- -------------------- ---------- ----- --------- ------ ----- --- --- --------- --------- ------ ----- ----- ----- ----- ----- ----- ----- ---- ---------------------
ketoacyl-synt PF00109.27 253 1_1 - 12600 0 2392.7 48.5 1 8 3.2e-77 2.8e-74 249.9 0.1 2 253 14 254 13 254 0.93 Beta-ketoacyl synthase, N-terminal domain
