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_1.hmm [HisKA]
Sequence database: CP001751.Gprot.fsa
per-sequence score cutoff: >= 26.8 [TC1]
per-domain score cutoff: >= 26.8 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HisKA
Accession: PF00512.12
Description: His Kinase A (phosphoacceptor) domain
[HMM has been calibrated; E-values are empirical estimates]
Scores for complete sequences (score includes all domains):
Sequence Description Score E-value N
-------- ----------- ----- ------- ---
CP001751_1207 CDS_1294593-1293553 76.3 1.8e-21 1
CP001751_1266 CDS_1391583-1389955 60.4 5.6e-17 1
CP001751_2203 CDS_2411500-2410181 55.4 1.5e-15 1
CP001751_0217 CDS_229289-231478 47.7 2.2e-13 1
CP001751_0215 CDS_226612-227784 46.7 4.2e-13 1
CP001751_1735 CDS_1911748-1910375 40.7 2.2e-11 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001751_1207 1/1 117 184 .. 1 87 [] 76.3 1.8e-21
CP001751_1266 1/1 321 386 .. 1 87 [] 60.4 5.6e-17
CP001751_2203 1/1 233 291 .. 1 87 [] 55.4 1.5e-15
CP001751_0217 1/1 498 569 .. 1 87 [] 47.7 2.2e-13
CP001751_0215 1/1 151 206 .. 1 78 [. 46.7 4.2e-13
CP001751_1735 1/1 238 303 .. 1 87 [] 40.7 2.2e-11
Alignments of top-scoring domains:
CP001751_1207: domain 1 of 1, from 117 to 184: score 76.3, E = 1.8e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++f+an+sHELR PLt+++g++e++ + ++ +++ +
CP001751_1 117 VRRDFVANVSHELRSPLTSLSGFIETMIN---------GDVTDDAMK 154
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++l+++++e++r+ rLi+dlL lsr e
CP001751_1 155 ----------LRFLKIMDEESKRMTRLIDDLLSLSRVEVY 184
CP001751_1266: domain 1 of 1, from 321 to 386: score 60.4, E = 5.6e-17
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a f+a+++HEL++PL+++r +ae++ + ++ q+
CP001751_1 321 ATAGFAADVAHELKNPLSSLRSAAETISR----------ISDPVQQ- 356
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++ ++il+++erl+rLi d+ + sr +++
CP001751_1 357 ----------QKLMNVILQDVERLNRLITDISQASRVDTE 386
CP001751_2203: domain 1 of 1, from 233 to 291: score 55.4, E = 1.5e-15
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++ +la++sH+LRTPLt +r+ + l+++
CP001751_2 233 ERTAMLAGVSHDLRTPLTRVRLQLALMED------------------ 261
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ l+ i+++++ ++++i+ L ++ e +
CP001751_2 262 ----------NDDLKAIRQDIDEMEDMIDGYLSFAAGEGE 291
CP001751_0217: domain 1 of 1, from 498 to 569: score 47.7, E = 2.2e-13
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a + ++ ++HE+++PLt I+++++ L++++ ++ +e
CP001751_0 498 AWADIARRIAHEIKNPLTPIQLASDRLRRKF----KPKDAIEAER-- 538
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++ey+ +i r+++ + r+++++ ++r+ ++
CP001751_0 539 ---------FEEYVTIISRQVDDIGRMVDEFSSFARMPQP 569
CP001751_0215: domain 1 of 1, from 151 to 206: score 46.7, E = 4.2e-13
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ + ++a+l+HE ++PL++I+g+a+lL++ +l + +
CP001751_0 151 SVTSMAAMLAHEVKNPLAGIKGAAQLLEA----------DLEASE-- 185
peLRledpeyreyletilreaerllrLindl<-*
+++ i +e++r +L++++
CP001751_0 186 ----------TNLTRMIVEETDRVTALLDRM 206
CP001751_1735: domain 1 of 1, from 238 to 303: score 40.7, E = 2.2e-11
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+k++f++n++H+LR+PL++I+++ae++ + +++ ++
CP001751_1 238 QKDRFISNAAHQLRNPLASILALAEVATD----------AKTIADA- 273
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
r +l ++++ rL ++lL ++r++
CP001751_1 274 ----------RHRSGELLTASRQAVRLSEQLLSYERLRGA 303
Histogram of all scores:
score obs exp (one = represents 9 sequences)
----- --- ---
-15 1 0|=
-14 4 0|=
-13 23 1|*==
-12 114 58|======*======
-11 230 301|========================== *
-10 447 559|================================================== *
-9 522 578|==========================================================*
-8 423 432|===============================================*
-7 317 272|==============================*=====
-6 222 157|=================*=======
-5 99 86|=========*=
-4 58 46|=====*=
-3 37 24|==*==
-2 21 12|=*=
-1 9 6|*
0 4 3|*
1 4 1|*
2 1 0|=
> 3 10 -|==
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 150.7989
P(chi-square) = 1.023e-26
Total sequences searched: 2546
Whole sequence top hits:
tophits_s report:
Total hits: 6
Satisfying E cutoff: 6
Total memory: 21K
Domain top hits:
tophits_s report:
Total hits: 6
Satisfying E cutoff: 6
Total memory: 23K
