nhmmer
Quick Start
- Command:
nhmmer - Local executable:
/home/vimalinx/miniforge3/envs/bio/bin/nhmmer - Version: HMMER 3.4
- Full reference: See
references/help.md
When To Use This Tool
- Search DNA or RNA queries against nucleotide sequence databases with HMMER sensitivity.
- Look for more remote nucleotide homologs than simple pairwise nucleotide aligners tend to recover.
- Restrict searches to the Watson or Crick strand when the biology justifies it.
- Prefer
nhmmscanwhen the target is a nucleotide HMM database rather than a plain sequence database.
Common Patterns
# 1) Search a nucleotide query against a genomic FASTA database
nhmmer \
--tblout hits.tbl \
--cpu 8 \
query.fa \
genome.fa
# 2) Search only one strand
nhmmer \
--watson \
--tblout watson_hits.tbl \
query.fa \
genome.fa
# 3) Use curated model thresholds and save Dfam-style output
nhmmer \
--cut_ga \
--dfamtblout dfam.tbl \
query.fa \
genome.fa
Recommended Workflow
- Start from a nucleotide query and a nucleotide target database in FASTA or another supported sequence format.
- Save parseable output with
--tblout, and use--dfamtbloutif you are working in a Dfam-like annotation workflow. - Apply strand restriction or curated thresholds only when the underlying assay or model supports those assumptions.
- Review significant hits before deciding whether to promote the query into a reusable profile model.
Guardrails
- In this workspace the binary currently fails to start because
libopenblas.so.0is missing. nhmmerexpects DNA or RNA queries and DNA/RNA targets; it is not the protein HMMER engine.- If the query comes from stdin, the local binary expects
--qformatto be set explicitly. - Use
--watsonor--crickdeliberately; do not silently halve the search space unless strand specificity is real.