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: AE008922.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
-------- ----------- ----- ------- ---
AE008922_2912 CDS_3521564-3520554 139.4 2.7e-47 1
AE008922_3637 CDS_4391309-4387752 130.4 3.1e-44 2
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
AE008922_2912 1/1 141 223 .. 1 85 [] 139.4 2.7e-47
AE008922_3637 2/2 989 1072 .. 1 85 [] 129.5 6.4e-44
Alignments of top-scoring domains:
AE008922_2912: domain 1 of 1, from 141 to 223: score 139.4, E = 2.7e-47
*->ELnHRVKNtLAvVQSiarQTlRnaasldefveaFsgRLqALsraHdL
ELnHRVKN+L +VQS+arQ+l+na+sld+++e++++RL+ALs+aH+
AE008922_2 141 ELNHRVKNSLVMVQSLARQSLKNAESLDDANEKIDARLMALSSAHNT 187
VpLtrseWagAdLreLleaeLaPyggedaeRitlsGPd<-*
Ltr++W +Ad++eL + +++ +++++++R++l+G +
AE008922_2 188 --LTRENWVSADIAELTRDAVELCESHQGQRFELHGSS 223
AE008922_3637: domain 2 of 2, from 989 to 1072: score 129.5, E = 6.4e-44
*->ELnHRVKNtLAvVQSiarQTlRnaasldefveaFsgRLqALsraHdL
EL+HR +N+++vVQS+a+QT ++++++def+++F++RL+ALsr+++L
AE008922_3 989 ELQHRTRNLITVVQSVATQTMQSSGTMDEFAGKFEDRLAALSRVQGL 1035
VpLtrseWagAdLreLleaeLaPygge.daeRitlsGPd<-*
L+++e a+A++++Ll++eL++ g+e+d++Ri l+GPd
AE008922_3 1036 --LSQAENAPATIGALLRMELDALGAEpDQRRIILDGPD 1072
Histogram of all scores:
score obs exp (one = represents 13 sequences)
----- --- ---
-14 10 0|=
-13 45 0|====
-12 279 0|======================
-11 429 34|==*==============================
-10 563 432|=================================*==========
-9 736 1070|========================================================= *
-8 726 1110|======================================================== *
-7 564 746|============================================ *
-6 297 407|======================= *
-5 215 201|===============*=
-4 139 95|=======*===
-3 78 44|===*==
-2 52 20|=*==
-1 19 9|*=
0 11 4|*
1 9 1|*
2 3 0|=
3 3 0|=
> 4 3 -|=
% Statistical details of theoretical EVD fit:
mu = -7.9989
lambda = 0.7843
chi-sq statistic = 5018.6582
P(chi-square) = 0
Total sequences searched: 4181
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
