1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
|
.TH "hmmer2" 1 "Oct 2003" "HMMER 2.3.2" "HMMER Manual"
.SH NAME
.TP
HMMER - profile hidden Markov model software
.SH SYNOPSIS
.TP
.B hmm2align
Align multiple sequences to a profile HMM.
.TP
.B hmm2build
Build a profile HMM from a given multiple sequence alignment.
.TP
.B hmm2calibrate
Determine appropriate statistical significance parameters
for a profile HMM prior to doing database searches.
.TP
.B hmm2convert
Convert HMMER profile HMMs to other formats, such as GCG profiles.
.TP
.B hmm2emit
Generate sequences probabilistically from a profile HMM.
.TP
.B hmm2fetch
Retrieve an HMM from an HMM database
.TP
.B hmm2index
Create a binary SSI index for an HMM database
.TP
.B hmm2pfam
Search a profile HMM database with a sequence (i.e., annotate various
kinds of domains in the query sequence).
.TP
.B hmm2search
Search a sequence database with a profile HMM (i.e., find additional
homologues of a modeled family).
.SH DESCRIPTION
These programs use profile hidden Markov models (profile HMMs) to
model the primary structure consensus of a family of protein or
nucleic acid sequences.
.SH OPTIONS
.PP
All
.B HMMER
programs give a brief summary of their command-line syntax and options
if invoked without any arguments.
When invoked with the single argument,
.B -h
(i.e., help), a program will report more verbose command-line usage
information, including rarely used, experimental, and expert options.
.B -h
will report version numbers which are useful if
you need to report a bug or problem to me.
.PP
Each
.B HMMER
program has its own man page briefly summarizing command line usage.
There is also a user's guide that came
with the software distribution, which includes a tutorial introduction
and more detailed descriptions of the programs.
See http://hmmer.janelia.org/ for on-line documentation and
the current HMMER release.
.PP
In general, no command line options should be needed by beginning users.
The defaults are set up for optimum performance in most situations.
Options that are single lowercase letters (e.g.
.B -a
) are "common" options that are expected to be frequently used
and will be important in many applications.
Options that are single uppercase letters (e.g.
.B -B
) are usually less common options, but also may be important
in some applications.
Options that are full words (e.g.
.B --verbose
) are either rarely used, experimental, or expert options.
Some experimental options are only there for my own ongoing experiments
with HMMER, and may not be supported or documented adequately.
.SH SEQUENCE FILE FORMATS
In general,
.B HMMER
attempts to read most common biological sequence file formats.
It autodetects the format of the file. It also autodetects
whether the sequences are protein or nucleic acid.
Standard IUPAC degeneracy codes are allowed in addition
to the usual 4-letter or 20-letter codes.
.TP
.B Unaligned sequences
Unaligned sequence files may be in FASTA, Swissprot, EMBL, GenBank,
PIR, Intelligenetics, Strider, or GCG format.
These formats
are documented in the User's Guide.
.TP
.B Sequence alignments
Multiple sequence alignments may be in CLUSTALW, SELEX, or GCG MSF
format. These formats
are documented in the User's Guide.
.SH ENVIRONMENT VARIABLES
For ease of using large stable sequence and HMM databases,
.B HMMER
looks for sequence files and HMM files in the current
working directory as well as in system directories specified
by environment variables.
.TP
.B BLASTDB
Specifies the directory location of sequence databases. Example:
.B /seqlibs/blast-db/.
In installations that use BLAST software, this environment variable
is likely to already be set.
.TP
.B HMMERDB
Specifies the directory location of HMM databases. Example:
.B /seqlibs/pfam/.
|
