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: CP000301.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
-------- ----------- ----- ------- ---
CP000301_3978 CDS_4481127-4479454 127.7 2.9e-43 1
CP000301_0609 CDS_668794-668018 119.3 2.2e-40 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000301_3978 1/1 349 427 .. 1 85 [] 127.7 2.9e-43
CP000301_0609 1/1 68 150 .. 1 85 [] 119.3 2.2e-40
Alignments of top-scoring domains:
CP000301_3978: domain 1 of 1, from 349 to 427: score 127.7, E = 2.9e-43
*->ELnHRVKNtLAvVQSiarQTlRnaasldefveaFsgRLqALsraHdL
ELnHRVKN+LAv+QSi++ TlR+ ++d ++e ++gR++A+s++ +L
CP000301_3 349 ELNHRVKNILAVIQSIVTRTLRHGSDVDVARELLIGRIHAMSNVVSL 395
VpLtrseWagAdLreLleaeLaPyggedaeRitlsGPd<-*
Lt+s+W+g+ Lr Llea P+ +eRi GPd
CP000301_3 396 --LTESQWQGVKLRGLLEARAIPH----SERIAANGPD 427
CP000301_0609: domain 1 of 1, from 68 to 150: score 119.3, E = 2.2e-40
*->ELnHRVKNtLAvVQSiarQTlRnaasldefveaFsgRLqALsraHdL
EL+HR+KNtLA+V++ia Q+lRna+s +++ +a++gRL AL+raHdL
CP000301_0 68 ELHHRIKNTLATVSAIASQSLRNATSTEHAQQAIEGRLLALGRAHDL 114
VpLtrseWagAdLreLleaeLaPyggedaeRitlsGPd<-*
L +++W +AdL+++++ + ++++ da + + GPd
CP000301_0 115 --LLQAQWTSADLGKVIRGATEAFDNRDAAKFSIAGPD 150
Histogram of all scores:
score obs exp (one = represents 15 sequences)
----- --- ---
-16 1 0|=
-15 2 0|=
-14 4 0|=
-13 61 0|=====
-12 395 0|===========================
-11 445 39|==*===========================
-10 786 505|=================================*===================
-9 835 1250|======================================================== *
-8 876 1298|==========================================================*
-7 581 872|======================================= *
-6 355 476|======================== *
-5 232 235|===============*
-4 139 111|=======*==
-3 73 51|===*=
-2 43 23|=*=
-1 19 10|*=
0 12 5|*
1 10 2|*
2 4 1|*
3 2 0|=
4 1 0|=
5 2 0|=
6 1 0|=
> 7 7 -|=
% Statistical details of theoretical EVD fit:
mu = -7.9989
lambda = 0.7843
chi-sq statistic = 4714.9331
P(chi-square) = 0
Total sequences searched: 4886
Whole sequence top hits:
tophits_s report:
Total hits: 2
Satisfying E cutoff: 2
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 2
Satisfying E cutoff: 2
Total memory: 21K
