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: CP000531.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
-------- ----------- ----- ------- ---
CP000531_0113 CDS_123563-122187 51.8 9.6e-16 1
CP000531_0034 CDS_41004-39859 42.3 4.8e-13 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000531_0113 1/1 234 300 .. 1 87 [] 51.8 9.6e-16
CP000531_0034 1/1 159 216 .. 1 80 [. 42.3 4.8e-13
Alignments of top-scoring domains:
CP000531_0113: domain 1 of 1, from 234 to 300: score 51.8, E = 9.6e-16
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++ f l+HE+RTP+++++g ++++l ++++ ++
CP000531_0 234 RLQAFSGELAHEMRTPVATLMGRTQVALS---------KTRTAAEL- 270
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
re le +e+erl rLi d+L ++ +e++
CP000531_0 271 ----------REVLEGNVEELERLTRLIADMLFIASAEHN 300
CP000531_0034: domain 1 of 1, from 159 to 216: score 42.3, E = 4.8e-13
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++ +a+l H++RTPL++I g+a+lL + + eq+
CP000531_0 159 FRERIAATLTHDMRTPLSVIVGAAHLLTR-----------VDKEQI- 193
peLRledpeyreyletilreaerllrLindlLd<-*
+ + ++i++ ++rl+++ n+ Ld
CP000531_0 194 ----------ATLAKKIIDNGRRLEAMFNEQLD 216
Histogram of all scores:
score obs exp (one = represents 1 sequences)
----- --- ---
-13 6 0|======
-12 3 3|==*
-11 20 19|==================*=
-10 28 35|============================ *
-9 34 36|================================== *
-8 22 27|====================== *
-7 27 17|================*==========
-6 8 9|========*
-5 3 5|=== *
-4 4 2|=*==
-3 1 1|*
-2 1 0|=
> -1 4 -|====
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 8.7088
P(chi-square) = 0.1213
Total sequences searched: 161
Whole sequence top hits:
tophits_s report:
Total hits: 2
Satisfying E cutoff: 2
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 2
Satisfying E cutoff: 2
Total memory: 21K
