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: CP000655.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
-------- ----------- ----- ------- ---
CP000655_1588 CDS_1667598-1666300 79.5 1.8e-22 1
CP000655_1800 CDS_1887252-1889672 77.6 5.9e-22 1
CP000655_1658 CDS_1732529-1729836 63.7 5.2e-18 1
CP000655_1246 CDS_1320410-1319277 61.9 1.7e-17 1
CP000655_0561 CDS_554638-552863 59.2 9.7e-17 1
CP000655_2064 CDS_2145368-2143056 57.1 3.8e-16 1
CP000655_0729 CDS_734283-731746 30.7 1.2e-08 1
CP000655_1079 CDS_1090996-1089554 27.7 8.5e-08 2
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000655_1588 1/1 206 272 .. 1 87 [] 79.5 1.8e-22
CP000655_1800 1/1 308 374 .. 1 87 [] 77.6 5.9e-22
CP000655_1658 1/1 666 733 .. 1 87 [] 63.7 5.2e-18
CP000655_1246 1/1 147 212 .. 1 87 [] 61.9 1.7e-17
CP000655_0561 1/1 218 283 .. 1 87 [] 59.2 9.7e-17
CP000655_2064 1/1 535 607 .. 1 87 [] 57.1 3.8e-16
CP000655_0729 1/1 599 671 .. 1 87 [] 30.7 1.2e-08
CP000655_1079 1/2 257 327 .. 1 87 [] 27.4 1e-07
Alignments of top-scoring domains:
CP000655_1588: domain 1 of 1, from 206 to 272: score 79.5, E = 1.8e-22
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++++f+an+sHE+RTP+t+++g++e+ +l eq+
CP000655_1 206 MRRDFVANVSHEMRTPITVLMGFLETVHS---------LDLKKEQQ- 242
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+y++ + +a+r+ +L++dlL+l+ +e++
CP000655_1 243 ----------DQYFDMMMAQAQRMKSLVEDLLTLANLESN 272
CP000655_1800: domain 1 of 1, from 308 to 374: score 77.6, E = 5.9e-22
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
aks fla++sHE+RTP++aI+g elLl+ ++ + +
CP000655_1 308 AKSAFLATMSHEIRTPMNAILGVQELLLN--------SPQFPNRE-- 344
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++l++++ +ae+ll ++n++Ldls+ieag
CP000655_1 345 ----------KRLLKSAHASAESLLGILNQVLDLSKIEAG 374
CP000655_1658: domain 1 of 1, from 666 to 733: score 63.7, E = 5.2e-18
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ +l+++sH+LRTPLta+ g++++L++ d++l++eq
CP000655_1 666 LRNSLLSAISHDLRTPLTALVGLTDALEM-------LDAPLTSEQ-- 703
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+e+ +++ a +++++++lLd++r+++g
CP000655_1 704 ----------KEIAVLLREKALKMSSQVSNLLDMARLQSG 733
CP000655_1246: domain 1 of 1, from 147 to 212: score 61.9, E = 1.7e-17
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a ++++nl+HE+++PL +Irg+a+lL+ else+ +
CP000655_1 147 ANKQLMRNLAHEIKNPLGGIRGAAQLLEF----------ELSEKGL- 182
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ey+++i++e++rl++L+++lL r +
CP000655_1 183 ----------HEYTQVIIKESDRLQSLVDRLLAPHRKAHV 212
CP000655_0561: domain 1 of 1, from 218 to 283: score 59.2, E = 9.7e-17
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
aks+f+a++sH+LR+PL+a+ y++l + +l+ ++
CP000655_0 218 AKSKFIATASHDLRQPLHAVNIYLDLFEP---------NKLNQKE-- 253
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
r + i+++++ l+++ lLdls+++ag
CP000655_0 254 ----------RINFLQIKKSIQTLNSMFKSLLDLSKLDAG 283
CP000655_2064: domain 1 of 1, from 535 to 607: score 57.1, E = 3.8e-16
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a se ++ l+HE+++PLt I++ ae L+ +l ++l++eq
CP000655_2 535 AWSEVARRLAHEIKNPLTPIQLSAERLQHKLA------GKLTPEQ-- 573
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++ + +eti+ +++ + +++nd++d+++ ++
CP000655_2 574 ------EEMINRSTETIIGQVQAMKEMVNDFRDFAKTPTP 607
CP000655_0729: domain 1 of 1, from 599 to 671: score 30.7, E = 1.2e-08
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++e++++l+HEL++PL+aI+ y+ ++l+ ++ ++
CP000655_0 599 TMGEMASSLAHELNQPLSAISNYCMGVAKRLE------GNIDPALGK 639
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ L +le++ ++a+r +i++++ + + ++
CP000655_0 640 DIL--------PALEKASDQAHRAGMIIQRIRGFVKRSEP 671
CP000655_1079: domain 1 of 2, from 257 to 327: score 27.4, E = 1e-07
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ + + a+++HEL++PL a +++++ L+++l d++l++ +
CP000655_1 257 STGALSASIAHELNQPLGASSLNIQFLQKKLA-----DGDLDPSVQ- 297
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
e+l+t+l + +r ++i+ l+ + +
CP000655_1 298 ----------SEILNTLLVDNNRAANIIRSLRSIFSEKHV 327
Histogram of all scores:
score obs exp (one = represents 8 sequences)
----- --- ---
-14 3 0|=
-13 21 1|*==
-12 86 47|=====*=====
-11 217 246|============================ *
-10 373 456|=============================================== *
-9 421 471|===================================================== *
-8 374 353|============================================*==
-7 241 222|===========================*===
-6 141 128|===============*==
-5 83 70|========*==
-4 50 37|====*==
-3 31 19|==*=
-2 15 10|=*
-1 7 5|*
0 2 2|*
> 1 12 -|==
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 74.6964
P(chi-square) = 1.55e-11
Total sequences searched: 2077
Whole sequence top hits:
tophits_s report:
Total hits: 8
Satisfying E cutoff: 8
Total memory: 21K
Domain top hits:
tophits_s report:
Total hits: 8
Satisfying E cutoff: 8
Total memory: 24K