1hmmstat(1) HMMER Manual hmmstat(1)
2
3
4
6 hmmstat - display summary statistics for a profile file
7
8
9
11 hmmstat [options] <hmmfile>
12
13
14
16 The hmmstat utility prints out a tabular file of summary statistics for
17 each profile in <hmmfile>.
18
19
20
21 <hmmfile> may be '-' (a dash character), in which case profiles are
22 read from a <stdin> pipe instead of from a file.
23
24
25 The columns are:
26
27
28 idx The index of this profile, numbering each on in the file start‐
29 ing from 1.
30
31
32 name The name of the profile.
33
34
35 accession
36 The optional accession of the profile, or "-" if there is none.
37
38
39 nseq The number of sequences that the profile was estimated from.
40
41
42 eff_nseq
43 The effective number of sequences that the profile was estimated
44 from, after HMMER applied an effective sequence number calcula‐
45 tion such as the default entropy weighting.
46
47
48 M The length of the model in consensus residues (match states).
49
50
51 relent Mean relative entropy per match state, in bits. This is the
52 expected (mean) score per consensus position. This is what the
53 default entropy-weighting method for effective sequence number
54 estimation focuses on, so for default HMMER3 models, you expect
55 this value to reflect the default target for entropy-weighting.
56
57
58 info Mean information content per match state, in bits. Probably not
59 useful. Information content is a slightly different calculation
60 than relative entropy.
61
62
63 p relE Mean positional relative entropy, in bits. This is a fancier
64 version of the per-match-state relative entropy, taking into
65 account the transition (insertion/deletion) probabilities; it
66 may be a more accurate estimation of the average score contrib‐
67 uted per model consensus position.
68
69
70 compKL Kullback-Leibler distance between the model's overall average
71 residue composition and the default background frequency distri‐
72 bution. The higher this number, the more biased the residue
73 composition of the profile is. Highly biased profiles can slow
74 the HMMER3 acceleration pipeline, by causing too many nonhomolo‐
75 gous sequences to pass the filters.
76
77
78
80 -h Help; print a brief reminder of command line usage and all
81 available options.
82
83
84
86 See hmmer(1) for a master man page with a list of all the individual
87 man pages for programs in the HMMER package.
88
89
90 For complete documentation, see the user guide that came with your
91 HMMER distribution (Userguide.pdf); or see the HMMER web page ().
92
93
94
95
97 Copyright (C) 2015 Howard Hughes Medical Institute.
98 pFreely distributed under the GNU General Public License (GPLv3).
99
100 For additional information on copyright and licensing, see the file
101 called COPYRIGHT in your HMMER source distribution, or see the HMMER
102 web page ().
103
104
105
107 Eddy/Rivas Laboratory
108 Janelia Farm Research Campus
109 19700 Helix Drive
110 Ashburn VA 20147 USA
111 http://eddylab.org
112
113
114
115
116
117
118
119HMMER 3.1b2 February 2015 hmmstat(1)