|
|
ehmmbuild |
Please help by correcting and extending the Wiki pages.
Usage:
ehmmbuild [options] alignfile hmmfile
Important note: the alignfile (input) and hmmfile (output) parameters are specified in the reverse order in the original HMMER.
hmmbuild reads a multiple sequence alignment file
Go to the input files for this example
Example 2
Go to the input files for this example
Example 3
Go to the input files for this example
Example 4
Go to the input files for this example
More or less all options documented as "expert" in the original hmmer user guide are given in ACD as "advanced" options (-options must be specified on the command-line in order to be prompted for a value for them).
ehmmbuild reads any normal sequence USAs.
The following additional options are provided:
Please read the 'Notes' section below for a description of the differences between the original and EMBASSY HMMER, particularly which application command line options are supported.
Although we take every care to ensure that the results of the EMBOSS
version are identical to those from the original package, we recommend
that you check your inputs give the same results in both versions
before publication.
Please report all bugs in the EMBOSS version to the EMBOSS bug team,
not to the original author.
Jon Ison (jison © ebi.ac.uk)
This program is an EMBASSY wrapper to a program written by Sean Eddy as part of his hmmer package.
Please report any bugs to the EMBOSS bug team in the first instance, not to Sean Eddy.
Algorithm
Please read the Userguide.pdf distributed with the original HMMER and included in the EMBASSY HMMER distribution under the DOCS directory.
Usage
Here is a sample session with ehmmbuild
% ehmmbuild globins50.msf globin.hmm -nhmm globins50 -strategy D
Build a profile HMM from an alignment.
hmmbuild - build a hidden Markov model from an alignment
HMMER 2.3.2 (Oct 2003)
Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment file: ../../data/hmmnew/globins50.msf
File format: MSF
Search algorithm configuration: Multiple domain (hmmls)
Model construction strategy: MAP (gapmax hint: 0.50)
Null model used: (default)
Prior used: (default)
Sequence weighting method: G/S/C tree weights
New HMM file: globin.hmm [appending]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment: #1
Number of sequences: 50
Number of columns: 308
Determining effective sequence number ... done. [2]
Weighting sequences heuristically ... done.
Constructing model architecture ... done.
Converting counts to probabilities ... done.
Setting model name, etc. ... done. [globins50]
Constructed a profile HMM (length 143)
Average score: 189.04 bits
Minimum score: -17.62 bits
Maximum score: 234.09 bits
Std. deviation: 53.18 bits
Finalizing model configuration ... done.
Saving model to file ... done.
//
hmmbuild -n globins50 --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F globin.hmm ../../data/hmmnew/globins50.msf
Go to the output files for this example
% ehmmbuild rrm.sto myhmms -nhmm rrm -strategy D
Build a profile HMM from an alignment.
hmmbuild - build a hidden Markov model from an alignment
HMMER 2.3.2 (Oct 2003)
Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment file: ../../data/hmmnew/rrm.sto
File format: Stockholm
Search algorithm configuration: Multiple domain (hmmls)
Model construction strategy: MAP (gapmax hint: 0.50)
Null model used: (default)
Prior used: (default)
Sequence weighting method: G/S/C tree weights
New HMM file: myhmms [appending]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment: #1
Number of sequences: 90
Number of columns: 104
Determining effective sequence number ... done. [88]
Weighting sequences heuristically ... done.
Constructing model architecture ... done.
Converting counts to probabilities ... done.
Setting model name, etc. ... done. [rrm]
Constructed a profile HMM (length 77)
Average score: 75.37 bits
Minimum score: 26.86 bits
Maximum score: 119.83 bits
Std. deviation: 21.78 bits
Finalizing model configuration ... done.
Saving model to file ... done.
//
hmmbuild -n rrm --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F myhmms ../../data/hmmnew/rrm.sto
Go to the output files for this example
% ehmmbuild fn3.sto myhmms -nhmm fn3 -strategy D
Build a profile HMM from an alignment.
hmmbuild - build a hidden Markov model from an alignment
HMMER 2.3.2 (Oct 2003)
Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment file: ../../data/hmmnew/fn3.sto
File format: Stockholm
Search algorithm configuration: Multiple domain (hmmls)
Model construction strategy: MAP (gapmax hint: 0.50)
Null model used: (default)
Prior used: (default)
Sequence weighting method: G/S/C tree weights
New HMM file: myhmms [appending]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment: fn3
Number of sequences: 108
Number of columns: 119
Determining effective sequence number ... done. [108]
Weighting sequences heuristically ... done.
Constructing model architecture ... done.
Converting counts to probabilities ... done.
Setting model name, etc. ... done. [fn3]
Constructed a profile HMM (length 84)
Average score: 67.54 bits
Minimum score: 27.53 bits
Maximum score: 105.77 bits
Std. deviation: 17.07 bits
Finalizing model configuration ... done.
Saving model to file ... done.
//
hmmbuild -n fn3 --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F myhmms ../../data/hmmnew/fn3.sto
Go to the output files for this example
% ehmmbuild pkinase.sto myhmms -nhmm pkinase -strategy D
Build a profile HMM from an alignment.
hmmbuild - build a hidden Markov model from an alignment
HMMER 2.3.2 (Oct 2003)
Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment file: ../../data/hmmnew/pkinase.sto
File format: Stockholm
Search algorithm configuration: Multiple domain (hmmls)
Model construction strategy: MAP (gapmax hint: 0.50)
Null model used: (default)
Prior used: (default)
Sequence weighting method: G/S/C tree weights
New HMM file: myhmms [appending]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Alignment: pkinase
Number of sequences: 67
Number of columns: 471
Determining effective sequence number ... done. [67]
Weighting sequences heuristically ... done.
Constructing model architecture ... done.
Converting counts to probabilities ... done.
Setting model name, etc. ... done. [pkinase]
Constructed a profile HMM (length 294)
Average score: 312.79 bits
Minimum score: 176.48 bits
Maximum score: 390.44 bits
Std. deviation: 44.60 bits
Finalizing model configuration ... done.
Saving model to file ... done.
//
hmmbuild -n pkinase --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F myhmms ../../data/hmmnew/pkinase.sto
Go to the output files for this example
Command line arguments
Where possible, the same command-line qualifier names and parameter order is used as in the original hmmer. There are however several unavoidable differences and these are clearly documented in the "Notes" section below.
Standard (Mandatory) qualifiers:
[-alignfile] seqset (Aligned) protein sequence set filename and
optional format, or reference (input USA)
-nhmm string Name for this HMM. The name can be any
string of non-whitespace characters (e.g.
one 'word'). There is no length limit (at
least not one imposed by HMMER; your shell
will complain about command line lengths
frst). (Any string without whitespace is
accepted)
-strategy menu [D] All alignments are local with respect to
the sequence and are configured to be local
(fragmentary) or global with respect to the
HMM. The model is also configured to find a
single or multiple domains (matches) to a
sequence. The options for configuring the
model are as follows: (D): The default
setting. Multiple domains per sequence,
global alignments with respect to the HMM.
(F): Multiple domains per sequence, local
alignments with respect to the HMM.
Analogous to the old hmmfs program of HMMER
1. (G) Single domain per sequence, global
alignment with respect to the HMM. Analogous
to the old hmms program of HMMER 1. (S)
Single domain per sequence, local alignments
with respect to the HMM. Analogous to the
old hmmsw program of HMMER 1. (Values: D
(global-multidomain); F (local-multidomain);
G (global-singledomain); S
(local-singledomain))
[-hmmfile] outfile [*.ehmmbuild] HMMER hidden markov model
output file
Additional (Optional) qualifiers: (none)
Advanced (Unprompted) qualifiers:
-prior infile Read a Dirichlet prior from file, replacing
the default mixture Dirichlet. The format of
prior files is documented in the User's
Guide, and an example is given in the Demos
directory of the HMMER distribution.
-null infile Read a null model from file. The default for
protein is to use average amino acid
frequencies from Swissprot 34 and p1 =
350/351; for nucleic acid, the default is to
use 0.25 for each base and p1 = 1000/1001.
For documentation of the format of the null
model file and further explanation of how
the null model is used, see the User's
Guide.
-pam infile Apply a heuristic PAM- (substitution
matrix-) based prior on match emission
probabilities instead of the default mixture
Dirichlet. The substitution matrix is read
from file. See -pamwgt. The default
Dirichlet state transition prior and insert
emission prior are unaffected. Therefore in
principle you could combine -prior with -pam
but this isn't recommended, as it hasn't
been tested. ( -pam itself hasn't been
tested much!)
-pamwgt float [20.0] Controls the weight
Standard (Mandatory) qualifiers
Allowed values
Default
[-alignfile]
(Parameter 1)(Aligned) protein sequence set filename and optional format, or reference (input USA)
Readable set of sequences
Required
-nhmm
Name for this HMM. The name can be any string of non-whitespace characters (e.g. one 'word'). There is no length limit (at least not one imposed by HMMER; your shell will complain about command line lengths frst).
Any string without whitespace is accepted
An empty string is accepted
-strategy
All alignments are local with respect to the sequence and are configured to be local (fragmentary) or global with respect to the HMM. The model is also configured to find a single or multiple domains (matches) to a sequence. The options for configuring the model are as follows: (D): The default setting. Multiple domains per sequence, global alignments with respect to the HMM. (F): Multiple domains per sequence, local alignments with respect to the HMM. Analogous to the old hmmfs program of HMMER 1. (G) Single domain per sequence, global alignment with respect to the HMM. Analogous to the old hmms program of HMMER 1. (S) Single domain per sequence, local alignments with respect to the HMM. Analogous to the old hmmsw program of HMMER 1.
D (global-multidomain) F (local-multidomain) G (global-singledomain) S (local-singledomain) D
[-hmmfile]
(Parameter 2)HMMER hidden markov model output file
Output file
<*>.ehmmbuild
Additional (Optional) qualifiers
Allowed values
Default
(none)
Advanced (Unprompted) qualifiers
Allowed values
Default
-prior
Read a Dirichlet prior from file, replacing the default mixture Dirichlet. The format of prior files is documented in the User's Guide, and an example is given in the Demos directory of the HMMER distribution.
Input file
Required
-null
Read a null model from file. The default for protein is to use average amino acid frequencies from Swissprot 34 and p1 = 350/351; for nucleic acid, the default is to use 0.25 for each base and p1 = 1000/1001. For documentation of the format of the null model file and further explanation of how the null model is used, see the User's Guide.
Input file
Required
-pam
Apply a heuristic PAM- (substitution matrix-) based prior on match emission probabilities instead of the default mixture Dirichlet. The substitution matrix is read from file. See -pamwgt. The default Dirichlet state transition prior and insert emission prior are unaffected. Therefore in principle you could combine -prior with -pam but this isn't recommended, as it hasn't been tested. ( -pam itself hasn't been tested much!)
Input file
Required
-pamwgt
Controls the weight <x> on a PAM-based prior. Only has effect if -pam option is also in use. <x> is a positive real number, 20.0 by default. <x> is the number of 'pseudocounts' contriubuted by the heuristic prior. Very high values of <x> can force a scoring system that is entirely driven by the substitution matrix, making HMMER somewhat approximate Gribskov profiles.
Any numeric value
20.0
-pbswitch
For alignments with a very large number of sequences, the GSC, BLOSUM, and Voronoi weighting schemes are slow; they're O(N^2) for N sequences. Henikoff position-based weights (PB weights) are more effcient. At or above a certain threshold sequence number <n> hmmbuild will switch from GSC, BLOSUM, or Voronoi weights to PB weights. To disable this switching behavior (at the cost of compute time, set <n> to be something larger than the number of sequences in your alignment. <n> is a positive integer; the default is 1000.
Any integer value
1000
-archpri
The value of the 'architecture prior' used by MAP architecture construction. This value is a probability between 0 and 1. This parameter governs a geometric prior distribution over model lengths. As 'archpri' increases, longer models are favored a priori. As 'archpri' decreases, it takes more residue conservation in a column to make a column a 'consensus' match column in the model architecture. The 0.85 default has been chosen empirically as a reasonable setting.
Any numeric value
0.85
-binary
Write the HMM to file in HMMER binary format instead of readable ASCII text.
Boolean value Yes/No
No
-fast
Quickly and heuristically determine the architecture of the model by assigning all columns with more than a certain fraction of gap characters to insert states. By default this fraction is 0.5, and it can be changed using the --gapmax option. The default construction algorithm is a maximum a posteriori (MAP) algorithm, which is slower.
Boolean value Yes/No
No
-gapmax
Controls the -fast model construction algorithm, but if -fast is not being used, has no effect. If a column has more than a fraction <x> of gap symbols in it, it gets assigned to an insert column. <x> is a frequency from 0 to 1, and by default is set to 0.5. Higher values of <x> mean more columns get assigned to consensus, and models get longer; smaller values of <x> mean fewer columns get assigned to consensus, and models get smaller.
Any numeric value
0.5
-hand
Specify the architecture of the model by hand: the alignment file must be in SELEX or Stockholm format, and the reference annotation line (RF in SELEX, GC RF in Stockholm) is used to specify the architecture. Any column marked with a non-gap symbol (such as an 'x', for instance) is assigned as a consensus (match) column in the model.
Boolean value Yes/No
No
-sidlevel
Controls both the determination of effective sequence number and the behavior of the -wblosum weighting option. The sequence alignment is clustered by percent identity, and the number of clusters at a cutoff threshold of <x> is used to determine the effective sequence number. Higher values of <x> give more clusters and higher effective sequence numbers; lower values of <x> give fewer clusters and lower effective sequence numbers. <x> is a fraction from 0 to 1, and by default is set to 0.62 (corresponding to the clustering level used in constructing the BLOSUM62 substitution matrix).
Any numeric value
0.62
-noeff
Turn off the effective sequence number calculation, and use the true number of sequences instead. This will usually reduce the sensitivity of the final model (so don't do it without good reason!)
Boolean value Yes/No
No
-swentry
Controls the total probability that is distributed to local entries into the model, versus starting at the beginning of the model as in a global alignment. <x> is a probability from 0 to 1, and by default is set to 0.5. Higher values of <x> mean that hits that are fragments on their left (N or 5'-terminal) side will be penalized less, but complete global alignments will be penalized more. Lower values of <x> mean that fragments on the left will be penalized more, and global alignments on this side will be favored. This option only affects the confgurations that allow local alignments, e.g. -s and -f; unless one of these options is also activated, this option has no effect. You have independent control over local/global alignment behavior for the N/C (5'/3') termini of your target sequences using --swentry and --swexit.
Any numeric value
0.5
-swexit
Controls the total probability that is distributed to local exits from the model, versus ending an alignment at the end of the model as in a global alignment. <x> is a probability from 0 to 1, and by default is set to 0.5. Higher values of <x> mean that hits that are fragments on their right (C or 3'-terminal) side will be penalized less, but complete global alignments will be penalized more. Lower values of <x> mean that fragments on the right will be penalized more, and global alignments on this side will be favored. This option only affects the confgurations that allow local alignments, e.g. -s and -f; unless one of these options is also activated, this option has no effect. You have independent control over local/global alignment behavior for the N/C (5'/3') termini of your target sequences using -swentry and -swexit.
Any numeric value
0.5
-verbosity
Print more possibly useful stuff, such as the individual scores for each sequence in the alignment.
Boolean value Yes/No
No
-weighting
Values (B)(-wblosum in HMMER) Use the BLOSUM filtering algorithm to weight the sequences. Cluster the sequences at a given percentage identity (see -idlevel); assign each cluster a total weight of 1.0, distributed equally amongst the members of that cluster. (G)(-wgsc in HMMER) Use the Gerstein/Sonnhammer/Chothia ad hoc sequence weighting algorithm. This is the default. (K)(-wme in HMMER) Use the Krogh/Mitchison maximum entropy algorithm to 'weight' the sequences. This supercedes the Eddy/Mitchison/Durbin maximum discrimination algorithm, which gives almost identical weights but is less robust. ME weighting seems to give a marginal increase in sensitivity over the default GSC weights, but takes a fair amount of time. (W) (-wpb in HMMER) Use the Henikoff position-based weighting scheme. (V) (-wvoronoi in HMMER) Use the Sibbald/Argos Voronoi sequence weighting algorithm in place of the default GSC weighting. (N) (-wnone in HMMER) Turn off all sequence weighting.
B (Blosum) G (Gerstein/Sonnhammer/Chothia) K (Krogh/Mitchison) W (Henikoff) V (Sibbald/Argos Voronoi) N (None) G
-o
Re-save the starting alignment to file, in Stockholm format. The columns which were assigned to match states will be marked with x's in an RF annotation line. If either the -hand or -fast construction options were chosen, the alignment may have been slightly altered to be compatible with Plan 7 transitions, so saving the final alignment and comparing to the starting alignment can let you view these alterations. See the User's Guide for more information on this arcane side effect.
Output file
<*>.ehmmbuild
-cfile
Save the observed emission and transition counts to file after the architecture has been determined (e.g. after residues/gaps have been assigned to match, delete, and insert states). This option is used in HMMER development for generating data files useful for training new Dirichlet priors. The format of count files is documented in the User's Guide.
Output file
<*>.ehmmbuild
Input file format
Alignment and sequence formats
Input and output of alignments and sequences is limited to the formats that the original hmmer supports. These include stockholm, SELEX, MSF, Clustal, Phylip and A2M /aligned FASTA (alignments) and FASTA, GENBANK, EMBL, GCG, PIR (sequences). It would be fairly straightforward to adapt the code to support all EMBOSS-supported formats.
Compressed input files
Automatic processing of gzipped files is not supported.
Input files for usage example
File: globins50.msf
!!AA_MULTIPLE_ALIGNMENT 1.0
PileUp of: *.pep
Symbol comparison table: GenRunData:blosum62.cmp CompCheck: 6430
GapWeight: 12
GapLengthWeight: 4
pileup.msf MSF: 308 Type: P August 16, 1999 09:09 Check: 9858 ..
Name: lgb1_pea Len: 308 Check: 2200 Weight: 1.00
Name: lgb1_vicfa Len: 308 Check: 214 Weight: 1.00
Name: myg_escgi Len: 308 Check: 3961 Weight: 1.00
Name: myg_horse Len: 308 Check: 5619 Weight: 1.00
Name: myg_progu Len: 308 Check: 6401 Weight: 1.00
Name: myg_saisc Len: 308 Check: 6606 Weight: 1.00
Name: myg_lycpi Len: 308 Check: 6090 Weight: 1.00
Name: myg_mouse Len: 308 Check: 6613 Weight: 1.00
Name: myg_musan Len: 308 Check: 3942 Weight: 1.00
Name: hba_ailme Len: 308 Check: 4558 Weight: 1.00
Name: hba_prolo Len: 308 Check: 5054 Weight: 1.00
Name: hba_pagla Len: 308 Check: 5383 Weight: 1.00
Name: hba_macfa Len: 308 Check: 5135 Weight: 1.00
Name: hba_macsi Len: 308 Check: 5198 Weight: 1.00
Name: hba_ponpy Len: 308 Check: 5050 Weight: 1.00
Name: hba2_galcr Len: 308 Check: 5609 Weight: 1.00
Name: hba_mesau Len: 308 Check: 4702 Weight: 1.00
Name: hba2_bosmu Len: 308 Check: 4241 Weight: 1.00
Name: hba_erieu Len: 308 Check: 4680 Weight: 1.00
Name: hba_frapo Len: 308 Check: 3549 Weight: 1.00
Name: hba_phaco Len: 308 Check: 4440 Weight: 1.00
Name: hba_trioc Len: 308 Check: 5465 Weight: 1.00
Name: hba_ansse Len: 308 Check: 3300 Weight: 1.00
Name: hba_colli Len: 308 Check: 3816 Weight: 1.00
Name: hbad_chlme Len: 308 Check: 4571 Weight: 1.00
Name: hbad_pasmo Len: 308 Check: 6777 Weight: 1.00
Name: hbaz_horse Len: 308 Check: 7187 Weight: 1.00
Name: hba4_salir Len: 308 Check: 7329 Weight: 1.00
Name: hbb_ornan Len: 308 Check: 2667 Weight: 1.00
Name: hbb_tacac Len: 308 Check: 4356 Weight: 1.00
Name: hbe_ponpy Len: 308 Check: 3827 Weight: 1.00
Name: hbb_speci Len: 308 Check: 1556 Weight: 1.00
Name: hbb_speto Len: 308 Check: 2051 Weight: 1.00
Name: hbb_equhe Len: 308 Check: 3414 Weight: 1.00
Name: hbb_sunmu Len: 308 Check: 2927 Weight: 1.00
Name: hbb_calar Len: 308 Check: 3836 Weight: 1.00
Name: hbb_mansp Len: 308 Check: 4322 Weight: 1.00
Name: hbb_ursma Len: 308 Check: 4428 Weight: 1.00
Name: hbb_rabit Len: 308 Check: 4190 Weight: 1.00
Name: hbb_tupgl Len: 308 Check: 4185 Weight: 1.00
[Part of this file has been deleted for brevity]
lgb1_pea ~~~~~~~~
lgb1_vicfa ~~~~~~~~
myg_escgi ~~~~~~~~
myg_horse ~~~~~~~~
myg_progu ~~~~~~~~
myg_saisc ~~~~~~~~
myg_lycpi ~~~~~~~~
myg_mouse ~~~~~~~~
myg_musan ~~~~~~~~
hba_ailme ~~~~~~~~
hba_prolo ~~~~~~~~
hba_pagla ~~~~~~~~
hba_macfa ~~~~~~~~
hba_macsi ~~~~~~~~
hba_ponpy ~~~~~~~~
hba2_galcr ~~~~~~~~
hba_mesau ~~~~~~~~
hba2_bosmu ~~~~~~~~
hba_erieu ~~~~~~~~
hba_frapo ~~~~~~~~
hba_phaco ~~~~~~~~
hba_trioc ~~~~~~~~
hba_ansse ~~~~~~~~
hba_colli ~~~~~~~~
hbad_chlme ~~~~~~~~
hbad_pasmo ~~~~~~~~
hbaz_horse ~~~~~~~~
hba4_salir ~~~~~~~~
hbb_ornan ~~~~~~~~
hbb_tacac ~~~~~~~~
hbe_ponpy ~~~~~~~~
hbb_speci ~~~~~~~~
hbb_speto ~~~~~~~~
hbb_equhe ~~~~~~~~
hbb_sunmu ~~~~~~~~
hbb_calar ~~~~~~~~
hbb_mansp ~~~~~~~~
hbb_ursma ~~~~~~~~
hbb_rabit ~~~~~~~~
hbb_tupgl ~~~~~~~~
hbb_triin ~~~~~~~~
hbb_colli ~~~~~~~~
hbb_larri ~~~~~~~~
hbb1_varex ~~~~~~~~
hbb2_xentr ~~~~~~~~
hbbl_ranca ~~~~~~~~
hbb2_tricr ~~~~~~~~
glb2_mormr ~~~~~~~~
glbz_chith FGAVFAKM
hbf1_ureca VAAMK~~~
Input files for usage example 2
File: rrm.sto
# STOCKHOLM 1.0
YIS5_YEAST/33-104 IYIGNL......NRELTEGDILTVFS.....E.YGVP..VDVILSRD...
YIS1_YEAST/66-136 IFVGNI......TPDVTPEQIEDHFK.....D.CGQI..KRITLLYD...
RN15_YEAST/20-91 VYLGSI......PYDQTEEQILDLCS.....N.VGPV..INLKMMFD...
CST2_HUMAN/18-89 VFVGNI......PYEATEEQLKDIFS.....E.VGPV..VSFRLVYD...
NSR1_YEAST/269-340 LFLGNL......SFNADRDAIFELFA.....K.HGEV..VSVRIPTH...
IF4B_HUMAN/98-168 AFLGNL......PYDVTEESIKEFFR.....G.LNIS...AVRLPR....
TIA1_HUMAN/108-179 VFVGDL......SPQITTEDIKAAFA.....P.FGRI..SDARVVKD...
PUB1_YEAST/163-234 LFVGDL......NVNVDDETLRNAFK.....D.FPSY..LSGHVMWD...
NAM8_YEAST/165-237 IFVGDL......APNVTESQLFELFI.....NRYAST..SHAKIVHD...
RT19_ARATH/33-104 LYIGGL......SPGTDEHSLKDAFS.....S.FNGV..TEARVMTN...
NSR1_YEAST/170-241 IFVGRL......SWSIDDEWLKKEFE.....H.IGGV..IGARVIYE...
PES4_YEAST/93-164 LFIGDL......HETVTEETLKGIFK.....K.YPSF..VSAKVCLD...
PABP_DROME/4-75 LYVGDL......PQDVNESGLFDKFS.....S.AGPV..LSIRVCRD...
NOP4_YEAST/28-98 LFVRSI......PQDVTDEQLADFFS.....N.FAPI..KHAVVVKD...
ROC5_NICSY/219-290 LYVANL......SWALTSQGLRDAFA.....D.QPGF..MSAKVIYD...
ROC5_NICSY/116-187 LYVGNL......PFSMTSSQLSEIFA.....E.AGTV..ANVEIVYD...
ROC3_NICSY/99-170 LFVGNL......PYDIDSEGLAQLFQ.....Q.AGVV..EIAEVIYN...
GBP2_YEAST/221-291 VFIINL......PYSMNWQSLKDMFK.....E.CGHV..LRADVELD...
PABP_SCHPO/263-333 VYIKNL......DTEITEQEFSDLFG.....Q.FGEI..TSLSLVKD...
PABP_DROME/183-254 VYVKNF......TEDFDDEKLKEFFE.....P.YGKI..TSYKVMS....
PABP_DROME/92-162 VFIKNL......DRAIDNKAIYDTFS.....A.FGNI..LSCKVATD...
SFR3_HUMAN/12-78 VYVGNL......GNNGNKTELERAFG.....Y.YGPL..RSVWVARN...
SR55_DROME/5-68 VYVGGL......PYGVRERDLERFFK.....G.YGRT..RDILIKN....
SFR1_HUMAN/17-85 IYVGNL......PPDIRTKDIEDVFY.....K.YGAI..RDIDLKNR...
SFR2_CHICK/16-87 LKVDNL......TYRTSPDTLRRVFE.....K.YGRV..GDVYIPRD...
ROG_HUMAN/10-81 LFIGGL......NTETNEKALEAVFG.....K.YGRI..VEVLLMKD...
TRA2_DROME/99-170 IGVFGL......NTNTSQHKVRELFN.....K.YGPI..ERIQMVID...
ROC_HUMAN/18-82 VFIGNL.....NTLVVKKSDVEAIFS.....K.YGKI..VGCSVHK....
NONA_DROME/304-369 LYVGNL......TNDITDDELREMFK.....P.YGEI..SEIFSNLD...
RU17_DROME/104-175 LFIARI......NYDTSESKLRREFE.....F.YGPI..KKIVLIHD...
TIA1_HUMAN/216-281 VYCGGV......TSGLTEQLMRQTFS.....P.FGQI..MEIRVFPD...
TIA1_HUMAN/9-78 LYVGNL......SRDVTEALILQLFS.....Q.IGPC..KNCKMIMD...
PUB1_YEAST/76-146 LYVGNL......DKAITEDILKQYFQ.....V.GGPI..ANIKIMID...
SXL_DROME/213-285 LYVTNL......PRTITDDQLDTIFG.....K.YGSI..VQKNILRD...
ELAV_DROME/250-322 LYVSGL......PKTMTQQELEAIFA.....P.FGAI..ITSRILQN...
SXL_DROME/127-198 LIVNYL......PQDMTDRELYALFR.....A.IGPI..NTCRIMRD...
ELV4_HUMAN/48-119 LIVNYL......PQNMTQEEFRSLFG.....S.IGEI..ESCKLVRD...
MSSP_HUMAN/31-102 LYIRGL......PPHTTDQDLVKLCQ.....P.YGKI..VSTKAILD...
ELAV_DROME/404-475 IFIYNL......APETEEAALWQLFG.....P.FGAV..QSVKIVKD...
SQD_DROME/138-208 IFVGGL......TTEISDEEIKTYFG.....Q.FGNI..VEVEMPLD...
ROA1_BOVIN/106-176 IFVGGI......KEDTEEHHLRDYFE.....Q.YGKI..EVIEIMTD...
CABA_MOUSE/161-231 IFVGGL......NPEATEEKIREYFG.....Q.FGEI..EAIELPID...
SQD_DROME/58-128 LFVGGL......SWETTEKELRDHFG.....K.YGEI..ESINVKTD...
RB97_DROME/34-104 LFIGGL......APYTTEENLKLFYG.....Q.WGKV..VDVVVMRD...
GR10_BRANA/8-79 CFVGGL......AWATGDAELERTFS.....Q.FGEV..IDSKIIND...
CABA_MOUSE/77-147 MFVGGL......SWDTSKKDLKDYFT.....K.FGEV..VDCTIKMD...
[Part of this file has been deleted for brevity]
ROA1_BOVIN/106-176 NCEV
CABA_MOUSE/161-231 KCEI
SQD_DROME/58-128 KVDP
RB97_DROME/34-104 TVEA
GR10_BRANA/8-79 TITV
CABA_MOUSE/77-147 VIDP
SFPQ_HUMAN/373-443 TTPR
NONA_DROME/378-448 ASLR
NOP3_YEAST/127-190 PLEV
GBP2_YEAST/124-193 KLMV
U2AF_SCHPO/312-383 KLHA
U2AF_HUMAN/261-332 KLLV
NUCL_CHICK/283-352 EIKL
MODU_DROME/342-410 EWKL
RU2B_HUMAN/9-81 PMRI
RU1A_HUMAN/12-84 PMRI
PUB1_YEAST/342-407 NLRT
PTB_HUMAN/186-253 CCTL
EWS_HUMAN/363-442 KLKV
ARP2_PLAFA/364-438 YLKV
YHH5_YEAST/315-384 TIWV
PES4_YEAST/305-374 KLLV
YHC4_YEAST/348-415 TMHV
PR24_YEAST/212-284 EISV
MODU_DROME/422-484 PIFI
RNP1_YEAST/37-109 KIVI
GBP2_YEAST/351-421 SLQI
NOP4_YEAST/292-363 LIGD
MEI2_SCHPO/197-265 LLYF
IF39_YEAST/79-157 RLFL
ROF_HUMAN/113-183 YIEV
MODU_DROME/260-326 FITV
WHI3_YEAST/540-614 TVSS
CPO_DROME/453-526 PQTI
RU2B_HUMAN/153-220 AMKI
RU1A_HUMAN/210-276 NAMK
RU1A_YEAST/229-293 NDVT
SFR1_HUMAN/122-186 EGET
NOP3_YEAST/202-270 VITV
U2AG_HUMAN/67-142 PIHA
D111_ARATH/285-364 TVRA
PR24_YEAST/43-111 EIIV
MODU_DROME/177-246 VLSV
SSB1_YEAST/39-114 EIHI
RN12_YEAST/200-267 VLHI
PTB_HUMAN/61-128 PIYI
LA_HUMAN/113-182 DLLI
LA_DROME/151-225 ELLR
PR24_YEAST/119-190 TLVT
IF39_SCHPO/41-124 TFVV
//
Input files for usage example 3
File: fn3.sto
# STOCKHOLM 1.0
#=GF ID fn3
#=GF AC PF00041
#=GF DE Fibronectin type III domain
#=GF AU Sonnhammer ELL
#=GF GA 7.8 -1.0
#=GF TC 7.7 7.7
#=GF TC 7.9 0.1
#=GF AL HMM_simulated_annealing
#=GF SE Swissprot_feature_table
#=GF TP Domain
#=GF BM hmmbuild -F HMM_ls SEED
#=GF BM hmmcalibrate --seed 0 HMM_ls
#=GF BM hmmbuild -f -F HMM_fs SEED
#=GF BM hmmcalibrate --seed 0 HMM_fs
#=GF AM globalfirst
#=GF RN [1]
#=GF RM 90384924
#=GF RT Structural design and molecular evolution of a cytokine
#=GF RT receptor superfamily.
#=GF RA Bazan JF;
#=GF RL PNAS USA 1990;87:6934-6938.
#=GF RN [2]
#=GF RM 95106303
#=GF RT Tracing the spread of fibronectin type III domains in
#=GF RT bacterial glycohydrolases.
#=GF RA Little E, Bork P, Doolittle R;
#=GF RL J Mol Evol 1994;39:631-643.
#=GF RC Definition of fibronectin domains
#=GF RN [3]
#=GF RM 85284965
#=GF RT Primary structure of human fibronectin: differential
#=GF RT splicing may generate at least 10 polypeptides from a
#=GF RT single gene.
#=GF RA Kornblihtt AR, Umezawa K, Vibe-Pedersen K, Baralle FE;
#=GF RL EMBO J. 1985;4:1755-1759.
#=GF DR PROSITE; PDOC00214;
#=GF DR PRINTS; PR00014;
#=GF DR SCOP; 1ttf; sf;
#=GF DR SMART; FN3;
#=GF DR HOMSTRAD; fn3;
#=GF DC The following Pfam-B families contain sequences that according to Prodom
#=GF DC are members of this Pfam-A family.
#=GF DR PFAMB; PB001437;
#=GF DR PFAMB; PB004758;
#=GF DR PFAMB; PB007041;
#=GF DR PFAMB; PB007646;
#=GF DR PFAMB; PB012739;
#=GF DR PFAMB; PB016725;
#=GF DR PFAMB; PB016744;
[Part of this file has been deleted for brevity]
FAS2_SCHAM/530-616 RFA.AQNEVGF....GPWS
PTPZ_HUMAN/312-401 QIV.AICTNGL...YGKYS
NCA1_BOVIN/610-691 YVV.AENQQ.......GKS
FAS2_SCHAM/642-735 EVR.ATNAIGN....SVPG
TIE2_HUMAN/543-626 RAR..VNTKAQ....GEWS
TENA_CHICK/682-767 SLH.IVKNNTR...GPGLS
UFO_HUMAN/327-411 CVA.AYTAAGD....GPWS
TIE1_HUMAN/545-632 DVQ.LYHCTLL....GPAS
7LES_DROME/1799-1891 WVQ.AHATPTK....SNSS
NRG_DROME/717-799 KVV.AINDR.......GES
#=GR NRG_DROME/717-799 SS EEE.EEEEE.......EEE
#=GR NRG_DROME/717-799 SA 030.10200.......020
NRCA_CHICK/726-810 KVQ.ALNDLGY...APEPS
NGCA_CHICK/700-794 RVQ.AVNGAGK....GPEA
TIE2_HUMAN/639-724 DIF.AENNIGS....SNPA
LAR_DROME/515-598 WLA.ARSQRGE....GATT
PTP6_DROME/236-321 RIV.GKNSIGN....GQPT
TIE2_HUMAN/444-529 CVQ.LVRRGEG....GEGH
TIE1_HUMAN/446-533 RVQLSRPGEGG...EGAWG
PTPK_MOUSE/290-376 RVLLTRPGEGG...TGLPG
PTPB_HUMAN/1355-1434 SIK.VQSAG.......MTS
PTP1_DROME/583-661 QVN.TVSFG.......VES
PTP1_DROME/312-394 IVK.TVSGK.......VTS
PTP1_DROME/217-301 SVQ.TMSED.......EIS
PTP1_DROME/123-205 QAY.TIYDG.......KES
PTPB_HUMAN/1261-1344 NVK.TVSGDSW....KTYS
PTP1_DROME/864-944 TVV.VRSGTESS..VLRSS
PTPB_HUMAN/1173-1250 WVV.THSGD........LS
PTPB_HUMAN/1085-1162 VIV.THSGE........LS
PTPB_HUMAN/467-543 TVT.SISGD........LK
PTPB_HUMAN/995-1074 QIL.TVSGG.......LFS
PTPB_HUMAN/731-808 TVT.TKSGQ........YE
PTPB_HUMAN/643-725 TIT.TRSGKYE...NHSFS
PTPB_HUMAN/554-632 VVT.TVSGG.......ISS
PTPB_HUMAN/907-984 MIA.SVSGS........LK
PTPB_HUMAN/112-192 AIT.AVSGG.......KRS
IL7R_HUMAN/129-221 KVR.SIPDHYFKGFWSEWS
FINC_BOVIN/689-768 NVY.EISEE.......GEQ
KALM_CHICK/178-271 RVA.AVNVHGT...RGFTA
TENA_HUMAN/1529-1608 MVS.GFTQG.......HQT
TENA_CHICK/1227-1306 ELY.GVSSG.......RRS
TENA_CHICK/1136-1215 YLY.GISHG.......FRT
TENA_HUMAN/1255-1335 TLH.GEVRG.......HST
TENA_CHICK/1045-1124 TLY.GVIRG.......YRT
7LES_DROVI/1917-1997 RLA.LAYAATP....GAPI
PTPB_HUMAN/22-103 KII.SLDEE........RT
KALM_CHICK/544-642 EVQ.VLTTGGE....GPAT
#=GC SS_cons EEE.EEECCCC....CCCC
#=GC SA_cons 221.2123454....3232
#=GC seq_cons pVt.Ahsst.......t.S
//
Input files for usage example 4
File: pkinase.sto
# STOCKHOLM 1.0
#=GF ID pkinase
#=GF AC PF00069
#=GF DE Protein kinase domain
#=GF AU Sonnhammer ELL
#=GF GA -57.9 -57.9
#=GF TC -58.0 -58.0
#=GF TC -57.7 -57.7
#=GF AL Clustalw
#=GF SE Unknown
#=GF TP Domain
#=GF BM hmmbuild -F HMM_ls SEED
#=GF BM hmmcalibrate --seed 0 HMM_ls
#=GF BM hmmbuild -f -F HMM_fs SEED
#=GF BM hmmcalibrate --seed 0 HMM_fs
#=GF AM byscore
#=GF RN [1]
#=GF RM 92065863
#=GF RT Protein kinase catalytic domain sequence database:
#=GF RT identification of conserved features of primary structure
#=GF RT and classification of family members.
#=GF RA Hanks SK, Quinn AM;
#=GF RL Methods Enzymol 1991;200:38-62.
#=GF RN [2]
#=GF RM 95285959
#=GF RT Protein kinases 6. The eukaryotic protein kinase
#=GF RT superfamily: kinase (catalytic) domain structure and
#=GF RT classification.
#=GF RA Hanks SK, Hunter T;
#=GF RL FASEB J 1995;9:576-596.
#=GF RN [3]
#=GF RM 97172697
#=GF RT The protein kinases of budding yeast: six score and more.
#=GF RA Hunter T, Plowman GD;
#=GF RL Trends Biochem Sci 1997;22:18-22.
#=GF DR PROSITE; PDOC00100;
#=GF DR PROSITE; PDOC00212;
#=GF DR PROSITE; PDOC00213;
#=GF DR PROSITE; PDOC00629;
#=GF DR PROSITE_PROFILE; PS50011;
#=GF DR PRINTS; PR00109;
#=GF DC Protein kinase; unclassified specificity.
#=GF DR SMART; STYKc;
#=GF DC Serine/Threonine protein kinases, catalytic domain
#=GF DR SMART; S_TKc;
#=GF DC Tyrosine kinase, catalytic domain
#=GF DR SMART; TyrKc;
#=GF DR SCOP; 1apm; sf;
#=GF DR URL; http://www.sdsc.edu/Kinases/pk_home.html;
#=GF DR HOMSTRAD; kinase;
[Part of this file has been deleted for brevity]
#=GR SAT4_YEAST/316-590 AS .....................
KKL6_YEAST/192-508 TLTRKN.QRWSINEIYESPFV
#=GR KKL6_YEAST/192-508 AS .....................
RAN1_SCHPO/18-295 PKTRI.....TLPEL..STLV
#=GR RAN1_SCHPO/18-295 AS .....................
TRKA_HUMAN/510-781 PQQRH.....SIKDVHARL..
#=GR TRKA_HUMAN/510-781 AS .....................
KSYK_PIG/364-619 VENRP.....GFVAVELRL..
#=GR KSYK_PIG/364-619 AS .....................
BFR2_HUMAN/367-643 PSQRP.....TFKQLVEDL..
#=GR BFR2_HUMAN/367-643 AS .....................
ABL1_CAEEL/311-562 PSDRP.....RFRDIHFNL..
#=GR ABL1_CAEEL/311-562 AS .....................
7LES_DROME/2209-2481 PWERP.....SFRRCYNTL..
JAK1_HUMAN/864-1137 PSNRT.....SFQNLIEGF..
#=GR JAK1_HUMAN/864-1137 AS .....................
EGFR_DROME/938-1194 AAMRP.....TFKQLTTVF..
#=GR EGFR_DROME/938-1194 AS .....................
RYK_HUMAN/327-593 PEERP.....KFQQLVQCL..
#=GR RYK_HUMAN/327-593 AS .....................
MET_HUMAN/1078-1337 AEMRP.....SFSELVSRI..
#=GR MET_HUMAN/1078-1337 AS .....................
PHY1_CERPU/1004-1282 PLKRP......SFAVICQKL.
#=GR PHY1_CERPU/1004-1282 AS .....................
M3K9_HUMAN/3-262 PHSRP.....SFTNILDQL..
#=GR M3K9_HUMAN/3-262 AS .....................
KYK2_DICDI/108-364 PNDRP.....PFTYIVNKL..
#=GR KYK2_DICDI/108-364 AS .....................
KYK1_DICDI/1289-1559 PDSRP.....TFKQIIVHL..
#=GR KYK1_DICDI/1289-1559 AS .....................
KMIL_AVIMH/82-339 REERP.....LFPQILSSI..
#=GR KMIL_AVIMH/82-339 AS .....................
KPRO_MAIZE/534-812 RSKRP.....TMEHAVQTL..
#=GR KPRO_MAIZE/534-812 AS .....................
AVR1_HUMAN/208-495 PSARL.....TALRIKKTL..
#=GR AVR1_HUMAN/208-495 AS .....................
AVR2_HUMAN/192-479 AEARL.....SAGCVGERI..
#=GR AVR2_HUMAN/192-479 AS .....................
PKN1_MYXXA/59-320 PEDRY....ASIAAFRNALQV
#=GR PKN1_MYXXA/59-320 AS .....................
KMOS_CERAE/60-338 AAQRP.....SARPLLVDL..
#=GR KMOS_CERAE/60-338 AS .....................
KR1_PRVKA/53-332 PVRRP.....SADEILNFGMW
#=GR KR1_PRVKA/53-332 AS .....................
KR1_HSV11/191-478 GALRP.....SAAELLCLPLF
#=GR KR1_HSV11/191-478 AS .....................
JAK1_HUMAN/571-833 PNQRP.....FFRAIMRDI..
HR25_YEAST/9-273 FDEKP......DYLFLARLF.
#=GR HR25_YEAST/9-273 AS .....................
#=GC seq_cons PppR......ohppllpc.ah
//
Output file format
ehmmbuild
outputs a graph to the specified graphics device.
outputs a report format file. The default format is ...
Output files for usage example
File: globin.hmm
HMMER2.0 [2.3.2]
NAME globins50
LENG 143
ALPH Amino
RF no
CS no
MAP yes
COM hmmbuild -n globins50 --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F globin.hmm ../../data/hmmnew/globins50.msf
NSEQ 50
DATE Tue Jul 15 12:00:00 2008
CKSUM 9858
XT -8455 -4 -1000 -1000 -8455 -4 -8455 -4
NULT -4 -8455
NULE 595 -1558 85 338 -294 453 -1158 197 249 902 -1085 -142 -21 -313 45 531 201 384 -1998 -644
HMM A C D E F G H I K L M N P Q R S T V W Y
m->m m->i m->d i->m i->i d->m d->d b->m m->e
-450 * -1900
1 591 -1587 159 1351 -1874 -201 151 -1600 998 -1591 -693 389 -1272 595 42 -31 27 -693 -1797 -1134 14
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 -450 *
2 -926 -2616 2221 2269 -2845 -1178 -325 -2678 -300 -2596 -1810 220 -1592 939 -974 -671 -939 -2204 -2785 -1925 15
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
3 -638 -1715 -680 497 -2043 -1540 23 -1671 2380 -1641 -840 -222 -1595 437 1040 -564 -523 -1363 2124 -1313 16
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
4 829 -1571 -37 660 -1856 -873 152 -1578 894 -1573 -678 769 -1273 1284 58 224 447 -1175 -1782 -1125 17
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
5 369 -433 -475 286 -974 -1312 -19 -412 664 398 406 1030 -1394 388 -214 -261 85 -166 -1227 -725 18
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
6 -1291 -884 -3696 -3261 -1137 -3425 -2802 2322 -3066 111 19 -3028 -3275 -2855 -3100 -2670 -1269 2738 -2450 -2062 19
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
7 157 -413 -236 316 -1387 -1231 89 -863 1084 -431 -348 910 -1319 635 297 15 704 -483 -1497 -922 20
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
8 770 -1431 -43 459 -1751 -340 78 -1449 440 -1497 -631 866 -1302 825 -51 953 364 -1076 -1750 -1121 21
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
9 420 -186 -2172 -1577 8 -1818 -694 1477 -1281 760 614 -1299 -1867 -1001 -1262 -189 -12 1401 -722 -364 22
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
10 -961 -879 -2277 -1821 1366 -2213 -204 -399 -1500 -130 -39 -1427 -2266 -1186 -1511 -159 -913 -367 4721 1177 23
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
11 -48 -1782 809 844 -2073 1456 8 -1811 315 -1803 -932 180 -1365 921 -218 173 -115 -1399 -2018 -1327 24
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -68 -6528 -4832 -894 -1115 -701 -1378 * *
[Part of this file has been deleted for brevity]
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
128 -415 -1926 1575 1399 -2219 -1163 17 -1983 527 -1929 -1039 341 -1367 1597 -212 257 -222 -1536 -2109 -1387 144
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
129 -529 -1434 -629 -143 -1926 -626 -171 -1460 2679 -1597 -839 -309 -1599 207 317 -530 -510 -130 -1840 -1369 145
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
130 811 -397 -2389 -1807 1883 -2039 -907 594 -1512 1077 687 -1532 -2065 -1201 -1483 -1125 -465 1067 -843 -472 146
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
131 -241 -102 -2327 -1710 724 -1767 -616 650 -1363 1074 1765 -718 -1809 -1026 -1252 -842 -181 1331 -541 695 147
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
132 723 95 385 823 -1820 -1168 167 -1540 875 -1362 -644 320 -1261 810 246 693 -67 -1141 -1753 -1098 148
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
133 551 -430 -1049 -481 -442 469 -241 465 -313 133 947 -411 -1543 197 -587 -146 202 522 -843 -429 149
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
134 -1086 -777 -3351 -2800 816 -2898 -1861 1501 -2515 1149 586 -2483 -2775 -2108 -2400 -2046 -1030 2380 -1511 -1216 150
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
135 1393 1409 -876 -345 -997 -525 -315 -590 -198 -847 -109 -420 -1441 -97 412 766 -130 139 -1306 -858 151
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
136 98 -1299 36 365 -1495 -1211 1241 -404 523 -952 -426 1174 -1303 511 -18 347 882 -853 -1566 -970 152
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
137 1308 -787 564 -132 -966 -1332 -203 -362 -49 -395 -57 -305 -1481 49 -437 -190 -182 1020 -1282 -802 153
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
138 -1746 -1358 -3897 -3341 -216 -3621 -2478 1774 -3040 2442 1158 -3189 -3229 -2422 -2853 -2824 -1659 392 -1720 -1647 154
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
139 1176 -1289 -179 534 -1606 -607 34 -1278 734 -1372 -534 44 -1325 433 -89 521 826 -941 -1666 -1072 155
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6528 -7571 -894 -1115 -701 -1378 * *
140 602 -1500 -135 850 -1753 -1214 1951 -1452 838 -1484 431 118 -1306 555 347 489 -153 -1085 -1723 -1092 156
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -22 -6602 -7644 -894 -1115 -701 -1378 * *
141 351 -1646 -165 546 -1976 -498 46 -1667 2193 -1662 -798 35 -1405 476 311 -73 -306 -1287 -1859 -1254 157
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -6561 -7603 -894 -1115 -701 -1378 * *
142 -1995 -1606 -3095 -2870 1739 -3015 -98 -1012 -2520 -730 655 -1990 -2962 -1884 -2326 -2167 -1915 -1128 548 4089 158
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -25 -6455 -7497 -894 -1115 -701 -1378 * *
143 -253 -1373 -267 301 -911 -565 1956 -450 1188 -1330 -497 33 -1352 502 1358 -205 -184 -941 -1604 -1026 159
- * * * * * * * * * * * * * * * * * * * *
- * * * * * * * * 0
//
Output files for usage example 2
File: myhmms
HMMER2.0 [2.3.2]
NAME rrm
LENG 77
ALPH Amino
RF no
CS no
MAP yes
COM hmmbuild -n rrm --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F myhmms ../../data/hmmnew/rrm.sto
NSEQ 90
DATE Tue Jul 15 12:00:00 2008
CKSUM 8325
XT -8455 -4 -1000 -1000 -8455 -4 -8455 -4
NULT -4 -8455
NULE 595 -1558 85 338 -294 453 -1158 197 249 902 -1085 -142 -21 -313 45 531 201 384 -1998 -644
HMM A C D E F G H I K L M N P Q R S T V W Y
m->m m->i m->d i->m i->i d->m d->d b->m m->e
-16 * -6492
1 -1084 390 -8597 -8255 -5793 -8424 -8268 2395 -8202 2081 -1197 -8080 -8115 -8020 -8297 -7789 -5911 1827 -7525 -7140 1
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 -16 *
2 -2140 -3785 -6293 -2251 3226 -2495 -727 -638 -2421 -545 -675 -5146 -5554 -4879 -1183 -2536 -1928 267 76 3171 2
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
3 -2542 458 -8584 -8273 -6055 -8452 -8531 2304 -8255 -324 101 -8104 -8170 -8221 -8440 -7840 -5878 3145 -7857 -7333 3
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
4 -1505 -5144 -1922 -558 -1842 2472 -3303 -2213 1099 -5160 -4233 372 -4738 -530 1147 168 498 -4766 -5327 -1476 4
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
5 -3724 -5184 300 -3013 -1655 1803 -3353 -5245 -1569 -2686 -4276 3495 -1963 -1331 -1054 -1472 -3664 -4803 -5369 -2 5
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
6 -1569 -6106 -8967 -8363 555 -8531 -7279 654 -8092 2953 -94 -8220 -7908 -1643 -7682 -7771 -6460 -59 -6191 -6284 6
- -151 -504 230 45 -380 399 101 -621 211 -470 -713 278 399 48 91 360 113 -364 -299 -254
- -178 -3113 -12684 -1600 -578 -701 -1378 * *
7 -409 -5130 -215 -2987 -1709 -956 690 -5188 -395 -5144 -4224 729 3054 -2862 -3409 354 1293 -1381 -5321 -4644 13
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
8 -3674 -5118 -1004 639 420 -4652 176 -2050 404 -1039 -935 16 1755 168 147 -275 198 -1472 1889 1977 14
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
9 -408 -5134 2415 1299 -950 -66 -767 -1296 -2889 -1843 -4224 1084 -968 -1439 -1854 540 -314 -2304 -5320 -60 15
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
10 586 1804 -6294 -631 -1627 -1671 -4374 1029 -2223 -162 1172 -5147 -5554 -1870 -5058 -2327 1741 1687 -4242 687 16
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
11 -2134 -5144 845 -1187 -1652 -1667 -3303 -5216 -513 -801 -4233 1026 -1873 -543 -619 575 2956 -4766 -5327 -4644 17
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
[Part of this file has been deleted for brevity]
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
62 -1537 -3785 -6300 -5664 -786 -5505 -4376 2171 652 1595 833 -5150 -5555 -869 -1636 -2574 -3898 1605 -4242 -1298 85
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
63 57 -5144 1387 1323 -5465 -2264 -3303 -780 812 -1898 -4233 131 -4738 1747 1061 108 -1766 -907 -5327 -793 86
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -310 -11642 -2372 -894 -1115 -701 -1378 * *
64 1233 705 139 1044 -1390 -413 -504 -4926 814 -1849 -873 617 -4450 473 -254 605 234 -2178 137 -4356 87
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11332 -12374 -894 -1115 -123 -3612 * *
65 -375 873 -6218 -1441 516 -5480 -4348 -3270 -294 2247 2716 -5104 -5531 -342 -2306 -4562 -1171 -333 -4233 564 88
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11626 -12668 -894 -1115 -382 -2105 * *
66 -2168 -5145 1778 -1715 -5466 -1556 1110 -5217 -347 -5161 -4234 3166 -102 1100 -928 -613 -1766 -4767 -5328 -411 89
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
67 -1458 -5144 -1823 270 -984 2973 159 -5216 311 -5160 -4233 841 -4738 -1432 -954 -817 -981 -4766 -5328 -609 90
- -148 -501 232 42 -379 397 104 -627 216 -467 -722 274 395 48 95 358 116 -365 -296 -251
- -157 -3284 -12684 -21 -6140 -701 -1378 * *
68 -674 -5080 -1101 -628 453 -1645 1133 200 1615 -1341 445 -3310 -4756 953 1117 -19 662 -1560 -5281 955 92
- -149 -500 232 43 -381 398 105 -627 210 -466 -721 275 396 45 95 359 120 -366 -295 -250
- -20 -6226 -12684 -1519 -619 -701 -1378 * *
69 -3696 -785 -17 1147 150 -1499 -3377 543 1044 -698 461 -736 -351 445 -584 -975 753 89 793 475 96
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
70 -3956 -3787 -1147 -2465 2321 -2554 -4371 1615 -2185 1495 -817 -5139 -2093 -4870 -2110 -2380 -1415 1282 -4244 800 97
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
71 -1015 -904 926 188 -880 1152 -650 -5215 677 -786 923 1723 -676 247 421 -1028 -1335 -4766 -5327 275 98
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
72 -2168 -5150 604 -1054 -5471 2761 -702 -5222 -918 -2141 -4239 1981 -4742 -142 -1740 262 -2027 -4772 -5333 -4650 99
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
73 -1412 853 -3523 -1169 -883 -4647 1925 -5199 1786 -1944 -937 588 -4740 633 2507 -312 -306 -4755 -5321 -528 100
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
74 -301 -498 -1749 1022 -866 -4653 -3313 -1356 1140 -1447 -4209 233 1336 -200 701 -1399 1158 451 -5307 -378 101
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
75 -1685 926 -1227 -5667 -1443 -2425 -4376 2317 -5261 1910 668 -5151 -2133 -1673 -5061 -2536 -2110 1189 -5 -3899 102
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
76 -2255 -5128 -302 363 -784 -2353 1398 -2058 1512 -314 -275 -1693 -1200 -571 1806 -524 252 -175 -255 793 103
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
77 -633 879 -2198 -5620 -1457 -5498 -4367 1846 -2705 -311 -2992 -5128 -1361 -1689 -695 -1593 53 2770 -4246 -3903 104
- * * * * * * * * * * * * * * * * * * * *
- * * * * * * * * 0
//
Output files for usage example 3
File: myhmms
HMMER2.0 [2.3.2]
NAME rrm
LENG 77
ALPH Amino
RF no
CS no
MAP yes
COM hmmbuild -n rrm --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F myhmms ../../data/hmmnew/rrm.sto
NSEQ 90
DATE Tue Jul 15 12:00:00 2008
CKSUM 8325
XT -8455 -4 -1000 -1000 -8455 -4 -8455 -4
NULT -4 -8455
NULE 595 -1558 85 338 -294 453 -1158 197 249 902 -1085 -142 -21 -313 45 531 201 384 -1998 -644
HMM A C D E F G H I K L M N P Q R S T V W Y
m->m m->i m->d i->m i->i d->m d->d b->m m->e
-16 * -6492
1 -1084 390 -8597 -8255 -5793 -8424 -8268 2395 -8202 2081 -1197 -8080 -8115 -8020 -8297 -7789 -5911 1827 -7525 -7140 1
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 -16 *
2 -2140 -3785 -6293 -2251 3226 -2495 -727 -638 -2421 -545 -675 -5146 -5554 -4879 -1183 -2536 -1928 267 76 3171 2
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
3 -2542 458 -8584 -8273 -6055 -8452 -8531 2304 -8255 -324 101 -8104 -8170 -8221 -8440 -7840 -5878 3145 -7857 -7333 3
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
4 -1505 -5144 -1922 -558 -1842 2472 -3303 -2213 1099 -5160 -4233 372 -4738 -530 1147 168 498 -4766 -5327 -1476 4
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
5 -3724 -5184 300 -3013 -1655 1803 -3353 -5245 -1569 -2686 -4276 3495 -1963 -1331 -1054 -1472 -3664 -4803 -5369 -2 5
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
6 -1569 -6106 -8967 -8363 555 -8531 -7279 654 -8092 2953 -94 -8220 -7908 -1643 -7682 -7771 -6460 -59 -6191 -6284 6
- -151 -504 230 45 -380 399 101 -621 211 -470 -713 278 399 48 91 360 113 -364 -299 -254
- -178 -3113 -12684 -1600 -578 -701 -1378 * *
7 -409 -5130 -215 -2987 -1709 -956 690 -5188 -395 -5144 -4224 729 3054 -2862 -3409 354 1293 -1381 -5321 -4644 13
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
8 -3674 -5118 -1004 639 420 -4652 176 -2050 404 -1039 -935 16 1755 168 147 -275 198 -1472 1889 1977 14
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
9 -408 -5134 2415 1299 -950 -66 -767 -1296 -2889 -1843 -4224 1084 -968 -1439 -1854 540 -314 -2304 -5320 -60 15
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
10 586 1804 -6294 -631 -1627 -1671 -4374 1029 -2223 -162 1172 -5147 -5554 -1870 -5058 -2327 1741 1687 -4242 687 16
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
11 -2134 -5144 845 -1187 -1652 -1667 -3303 -5216 -513 -801 -4233 1026 -1873 -543 -619 575 2956 -4766 -5327 -4644 17
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
[Part of this file has been deleted for brevity]
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11935 -12977 -894 -1115 -701 -1378 * *
69 -8583 -7493 -8946 -9308 -409 -8830 -5009 -7467 -8865 -3133 -6867 -2231 -8684 -7580 -8224 -8080 -8432 -7625 -4255 4885 98
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11935 -12977 -894 -1115 -701 -1378 * *
70 -1815 -5406 -2090 1238 -1053 -2697 642 -683 206 -981 -1214 1027 -5018 797 -56 458 1323 686 -5593 -187 99
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11935 -12977 -894 -1115 -701 -1378 * *
71 -1772 -576 -6691 -6059 2751 -2721 -4781 1861 -5659 107 276 -5548 -5946 -5283 -5462 -4990 -4280 2105 -4641 656 100
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11935 -12977 -894 -1115 -701 -1378 * *
72 -1774 157 -771 -126 -5740 -2696 -247 -1555 925 -5435 -361 -242 -5015 1187 2154 813 923 -1324 1141 185 101
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
E -1 -11935 -12977 -894 -1115 -701 -1378 * *
73 -1859 -5823 -9013 -8663 -1043 -8868 -8653 1700 -8616 373 -1223 -8529 -8476 -8301 -8657 -8252 -6326 3293 -7705 -7445 102
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
E -1 -11935 -12977 -894 -1115 -701 -1378 * *
74 339 -912 -2304 -1480 -161 -4982 242 328 485 -1007 -4340 -1086 -5073 1886 1206 -1252 451 493 -5463 1681 103
- -149 -500 233 43 -381 398 105 -626 210 -462 -721 275 394 45 96 359 117 -369 -295 -250
E -43 -5729 -6603 -86 -4119 -701 -1378 * *
75 3003 -4642 -7089 -2701 -4985 -369 -5478 -2243 -6266 -1473 -4220 -5882 -6177 -5887 -6084 548 1259 -953 -5472 -5145 105
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
E -1 -11920 -12962 -894 -1115 -357 -2190 * *
76 -890 -701 -6337 -456 -128 -5749 302 444 -435 294 -981 -5315 -5801 452 1013 -1812 -1348 2281 -4547 1032 106
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
E -1 -11935 -12977 -894 -1115 -701 -1378 * *
77 -2363 -782 199 -2083 -5740 -2662 -943 -1132 147 -3104 -4508 2805 -755 -1553 -765 1511 1258 -2796 -5603 -104 107
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
E -14 -11935 -6734 -894 -1115 -701 -1378 * *
78 349 -725 812 169 -2094 1817 -228 -2472 -192 -5424 -4497 321 -2426 905 63 -120 421 -272 -5591 -4908 108
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11921 -12963 -894 -1115 -373 -2135 * *
79 763 -866 480 290 -2215 944 426 219 379 -1687 -278 514 -2124 507 130 -1212 -544 72 -5596 -747 109
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -733 -11935 -1329 -894 -1115 -701 -1378 * *
80 -1113 -4530 -1626 -702 -1338 3243 -3172 -4385 -2855 -2220 -3677 -1218 -4564 -2783 -3335 -609 475 -4086 -4819 473 110
- -153 -500 231 50 -364 398 93 -624 206 -464 -733 283 383 46 94 356 116 -368 -250 -237
C -6367 -18 -12245 -485 -1807 -3912 -99 * *
81 -5380 -6068 -5956 -5194 -7389 3194 -752 -6799 -1418 -2482 -5840 -5063 -706 -4120 -1351 1839 -5344 -6386 -6417 -6291 116
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11248 -12290 -894 -1115 -277 -2516 * *
82 68 -5329 -386 943 -1053 -34 1308 -2479 103 -1629 757 -877 1926 1130 381 -1141 -878 -1038 -5513 -1828 117
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11840 -12882 -894 -1115 -101 -3893 * *
83 102 -5296 -3867 1691 21 -1212 -3624 -1510 -259 -602 -1204 -1862 1786 -1626 665 -353 -181 -4888 2391 466 118
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
C -1 -11935 -12977 -894 -1115 -701 -1378 * *
84 -1211 -5433 -3815 -1455 -5755 -850 -140 -2313 -1280 -5448 -4523 -3575 -2241 -1598 -412 3271 384 -5056 -5615 -4936 119
- * * * * * * * * * * * * * * * * * * * *
C * * * * * * * * 0
//
Output files for usage example 4
File: myhmms
HMMER2.0 [2.3.2]
NAME rrm
LENG 77
ALPH Amino
RF no
CS no
MAP yes
COM hmmbuild -n rrm --pbswitch 1000 --archpri 0.850000 --idlevel 0.620000 --swentry 0.500000 --swexit 0.500000 --wgsc -A -F myhmms ../../data/hmmnew/rrm.sto
NSEQ 90
DATE Tue Jul 15 12:00:00 2008
CKSUM 8325
XT -8455 -4 -1000 -1000 -8455 -4 -8455 -4
NULT -4 -8455
NULE 595 -1558 85 338 -294 453 -1158 197 249 902 -1085 -142 -21 -313 45 531 201 384 -1998 -644
HMM A C D E F G H I K L M N P Q R S T V W Y
m->m m->i m->d i->m i->i d->m d->d b->m m->e
-16 * -6492
1 -1084 390 -8597 -8255 -5793 -8424 -8268 2395 -8202 2081 -1197 -8080 -8115 -8020 -8297 -7789 -5911 1827 -7525 -7140 1
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 -16 *
2 -2140 -3785 -6293 -2251 3226 -2495 -727 -638 -2421 -545 -675 -5146 -5554 -4879 -1183 -2536 -1928 267 76 3171 2
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
3 -2542 458 -8584 -8273 -6055 -8452 -8531 2304 -8255 -324 101 -8104 -8170 -8221 -8440 -7840 -5878 3145 -7857 -7333 3
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
4 -1505 -5144 -1922 -558 -1842 2472 -3303 -2213 1099 -5160 -4233 372 -4738 -530 1147 168 498 -4766 -5327 -1476 4
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
5 -3724 -5184 300 -3013 -1655 1803 -3353 -5245 -1569 -2686 -4276 3495 -1963 -1331 -1054 -1472 -3664 -4803 -5369 -2 5
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
6 -1569 -6106 -8967 -8363 555 -8531 -7279 654 -8092 2953 -94 -8220 -7908 -1643 -7682 -7771 -6460 -59 -6191 -6284 6
- -151 -504 230 45 -380 399 101 -621 211 -470 -713 278 399 48 91 360 113 -364 -299 -254
- -178 -3113 -12684 -1600 -578 -701 -1378 * *
7 -409 -5130 -215 -2987 -1709 -956 690 -5188 -395 -5144 -4224 729 3054 -2862 -3409 354 1293 -1381 -5321 -4644 13
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
8 -3674 -5118 -1004 639 420 -4652 176 -2050 404 -1039 -935 16 1755 168 147 -275 198 -1472 1889 1977 14
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
9 -408 -5134 2415 1299 -950 -66 -767 -1296 -2889 -1843 -4224 1084 -968 -1439 -1854 540 -314 -2304 -5320 -60 15
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
10 586 1804 -6294 -631 -1627 -1671 -4374 1029 -2223 -162 1172 -5147 -5554 -1870 -5058 -2327 1741 1687 -4242 687 16
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
11 -2134 -5144 845 -1187 -1652 -1667 -3303 -5216 -513 -801 -4233 1026 -1873 -543 -619 575 2956 -4766 -5327 -4644 17
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11642 -12684 -894 -1115 -701 -1378 * *
[Part of this file has been deleted for brevity]
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11253 -12295 -894 -1115 -701 -1378 * *
279 -7207 -7306 -8076 -6588 -8459 -7223 -5448 -7982 -1500 -7531 -6953 -6369 -7277 -5081 4236 -7139 -6862 -7777 -7053 -7277 454
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11253 -12295 -894 -1115 -701 -1378 * *
280 -694 -163 -5922 -5286 -1204 -2048 -610 1082 -1800 1434 -2618 -4776 2951 -4509 -4688 -1216 -1648 -2829 202 21 455
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -3168 -11253 -171 -894 -1115 -701 -1378 * *
281 -1412 -2132 -2007 -2293 -4366 3113 -2847 -4225 -3107 -4377 -3503 1660 -2881 -2661 -3396 961 -1821 -3134 -4516 -4119 456
- -150 -489 232 42 -382 400 104 -627 211 -465 -722 274 393 51 95 359 116 -370 -296 -245
- -2121 -637 -2975 -831 -1191 -6099 -21 * *
282 -968 -1818 -1787 -1351 -3112 953 -1494 -2818 -1122 -2911 -2044 -1365 -2340 -1133 1510 1816 2121 -2205 -3137 -2649 459
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -228 -7899 -2816 -894 -1115 -4964 -47 * *
283 840 -1663 -994 969 1159 503 -604 -1413 -325 -1594 -814 -688 -1996 -267 1103 -851 -755 -1179 2900 -1437 460
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -9 -7938 -8980 -894 -1115 -89 -4060 * *
284 -3257 -4642 -697 -2590 -1218 -252 -2907 -4655 -1306 -2353 -529 482 -1607 -2459 -1398 2112 2745 -4246 -4848 -4187 461
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11183 -12226 -894 -1115 -186 -3045 * *
285 2163 763 -1619 -5296 2250 -2060 -4007 1241 -4891 -489 484 -4781 -226 -4515 -4692 -678 -1688 -813 264 -3530 462
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11253 -12295 -894 -1115 -701 -1378 * *
286 -268 -329 -158 917 -541 -1990 350 -4851 1273 -1075 388 -1130 233 840 993 -602 801 -595 -4964 -857 463
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11253 -12295 -894 -1115 -701 -1378 * *
287 109 -243 672 2304 -5103 -4283 488 -4854 -1317 -2269 -656 -492 -1519 2679 -655 -618 -3248 -4404 -4965 -1114 464
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11253 -12295 -894 -1115 -701 -1378 * *
288 1312 1294 -6215 -5593 -206 -1244 -4339 2188 -5201 1409 395 -5091 -5478 -4828 -5009 -4538 -3794 1162 -4188 -3846 465
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -23 -11253 -6022 -894 -1115 -701 -1378 * *
289 -3562 799 -5767 -2054 -1235 -2075 318 138 237 2164 1713 -1454 -5145 -1272 -730 -4172 -1640 1071 -3865 -34 466
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11231 -12273 -894 -1115 -1470 -646 * *
290 73 1351 -674 1236 -1549 -2008 1350 -4834 1049 -2498 -3851 1801 -4356 1813 -115 -223 -1582 -1052 -4945 -4262 467
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11231 -12273 -894 -1115 -369 -2147 * *
291 -1739 -320 777 -2654 -1419 -2051 4360 -4707 -1358 -2412 -689 -1300 -4399 -224 537 531 -289 -2010 -4905 -1057 468
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -11253 -12295 -894 -1115 -701 -1378 * *
292 -3345 -4494 -233 -332 -563 -1986 -3051 333 99 1063 -3616 -3072 2953 -1026 -1490 -943 -1528 -1070 -4753 -4151 469
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -10815 -11857 -894 -1115 -701 -1378 * *
293 -6409 -5751 -7614 -7636 2593 -7311 -4003 -5084 -7219 -150 -151 -6210 -7172 -849 -6723 -6510 -6299 -1387 4881 2807 470
- -149 -500 233 43 -381 399 106 -626 210 -466 -720 275 394 45 96 359 117 -369 -294 -249
- -1 -10749 -11791 -894 -1115 -701 -1378 * *
294 -4057 -3817 -6415 -5791 3203 -1638 -4541 1679 -5412 765 1434 -5333 -5617 -4930 -5182 -4791 -3987 1226 750 -3959 471
- * * * * * * * * * * * * * * * * * * * *
- * * * * * * * * 0
//
Data files
None.
Notes
1. Command-line arguments
The following original HMMER options are not supported:
-h : Use -help to get help information instead.
-f : Use -strategy option instead.
-g : Use -strategy option instead.
-s : Use -strategy option instead.
-A : Set append: "N" or append: "Y" for "hmmfile" in the ACD file instead.
-F : Always set (an existing hmmfile will be overwritten).
-amino : Sequence alignment type is specified via the ACD file.
-nucleic : Sequence alignment type is specified via the ACD file.
-informat : All common alignment formats are supported automatically.
-wblosum : Use -weighting option to specify the sequence weighting algorithm.
-wgsc : Use -weighting option to specify the sequence weighting algorithm.
-wme : Use -weighting option to specify the sequence weighting algorithm.
-wnone : Use -weighting option to specify the sequence weighting algorithm.
-wpb : Use -weighting option to specify the sequence weighting algorithm.
-wvoronoi : Use -weighting option to specify the sequence weighting algorithm.
-verbose : Use -verbosity instead.
-weighting : Sequence weighting algorithm.
-n : Use -nhmm instead (-n causes problems for GUI developers)
2. Installing EMBASSY HMMER
The EMBASSY HMMER package contains "wrapper" applications providing an EMBOSS-style interface to the applications in the original HMMER package version 2.3.2 developed by Sean Eddy. Please read the file INSTALL in the EMBASSY HMMER package distribution for installation instructions.
3. Installing original HMMER
To use EMBASSY HMMER, you will first need to download and install the original HMMER package. Please read the file 00README in the the original HMMER package distribution for installation instructions:
WWW home: http://hmmer.wustl.edu/
Distribution: ftp://ftp.genetics.wustl.edu/pub/eddy/hmmer/
4. Setting up HMMER
For the EMBASSY HMMER package to work, the directory containing the original HMMER executables *must* be in your path. For example if you executables were installed to "/usr/local/hmmer/bin", then type:
set path=(/usr/local/hmmer/bin/ $path)
rehash
5. Getting help
Please read the Userguide.pdf distributed with the original HMMER and included in the EMBASSY HMMER distribution under the DOCS directory. The first 3 chapters (Introduction, Installation and Tutorial) are particularly useful.
References
None.
Warnings
Types of input data
hmmer v3.2.1 and therefore EMBASSY HMMER is only recommended for use with protein sequences. If you provide a non-protein sequence you will be reprompted for a protein sequence. To accept nucleic acid sequences you must replace instances of < type: "protein" > in the application ACD files with Environment variables
The original hmmer uses BLAST environment variables (below), if defined, to locate files. The EMBASSY HMMER does not.
BLASTDB location of sequence databases to be searched
BLASMAT location of substitution matrices
HMMERDB location of HMMs
Alignment input
The user must provide the full filename of an alignment for the "alignfile" ACD option, not an indirect reference to a set of sequences, e.g. a USA is NOT acceptable. This is because hmmbuild (which ehmmbuild wraps) requires an alignment and does not support USAs.
Diagnostic Error Messages
None.
Exit status
It always exits with status 0.
Known bugs
None.
See also
Program name
Description
ehmmalign
Align sequences to an HMM profile
ehmmcalibrate
Calibrate HMM search statistics
ehmmconvert
Convert between profile HMM file formats
ehmmemit
Generate sequences from a profile HMM
ehmmfetch
Retrieve an HMM from an HMM database
ehmmindex
Create a binary SSI index for an HMM database
ehmmpfam
Search one or more sequences against an HMM database
ehmmsearch
Search a sequence database with a profile HMM
oalistat
Statistics for multiple alignment files
ohmmalign
Align sequences with an HMM
ohmmbuild
Build HMM
ohmmcalibrate
Calibrate a hidden Markov model
ohmmconvert
Convert between HMM formats
ohmmemit
Extract HMM sequences
ohmmfetch
Extract HMM from a database
ohmmindex
Index an HMM database
ohmmpfam
Align single sequence with an HMM
ohmmsearch
Search sequence database with an HMM
Author(s)
This program is an EMBOSS conversion of a program written by Sean Eddy
as part of his HMMER package.
European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK
History
Target users
This program is intended to be used by everyone and everything, from naive users to embedded scripts.