What it does
The cmstat utility prints out a tabular file of summary statistics for each given covariance model.
By default, cmstat prints general statistics of the model and the alignment it was built from, one line per model in a
The columns are:
- The index of this profile, numbering each on in the file starting from 1.
- The name of the profile.
- The optional accession of the profile, or ”-” if there is none.
- The number of sequences that the profile was estimated from.
- The effective number of sequences that the profile was estimated from, after Infernal applied an effective sequence number calculation such as the default entropy weighting.
- The length of the model in consensus residues (match states).
- The expected maximum length of a hit to the model.
- The number of basepairs in the model.
- The number of bifurcations in the model.
- What type of model will be used by default in cmsearch and cmscan for this profile, either ”cm” or ”hmm”. For profiles with 0 basepairs, this will be ”hmm” (unless the --nohmmonly option is used). For all other profiles, this will be ”cm”.
- Mean relative entropy per match state, in bits. This is the expected (mean) score per con-
sensus position. This is what the default entropy-weighting method for effective sequence
number estimation focuses on, so for default Infernal, this value will often reflect the default
target for entropy-weighting. If the ”model” field for this profile is ”hmm”, this field will be ”-”.
- Mean relative entropy per match state, in bits, if the CM were transformed into an HMM (information from structure is ignored). The larger the difference between the CM and HMM
relative entropy, the more the model will rely on structural conservation relative sequence conservation when identifying homologs.
For further questions please refere to the Infernal Userguide.