Scroll to navigation

HHSEARCH(1) User Commands HHSEARCH(1)

NAME

hhsearch - search a database of HMMs with a query alignment or query HMM

SYNOPSIS

hhsearch -i query -d database [options]

DESCRIPTION

HHsearch 3.3.0 Search a database of HMMs with a query alignment or query HMM (c) The HH-suite development team Steinegger M, Meier M, Mirdita M, V??hringer H, Haunsberger S J, and S??ding J (2019) HH-suite3 for fast remote homology detection and deep protein annotation. BMC Bioinformatics, doi:10.1186/s12859-019-3019-7

input/query multiple sequence alignment (a2m, a3m, FASTA) or HMM

<file> may be 'stdin' or 'stdout' throughout. Options:

database name (e.g. uniprot20_29Feb2012) Multiple databases may be specified with '-d <db1> -d <db2> ...'
[0,1] E-value cutoff for inclusion in result alignment (def=0.001)

Input alignment format:

use A2M/A3M (default): upper case = Match; lower case = Insert;
'-' = Delete; '.' = gaps aligned to inserts (may be omitted)
use FASTA: columns with residue in 1st sequence are match states
use FASTA: columns with fewer than X% gaps are match states
do NOT / do neutralize His-, C-myc-, FLAG-tags, and trypsin recognition sequence to background distribution (def=-notags)

Output options:

write results in standard format to file (default=<infile.hhr>)
write result MSA with significant matches in a3m format

-blasttab <name> write result in tabular BLAST format (compatible to -m 8 or -outfmt 6 output)

1
2 3 4 5 6 7 8 9 10 11 12
query target #match/tLen alnLen #mismatch #gapOpen qstart qend tstart tend eval score
write result MSA of significant matches in PSI-BLAST format
write HHM file for result MSA of significant matches
generate consensus sequence as master sequence of query MSA (default=don't)
don't show consensus sequence in alignments (default=show)
don't show predicted 2ndary structure in alignments (default=show)
don't show DSSP 2ndary structure in alignments (default=show)
show confidences for predicted 2ndary structure in alignments
write pairwise alignments in FASTA xor A2M (-Oa2m) xor A3M (-Oa3m) format
max. number of query/template sequences displayed (default=1)
number of columns per line in alignment list (default=80)
minimum probability in summary and alignment list (default=20)
maximum E-value in summary and alignment list (default=1E+06)
maximum number of lines in summary hit list (default=500)
minimum number of lines in summary hit list (default=10)
maximum number of alignments in alignment list (default=500)
minimum number of alignments in alignment list (default=10)

Filter options applied to query MSA, database MSAs, and result MSA

show all sequences in result MSA; do not filter result MSA
[0,100] maximum pairwise sequence identity (def=90)
filter MSAs by selecting most diverse set of sequences, keeping at least this many seqs in each MSA block of length 50 Zero and non-numerical values turn off the filtering. (def=100)
[0,100] minimum coverage with master sequence (%) (def=0)
[0,100] minimum sequence identity with master sequence (%) (def=0)
[0,100] minimum score per column with master sequence (default=-20.0)
target diversity of multiple sequence alignment (default=off)
do not filter out sequences marked by ">@"in their name line

HMM-HMM alignment options:

do NOT realign displayed hits with MAC algorithm (def=realign)
banded alignment: forbid <ovlp> largest diagonals |i-j| of DP matrix (def=0)
posterior prob threshold for MAC realignment controlling greediness at alignment ends: 0:global >0.1:local (default=0.35)
use global/local alignment mode for searching/ranking (def=local)
realign displayed hits with max. accuracy (MAC) algorithm
exclude query positions from the alignment, e.g. '1-33,97-168'
realign max. <int> hits (default=500)
show up to this many alternative alignments with raw score > smin(def=4)
minimum raw score for alternative alignments (def=20.0)
profile-profile score offset (def=-0.03)
weight of term for pair correlations (def=0.10)
<int> amino acid score (tja: template HMM at column j) (def=1)
0
= log2 Sum(tja*qia/pa) (pa: aa background frequencies)
1
= log2 Sum(tja*qia/pqa) (pqa = 1/2*(pa+ta) )
2
= log2 Sum(tja*qia/ta) (ta: av. aa freqs in template)
3
= log2 Sum(tja*qia/qa) (qa: av. aa freqs in query)
5
local amino acid composition correction
0: no ss scoring
1,2: ss scoring after or during alignment
[default=2]
3,4: ss scoring after or during alignment, predicted vs. predicted
weight of ss score (def=0.11)
SS substitution matrix = (1-ssa)*I + ssa*full-SS-substition-matrix [def=1.00)
use global sequence weighting for realignment!

Gap cost options:

Transition pseudocount admixture (def=1.00)
Transition pseudocount admixture for open gap (default=0.15)
Transition pseudocount admixture for extend gap (def=1.00)
factor to increase/reduce gap open penalty for deletes (def=0.60)
factor to increase/reduce gap open penalty for inserts (def=0.60)
factor to increase/reduce gap extend penalty for deletes(def=0.60)
factor to increase/reduce gap extend penalty for inserts(def=0.60)
[0,inf[ penalty (bits) for end gaps aligned to query residues (def=0.00)
[0,inf[ penalty (bits) for end gaps aligned to template residues (def=0.00)

Pseudocount (pc) options:

Context specific hhm pseudocounts:
position dependence of pc admixture 'tau' (pc mode, default=2)
0: no pseudo counts:
tau = 0
1: constant
tau = a
2: diversity-dependent: tau = a/(1+((Neff[i]-1)/b)^c) 3: CSBlast admixture: tau = a(1+b)/(Neff[i]+b) (Neff[i]: number of effective seqs in local MSA around column i)
[0,1] overall pseudocount admixture (def=0.9)
[1,inf[ Neff threshold value for mode 2 (def=4.0)
[0,3] extinction exponent c for mode 2 (def=1.0)
Context independent hhm pseudocounts (used for templates; used for query if contxt file is not available):
position dependence of pc admixture 'tau' (pc mode, default=2)
0: no pseudo counts:
tau = 0
1: constant
tau = a
2: diversity-dependent: tau = a/(1+((Neff[i]-1)/b)^c) (Neff[i]: number of effective seqs in local MSA around column i)
[0,1] overall pseudocount admixture (def=1.0)
[1,inf[ Neff threshold value for mode 2 (def=1.5)
[0,3] extinction exponent c for mode 2 (def=1.0)
Context-specific pseudo-counts:
use substitution-matrix instead of context-specific pseudocounts

-contxt <file> context file for computing context-specific pseudocounts (default=)

[0,inf] weight of central position in cs pseudocount mode (def=1.6)
[0,1] weight decay parameter for positions in cs pc mode (def=0.9)

Other options:

verbose mode: 0:no screen output 1:only warnings 2: verbose (def=2)
number of CPUs to use (for shared memory SMPs) (default=2)

-scores <file> write scores for all pairwise comparisons to file

<file> write all alignments in tabular layout to file
max number of input rows (def=65535)
max number of HMM columns (def=20001)

-maxmem [1,inf[ limit memory for realignment (in GB) (def=3.0)

Example: hhsearch -i a.1.1.1.a3m -d scop70_1.71

Download databases from <http://wwwuser.gwdg.de/~compbiol/data/hhsuite/databases/hhsuite_dbs/>.

May 2022 hhsearch 3.3.0+ds