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: CP001584.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
-------- ----------- ----- ------- ---
CP001584_0228 CDS_288191-289684 86.4 8.9e-25 1
CP001584_0484 CDS_576567-579332 72.5 8e-21 1
CP001584_0442 CDS_524590-523277 64.2 1.8e-18 1
CP001584_0640 CDS_776847-775048 54.0 1.4e-15 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001584_0228 1/1 268 333 .. 1 87 [] 86.4 8.9e-25
CP001584_0484 1/1 695 763 .. 1 87 [] 72.5 8e-21
CP001584_0442 1/1 230 288 .. 1 87 [] 64.2 1.8e-18
CP001584_0640 1/1 371 441 .. 1 87 [] 54.0 1.4e-15
Alignments of top-scoring domains:
CP001584_0228: domain 1 of 1, from 268 to 333: score 86.4, E = 8.9e-25
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ ++f+an+sHELRTPL+aI+g++e+L + +e+
CP001584_0 268 SNTKFFANISHELRTPLNAIIGFSEILMS---------ERDTEKS-- 303
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ y++ i++++++ll++ind+Ldls++ a
CP001584_0 304 ----------KNYIKDIHDAGKHLLSMINDILDLSKASAD 333
CP001584_0484: domain 1 of 1, from 695 to 763: score 72.5, E = 8e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
k efl+n+ HE++TPLt+I+ + e+L +++ ++ +e+qr
CP001584_0 695 LKHEFLRNINHEINTPLTGIISLGETLWANY-------DKFNEDQR- 733
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
r ++e+i +++ rl++ in++Ld+s++ +
CP001584_0 734 ----------RNAVEVIAKSSIRLNSFINNILDFSKLSSL 763
CP001584_0442: domain 1 of 1, from 230 to 288: score 64.2, E = 1.8e-18
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ +la++sH+LRTPLt +++ +el++e
CP001584_0 230 KRTNMLAMISHDLRTPLTRMKLQLELIEE------------------ 258
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
e+++ +++++ ++++in Ld+++ e+
CP001584_0 259 ----------SEAIKGLKQDILTMQQMINSYLDFAKGEST 288
CP001584_0640: domain 1 of 1, from 371 to 441: score 54.0, E = 1.4e-15
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a s+ + ++HE+++PLt I++++e Ll+++ ++e +e+
CP001584_0 371 AWSDVAKKVAHEIKNPLTPILLASERLLKKF------SSEINEKS-- 409
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
e++ yl+ i+r ++ + ++++++ ++r+ a+
CP001584_0 410 --------EFESYLKMIIRHTNDIKNIVSEFVLFARLPAP 441
Histogram of all scores:
score obs exp (one = represents 4 sequences)
----- --- ---
-16 1 0|=
-15 5 0|==
-14 23 0|======
-13 36 0|=========
-12 60 21|=====*=========
-11 115 112|===========================*=
-10 185 209|=============================================== *
-9 166 216|========================================== *
-8 109 161|============================ *
-7 94 102|======================== *
-6 68 58|==============*==
-5 47 32|=======*====
-4 12 17|=== *
-3 11 9|==*
-2 11 4|*==
-1 1 2|*
0 1 1|*
1 1 0|=
2 0 0|
> 3 6 -|==
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 109.8807
P(chi-square) = 1.555e-19
Total sequences searched: 952
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
