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/HisKA_3.hmm [HisKA_3]
Sequence database: AE008923.Gprot.fsa
per-sequence score cutoff: >= 30.1 [TC1]
per-domain score cutoff: >= 30.1 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HisKA_3
Accession: PF07730.2
Description: 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
-------- ----------- ----- ------- ---
AE008923_3938 CDS_4686511-4685234 88.9 6.8e-26 1
AE008923_1208 CDS_1406493-1407689 72.3 5.4e-21 1
AE008923_0135 CDS_158735-159961 64.3 1.3e-18 1
AE008923_1963 CDS_2327087-2324955 62.6 3.9e-18 1
AE008923_2128 CDS_2536698-2539748 50.9 1.1e-14 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
AE008923_3938 1/1 223 290 .. 1 70 [] 88.9 6.8e-26
AE008923_1208 1/1 199 263 .. 1 70 [] 72.3 5.4e-21
AE008923_0135 1/1 223 288 .. 1 70 [] 64.3 1.3e-18
AE008923_1963 1/1 504 569 .. 1 70 [] 62.6 3.9e-18
AE008923_2128 1/1 792 859 .. 1 70 [] 50.9 1.1e-14
Alignments of top-scoring domains:
AE008923_3938: domain 1 of 1, from 223 to 290: score 88.9, E = 6.8e-26
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ER+RI R+LHD +g++Ls+++l+++la+rl+ ++dp++ar ++d++
AE008923_3 223 ERERIGRDLHDLLGHTLSLVALKSDLAARLI--ERDPQAARNEIDAV 267
relarealaevRrllgdLRpaal<-*
++ +reala+vRr+++++R a+l
AE008923_3 268 AQVSREALAQVRRAVTGIRAAGL 290
AE008923_1208: domain 1 of 1, from 199 to 263: score 72.3, E = 5.4e-21
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ER+RI+RELHD +g++L+a++l le+a +++ ++a+e+++++
AE008923_1 199 ERTRISRELHDLLGHHLTALSLNLEVAGHIT-----EGQAQEHVRQA 240
relarealaevRrllgdLRpaal<-*
++la+ l++vR+++++LR+++
AE008923_1 241 HTLAKLLLTDVREAVSHLRDSGA 263
AE008923_0135: domain 1 of 1, from 223 to 288: score 64.3, E = 1.3e-18
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ER R+ARELHD+ g+ L+a+++ l l+ ++++ + r++++ +
AE008923_0 223 ERLRLARELHDVAGHKLTAMRINLRLLS-----ADPALAQRDEVAVV 264
relarealaevRrllgdLR.paal<-*
++l+ e l+++R++++ LR++++l
AE008923_0 265 EQLSAELLSDIRNVVQSLRdDQGL 288
AE008923_1963: domain 1 of 1, from 504 to 569: score 62.6, E = 3.9e-18
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ER+ I+R+LHD +gQ+++a+kl +++a+ d + d+++ re l e+
AE008923_1 504 ERRAISRDLHDDIGQAITAMKLSAHAAL---D-EPDADARREDLLEV 546
relarealaevRrllgdLRpaal<-*
la+++ + +R+l LRp++l
AE008923_1 547 VSLADTTVTKLRNLSMLLRPPQL 569
AE008923_2128: domain 1 of 1, from 792 to 859: score 50.9, E = 1.1e-14
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ER+RIARELHD++ Q+ +++l l++a r l + ++ + r +l+++
AE008923_2 792 ERERIARELHDTLLQGTQGLILRLHAASRSL--PWND-PRRLELEQA 835
relarealaevRrllgdLR.paal<-*
+la+ alae R+ +++LR+ + l
AE008923_2 836 VDLAENALAEGRDRVSGLRdSQSL 859
Histogram of all scores:
score obs exp (one = represents 12 sequences)
----- --- ---
-14 6 0|=
-13 36 0|===
-12 65 48|===*==
-11 124 395|=========== *
-10 425 917|==================================== *
-9 697 1041|==========================================================*
-8 649 802|======================================================= *
-7 662 504|=========================================*==============
-6 530 285|=======================*=====================
-5 482 153|============*============================
-4 219 79|======*============
-3 154 41|===*=========
-2 90 21|=*======
-1 64 10|*=====
0 37 5|*===
1 24 2|*=
2 14 1|*=
3 8 0|=
4 10 0|=
5 5 0|=
6 2 0|=
7 2 0|=
> 8 7 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 2790.4707
P(chi-square) = 0
Total sequences searched: 4312
Whole sequence top hits:
tophits_s report:
Total hits: 5
Satisfying E cutoff: 5
Total memory: 21K
Domain top hits:
tophits_s report:
Total hits: 5
Satisfying E cutoff: 5
Total memory: 22K
