|
HMMER
User's Guide
|
|
Dept. of Genetics |
WashU |
Medical School |
Sequencing Center |
CGM |
IBC|
|
Eddy lab |
Internal (lab only) |
HMMER |
PFAM |
tRNAscan-SE |
Software |
Publications
|
Next: hmmbuild - build a
Up: Manual pages
Previous: HMMER - profile hidden
Subsections
hmmalign
[options] hmmfile seqfile
hmmalign reads an HMM file from
hmmfile and a set of sequences from seqfile, aligns the sequences to
the profile HMM, and outputs a multiple sequence alignment.
seqfile
may be in any unaligned or aligned file format accepted by HMMER. If it
is in a multiple alignment format (e.g. MSF, SELEX, ClustalW), the existing
alignment is ignored.
- [-h ] Print brief help; includes version
number and summary of all options, including expert options.
- [-m ] Include
in the alignment only those symbols aligned to match states. Do not show
symbols assigned to insert states.
- [-o <f> ] Save alignment to file <f> instead
of to standard output.
- [-q ] quiet; suppress all output except the alignment
itself. Useful for piping or redirecting the output.
- [-mapali
<f> ] Reads an alignment from file <f> and aligns it as a single object to
the HMM; e.g. the alignment in <f> is held fixed. This allows you to align
sequences to a model with hmmalign and view them in the context of an
existing trusted multiple alignment. The alignment to the alignment is
defined by a "map" kept in the HMM, and so is fast and guaranteed to be
consistent with the way the HMM was constructed from the alignment. The
alignment in the file <f> must be exactly the alignment that the HMM was
built from. Compare the -withali option.
- [-withali <f> ] Reads an alignment
from file <f> and aligns it as a single object to the HMM; e.g. the alignment
in <f> is held fixed. This allows you to align sequences to a model with
hmmalign and view them in the context of an existing trusted multiple
alignment. The alignment to the alignment is done with a heuristic (nonoptimal)
dynamic programming procedure, which may be somewhat slow and is not guaranteed
to be completely consistent with the way the HMM was constructed (though
it should be quite close). However, any alignment can be used, not just
the alignment that the HMM was built from. Compare the -mapali option.
Next: hmmbuild - build a
Up: Manual pages
Previous: HMMER - profile hidden
Direct comments and questions to <eddy@genetics.wustl.edu>