OnWorks favicon

hhsearch - Online in the Cloud

Run hhsearch in OnWorks free hosting provider over Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator

This is the command hhsearch that can be run in the OnWorks free hosting provider using one of our multiple free online workstations such as Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator



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


hhsearch -i query -d database [options]


HHsearch version 2.0.16 (January 2013) Search a database of HMMs with a query alignment or
query HMM (C) Johannes Soeding, Michael Remmert, Andreas Biegert, Andreas Hauser Soding,
J. Protein homology detection by HMM-HMM comparison. Bioinformatics 21:951-960 (2005).

-i <file>
input/query multiple sequence alignment (a2m, a3m, FASTA) or HMM

-d <file>
HMM database of concatenated HMMs in hhm, HMMER, or a3m format, OR, if file has
extension pal, list of HMM file names, one per line. Multiple dbs, HMMs, or pal
files with -d '<db1> <db2>...'

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

Output options:
-o <file>
write results in standard format to file (default=<infile.hhr>)

-Ofas <file>
write pairwise alignments of significant matches in FASTA format Analogous for
output in a3m, a2m, and psi format (e.g. -Oa3m)

-oa3m <file>
write MSA of significant matches in a3m format Analogous for output in a2m, psi,
and hhm format (e.g. -ohhm)

-e [0,1]
E-value cutoff for inclusion in multiple alignment (def=0.001)

-seq <int>
max. number of query/template sequences displayed (def=1) Beware of overflows! All
these sequences are stored in memory.

-cons show consensus sequence as master sequence of query MSA

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

-p <float>
minimum probability in summary and alignment list (def=20)

-E <float>
maximum E-value in summary and alignment list (def=1E+06)

-Z <int>
maximum number of lines in summary hit list (def=500)

-z <int>
minimum number of lines in summary hit list (def=10)

-B <int>
maximum number of alignments in alignment list (def=500)

-b <int>
minimum number of alignments in alignment list (def=10)

-aliw [40,..[
number of columns per line in alignment list (def=80)

max length of database string to be printed in hhr file

Filter query multiple sequence alignment

-id [0,100] maximum pairwise sequence identity (%) (def=90)

-diff [0,inf[
filter MSA by selecting most diverse set of sequences, keeping at least this many
seqs in each MSA block of length 50 (def=100)

-cov [0,100] minimum coverage with query (%) (def=0)

-qid [0,100] minimum sequence identity with query (%) (def=0)

-qsc [0,100] minimum score per column with query (def=-20.0)

-neff [1,inf]
target diversity of alignment (default=off)

Input alignment format:
-M a2m use A2M/A3M (default): upper case = Match; lower case = Insert; '-' = Delete; '.' =
gaps aligned to inserts (may be omitted)

-M first
use FASTA: columns with residue in 1st sequence are match states

-M [0,100]
use FASTA: columns with fewer than X% gaps are match states

-tags do NOT neutralize His-, C-myc-, FLAG-tags, and trypsin recognition sequence to
background distribution

HMM-HMM alignment options:
do NOT realign displayed hits with MAC algorithm (def=realign)

-mact [0,1[
posterior probability threshold for MAC re-alignment (def=0.350) Parameter controls
alignment greediness: 0:global >0.1:local

use global/local alignment mode for searching/ranking (def=local)

-alt <int>
show up to this many significant alternative alignments(def=2)

-vit use Viterbi algorithm for searching/ranking (default)

-mac use Maximum Accuracy (MAC) algorithm for searching/ranking

use Forward probability for searching

-excl <range>
exclude query positions from the alignment, e.g. '1-33,97-168'

-shift [-1,1]
score offset (def=-0.03)

-corr [0,1]
weight of term for pair correlations (def=0.10)

-sc <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

-ssm {0,..,4}
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

-ssw [0,1] weight of ss score compared to column score (def=0.11)

-ssa [0,1] SS substitution matrix = (1-ssa)*I + ssa*full-SS-substition-matrix

Gap cost options:
-gapb [0,inf[
Transition pseudocount admixture (def=1.00)

-gapd [0,inf[
Transition pseudocount admixture for open gap (default=0.15)

-gape [0,1.5]
Transition pseudocount admixture for extend gap (def=1.00)

-gapf ]0,inf]
factor to increase/reduce the gap open penalty for deletes (def=0.60)

-gapg ]0,inf]
factor to increase/reduce the gap open penalty for inserts (def=0.60)

-gaph ]0,inf]
factor to increase/reduce the gap extend penalty for deletes(def=0.60)

-gapi ]0,inf]
factor to increase/reduce the gap extend penalty for inserts(def=0.60)

-egq [0,inf[ penalty (bits) for end gaps aligned to query residues (def=0.00)

-egt [0,inf[ penalty (bits) for end gaps aligned to template residues (def=0.00)

Pseudocount (pc) options:
-pcm {0,..,3}
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)
3: constant diversity pseudocounts

-pca [0,1] overall pseudocount admixture (def=1.0)

-pcb [1,inf[ Neff threshold value for -pcm 2 (def=1.5)

-pcc [0,3] extinction exponent c for -pcm 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

<file> column state file for fast database prefiltering (default=./data/cs219.lib)

-csw [0,inf] weight of central position in cs pseudocount mode (def=1.6)

-csb [0,1] weight decay parameter for positions in cs pc mode (def=0.9)

Other options:
-cpu <int>
number of CPUs to use (for shared memory SMPs) (default=1)

-v <int>
verbose mode: 0:no screen output 1:only warings 2: verbose

-maxres <int>
max number of HMM columns (def=15002)

-maxmem [1,inf[ max available memory in GB (def=3.0)

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

-calm {0,..,3} empirical score calibration of 0:query 1:template 2:both

default 3: neural network-based estimation of EVD params

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

Use hhsearch online using onworks.net services

Free Servers & Workstations

Download Windows & Linux apps

  • 1
    Fork of TeamWinRecoveryProject(TWRP)
    with many additional functions, redesign
    and more Features:Supports Treble and
    non-Treble ROMsUp-to-date Oreo kernel,
    Download OrangeFox
  • 2
    itop - ITSM  CMDB OpenSource
    itop - ITSM CMDB OpenSource
    IT Operations Portal: a complete open
    source, ITIL, web based service
    management tool including a fully
    customizable CMDB, a helpdesk system and
    a document man...
    Download itop - ITSM CMDB OpenSource
  • 3
    Clementine is a multi-platform music
    player and library organizer inspired by
    Amarok 1.4. It has a fast and
    easy-to-use interface, and allows you to
    search and ...
    Download Clementine
  • 4
    ATTENTION: Cumulative update 2.4.3 has
    been released!! The update works for any
    previous 2.x.x version. If upgrading
    from version v1.x.x, please download and
    Download XISMuS
  • 5
    Modular headtracking program that
    supports multiple face-trackers, filters
    and game-protocols. Among the trackers
    are the SM FaceAPI, AIC Inertial Head
    Tracker ...
    Download facetracknoir
  • 6
    PHP QR Code
    PHP QR Code
    PHP QR Code is open source (LGPL)
    library for generating QR Code,
    2-dimensional barcode. Based on
    libqrencode C library, provides API for
    creating QR Code barc...
    Download PHP QR Code
  • 7
    Cuckoo Sandbox
    Cuckoo Sandbox
    Cuckoo Sandbox uses components to
    monitor the behavior of malware in a
    Sandbox environment; isolated from the
    rest of the system. It offers automated
    analysis o...
    Download Cuckoo Sandbox
  • More »

Linux commands