hmmsearch - search a sequence database with a profile HMM
HMMER 2.3.2 (Oct 2003)
Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file: /home/projects/MicrobialGenomicsGroup/Packages/hmmdata-1.0.0.package/Resources/noarch/TCS/HWE_HK.hmm [HWE_HK]
Sequence database: CP000828.Gprot.fsa
per-sequence score cutoff: >= 60.0 [TC1]
per-domain score cutoff: >= 60.0 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HWE_HK
Accession: PF07536.4
Description: HWE histidine kinase
[HMM has been calibrated; E-values are empirical estimates]
Scores for complete sequences (score includes all domains):
Sequence Description Score E-value N
-------- ----------- ----- ------- ---
CP000828_5747 CDS_5954224-5956767 92.6 3.5e-31 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000828_5747 1/1 521 603 .. 1 85 [] 92.6 3.5e-31
Alignments of top-scoring domains:
CP000828_5747: domain 1 of 1, from 521 to 603: score 92.6, E = 3.5e-31
*->ELnHRVKNtLAvVQSiarQTlRnaasldefveaFsgRLqALsraHdL
ELnHRVKN+LA++ Si rQT+++++s++e+ ++ R+ AL+ aHdL
CP000828_5 521 ELNHRVKNILALIRSISRQTRHSTTSVEEYTLLLDRRIEALAHAHDL 567
VpLtrseWagAdLreLleaeLaPyggedaeRitlsGPd<-*
++ + ++ L Ll+ L+Py ++ +++ +sGP+
CP000828_5 568 --VSGHGLEWPNLQNLLTTGLRPYLAKPKHQVSFSGPE 603
Histogram of all scores:
score obs exp (one = represents 18 sequences)
----- --- ---
-15 9 0|=
-14 115 0|=======
-13 275 0|================
-12 878 0|=================================================
-11 808 51|==*==========================================
-10 1019 646|===================================*=====================
-9 948 1600|===================================================== *
-8 912 1661|=================================================== *
-7 527 1116|============================== *
-6 261 609|=============== *
-5 235 301|============== *
-4 109 143|=======*
-3 67 66|===*
-2 42 30|=*=
-1 27 14|*=
0 13 6|*
1 5 2|*
2 1 1|*
3 0 0|
4 0 0|
> 5 3 -|=
% Statistical details of theoretical EVD fit:
mu = -7.9989
lambda = 0.7843
chi-sq statistic = 12584.0605
P(chi-square) = 0
Total sequences searched: 6254
Whole sequence top hits:
tophits_s report:
Total hits: 1
Satisfying E cutoff: 1
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 1
Satisfying E cutoff: 1
Total memory: 21K
