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: CP000687.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
-------- ----------- ----- ------- ---
CP000687_0897 CDS_1042196-1040346 81.4 5e-23 1
CP000687_1217 CDS_1395984-1394692 75.6 2.2e-21 1
CP000687_1030 CDS_1196930-1198306 55.6 1e-15 1
CP000687_0586 CDS_676847-675447 50.5 2.8e-14 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000687_0897 1/1 114 179 .. 1 87 [] 81.4 5e-23
CP000687_1217 1/1 203 268 .. 1 87 [] 75.6 2.2e-21
CP000687_1030 1/1 239 305 .. 1 87 [] 55.6 1e-15
CP000687_0586 1/1 245 307 .. 1 87 [] 50.5 2.8e-14
Alignments of top-scoring domains:
CP000687_0897: domain 1 of 1, from 114 to 179: score 81.4, E = 5e-23
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
ks ++a++sHELRTPL++I+g++ +Lle +el+e+q
CP000687_0 114 DKSVLMATISHELRTPLNGIIGLSRILLE---------GELTEQQ-- 149
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
reyl+ti +a +l ++ +d+ dl++i++
CP000687_0 150 ----------REYLKTINISAVSLGHIFSDIIDLEKIDSR 179
CP000687_1217: domain 1 of 1, from 203 to 268: score 75.6, E = 2.2e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++ fl n+ HELRTPLt++rgy+elL+ ++ +e q
CP000687_1 203 SRQMFLTNMNHELRTPLTVLRGYLELLEG---------QAETELQ-- 238
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++ ++ ++++a+r+ +L+++l l++ie+
CP000687_1 239 ----------QKSIHAMQSQAKRMGNLLDQLNLLAKIETS 268
CP000687_1030: domain 1 of 1, from 239 to 305: score 55.6, E = 1e-15
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ ++f+++++HELR PL+++r +e++++ t +++++
CP000687_1 239 RERRFTSDAAHELRSPLAGLRIQTEIAQM---------TLDDPDTH- 275
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+l+ + + ++r+ +Li++lL+lsr+e+
CP000687_1 276 ----------LSALNNLTHGIDRIAQLIEQLLTLSRLENL 305
CP000687_0586: domain 1 of 1, from 245 to 307: score 50.5, E = 2.8e-14
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++++++++sHEL+TPL+ + +++ l+++ + +++
CP000687_0 245 HQQRLVSDISHELKTPLARLHLATALIRR---------RNGESPE-- 280
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
l +i+++++ l+ ++ dlL lsr + +
CP000687_0 281 -------------LTRIENQIQKLDTMVHDLLVLSRQQVN 307
Histogram of all scores:
score obs exp (one = represents 7 sequences)
----- --- ---
-15 2 0|=
-14 5 0|=
-13 23 1|*===
-12 101 46|======*========
-11 205 241|============================== *
-10 339 447|================================================= *
-9 405 462|==========================================================*
-8 355 346|=================================================*=
-7 247 218|===============================*====
-6 146 125|=================*===
-5 85 68|=========*===
-4 51 36|=====*==
-3 35 19|==*==
-2 8 10|=*
-1 10 5|*=
0 10 2|*=
1 3 1|*
2 0 0|
3 0 0|
> 4 6 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 136.2437
P(chi-square) = 9.445e-24
Total sequences searched: 2036
Whole sequence top hits:
tophits_s report:
Total hits: 4
Satisfying E cutoff: 4
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 4
Satisfying E cutoff: 4
Total memory: 22K
