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: CP000746.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
-------- ----------- ----- ------- ---
CP000746_0848 CDS_924497-926332 84.0 9.2e-24 1
CP000746_0770 CDS_842506-843927 62.6 1.1e-17 1
CP000746_1342 CDS_1517925-1519298 57.0 4.1e-16 1
CP000746_1690 CDS_1888650-1890632 31.1 8.9e-09 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000746_0848 1/1 93 158 .. 1 87 [] 84.0 9.2e-24
CP000746_0770 1/1 250 312 .. 1 87 [] 62.6 1.1e-17
CP000746_1342 1/1 239 305 .. 1 87 [] 57.0 4.1e-16
CP000746_1690 1/1 444 512 .. 1 87 [] 31.1 8.9e-09
Alignments of top-scoring domains:
CP000746_0848: domain 1 of 1, from 93 to 158: score 84.0, E = 9.2e-24
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
k++++a++sHELRTPL++I g++ +Ll+ t+l+eeq
CP000746_0 93 DKTDLMATISHELRTPLNGIVGLTRILLD---------TPLTEEQ-- 128
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ yl+ti +a +l ++ nd+ dl++i++
CP000746_0 129 ----------KNYLNTINVSAVSLGNIFNDVIDLEKIDSR 158
CP000746_0770: domain 1 of 1, from 250 to 312: score 62.6, E = 1.1e-17
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++++++sHELRTPLt +++++ lL++++ +e +e
CP000746_0 250 THQRLVSSISHELRTPLTRLQLAVALLRRRV-------GEYPE---- 285
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++i+ e+erl+++indlL lsr +
CP000746_0 286 -------------TKRIETESERLDKMINDLLLLSRNQLR 312
CP000746_1342: domain 1 of 1, from 239 to 305: score 57.0, E = 4.1e-16
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ ++f+++++HELR PL+++r ael+++ + + + r
CP000746_1 239 RERRFTSDAAHELRSPLAGLRIQAELAQM---------AGDDVKLR- 275
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+e+l+ + + ++r +Li++lL+lsr+++
CP000746_1 276 ----------EEALNHLTDGIDRATQLIEQLLTLSRLDNL 305
CP000746_1690: domain 1 of 1, from 444 to 512: score 31.1, E = 8.9e-09
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ l ++sHE+ +PL+a+ +y+ +++ l ++++ e +
CP000746_1 444 VVGKTLTSISHEITQPLNAMNAYIFSAKRAL-------QKQNYEAT- 482
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+yl++i r +e +i +l+++s+ ++g
CP000746_1 483 ----------YDYLDKINRLTEKTAVIIKRLRQFSKQGNG 512
Histogram of all scores:
score obs exp (one = represents 7 sequences)
----- --- ---
-15 1 0|=
-14 1 0|=
-13 19 1|*==
-12 86 47|======*======
-11 187 246|=========================== *
-10 370 457|===================================================== *
-9 396 472|========================================================= *
-8 383 353|==================================================*====
-7 266 222|===============================*======
-6 178 128|==================*=======
-5 82 70|=========*==
-4 45 37|=====*=
-3 23 19|==*=
-2 12 10|=*
-1 13 5|*=
0 5 2|*
1 5 1|*
2 2 0|=
> 3 5 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 119.1486
P(chi-square) = 2.689e-20
Total sequences searched: 2079
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
