hmmer2(1)

🌸 hmmer2 β€” the HMMER manual, but make it a group chat 🌸

okay bestie, so this is the manual for HMMER (that's the profile hidden Markov model software), living in section 1 of the manual. buckle up, it's giving biology but make it statistical πŸ’…


✨ NAME ✨

HMMER β€” profile hidden Markov model software πŸ§¬πŸ’–


🌈 SYNOPSIS 🌈

so there's a whole squad of little programs and they all do their own thing. here's the roster, no cap:

hmm2align
  • hmm2align β€” Align multiple sequences to a profile HMM. 🫢
hmm2build
  • hmm2build β€” Build a profile HMM from a given multiple sequence alignment. πŸ—οΈβœ¨
hmm2calibrate
  • hmm2calibrate β€” Determine appropriate statistical significance parameters for a profile HMM prior to doing database searches. πŸ“Š
hmm2convert
  • hmm2convert β€” Convert HMMER profile HMMs to other formats, such as GCG profiles. πŸ”„
hmm2emit
  • hmm2emit β€” Generate sequences probabilistically from a profile HMM. 🎲
hmm2fetch
  • hmm2fetch β€” Retrieve an HMM from an HMM database. πŸ“₯
hmm2index
  • hmm2index β€” Create a binary SSI index for an HMM database. πŸ—‚οΈ
hmm2pfam
  • hmm2pfam β€” Search a profile HMM database with a sequence (i.e., annotate various kinds of domains in the query sequence). πŸ”
hmm2search
  • hmm2search β€” Search a sequence database with a profile HMM (i.e., find additional homologues of a modeled family). πŸ‘―

πŸ¦‹ DESCRIPTION πŸ¦‹

These programs use profile hidden Markov models (profile HMMs, iykyk) to model the primary structure consensus of a family of protein or nucleic acid sequences. literally that's the whole vibe. πŸ§¬πŸ’œ


πŸ’… OPTIONS πŸ’…

okay so here's the tea. every single HMMER program, if you just run it with no arguments at all, will give you a brief lil summary of its command-line syntax and options. super chill. πŸŽ€

But the moment you hit it with the single argument -h (that's help, bestie), the program gets verbose β€” it reports way more command-line usage info, including the rarely used, experimental, and expert options. πŸ™ˆ -h will also report version numbers, which are useful if you need to report a bug or a problem to me. πŸ‘‘

Each HMMER program has its own man page that briefly summarizes command line usage. There's also a user's guide that came with the software distribution β€” it 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. 🌐

In general, no command line options should be needed by beginning users. The defaults are set up for optimum performance in most situations. periodt. πŸ’…

here's how to read the option energy:

  • Options that are single lowercase letters (e.g. -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) are usually less common options, but also may be important in some applications. πŸ”₯
  • Options that are full words (e.g. --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. so like… use those at your own risk bestie 😭

πŸŽ€ SEQUENCE FILE FORMATS πŸŽ€

In general, 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. iconic behaviour honestly 🌈

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. πŸ“–

Sequence alignments

Multiple sequence alignments may be in CLUSTALW, SELEX, or GCG MSF format. These formats are documented in the User's Guide. πŸ“–


🧚 ENVIRONMENT VARIABLES 🧚

For ease of using large stable sequence and HMM databases, HMMER looks for sequence files and HMM files in the current working directory as well as in system directories specified by environment variables. 🫢

  • BLASTDB β€” Specifies the directory location of sequence databases. Example: /seqlibs/blast-db/. In installations that use BLAST software, this environment variable is likely to already be set. πŸ’Ύ
  • HMMERDB β€” Specifies the directory location of HMM databases. Example: /seqlibs/pfam/. πŸ—„οΈ

and that's the whole page bestie, you ate it πŸ’– this was HMMER 2.3.2, from Oct 2003, documented in the HMMER Manual, section hmmer2(1) β€” go forth and profile-hide those Markov models, mwah πŸŽ€βœ¨πŸ‘‘