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: CP001578.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
-------- ----------- ----- ------- ---
CP001578_1607 CDS_1613407-1612472 146.9 3.6e-50 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001578_1607 1/1 118 200 .. 1 85 [] 146.9 3.6e-50
Alignments of top-scoring domains:
CP001578_1607: domain 1 of 1, from 118 to 200: score 146.9, E = 3.6e-50
*->ELnHRVKNtLAvVQSiarQTlRnaasldefveaFsgRLqALsraHdL
E++HR+KN+LA+VQSia+QT+R+++s+d+f+ +F+gR+q+Ls ++dL
CP001578_1 118 EVSHRSKNLLAIVQSIATQTARFTDSIDDFLVKFRGRIQSLSYSQDL 164
VpLtrseWagAdLreLleaeLaPyggedaeRitlsGPd<-*
+t+s+W+gA + +L++++++ y + d+eR+t++G++
CP001578_1 165 --VTDSNWRGALFLDLVRSQVEKYIDIDDERLTIEGDN 200
Histogram of all scores:
score obs exp (one = represents 7 sequences)
----- --- ---
-15 1 0|=
-14 16 0|===
-13 69 0|==========
-12 256 0|=====================================
-11 242 17|==*================================
-10 320 219|===============================*==============
-9 317 542|============================================== *
-8 374 563|====================================================== *
-7 225 378|================================= *
-6 105 206|=============== *
-5 81 102|============ *
-4 50 48|======*=
-3 30 22|===*=
-2 17 10|=*=
-1 10 4|*=
0 1 2|*
1 1 0|=
2 1 0|=
3 1 0|=
4 0 0|
> 5 3 -|=
% Statistical details of theoretical EVD fit:
mu = -7.9989
lambda = 0.7843
chi-sq statistic = 3240.7139
P(chi-square) = 0
Total sequences searched: 2120
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