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: CP000679.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
-------- ----------- ----- ------- ---
CP000679_2466 CDS_2723434-2721737 95.6 6.3e-27 1
CP000679_1953 CDS_2177764-2176067 90.1 2.3e-25 1
CP000679_1922 CDS_2150980-2149577 75.4 3.3e-21 1
CP000679_2206 CDS_2410812-2409445 74.6 5.7e-21 1
CP000679_0393 CDS_465628-466881 57.8 3.2e-16 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000679_2466 1/1 342 407 .. 1 87 [] 95.6 6.3e-27
CP000679_1953 1/1 345 412 .. 1 87 [] 90.1 2.3e-25
CP000679_1922 1/1 246 312 .. 1 87 [] 75.4 3.3e-21
CP000679_2206 1/1 238 304 .. 1 87 [] 74.6 5.7e-21
CP000679_0393 1/1 210 277 .. 1 87 [] 57.8 3.2e-16
Alignments of top-scoring domains:
CP000679_2466: domain 1 of 1, from 342 to 407: score 95.6, E = 6.3e-27
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ +f+an+sHELRTP+t+I++y+e+Ll+ +ee++
CP000679_2 342 MRKQFVANVSHELRTPITTIKTYSETLLD----------VDNEETK- 377
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++l++i++e +r+ rL++dlL lsr+++g
CP000679_2 378 ----------KQFLSVIIKECDRMTRLVSDLLYLSRLDSG 407
CP000679_1953: domain 1 of 1, from 345 to 412: score 90.1, E = 2.3e-25
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+s+f+an+sHEL+TPLt+I+g++e+L++ ++ +++++
CP000679_1 345 IRSDFVANVSHELKTPLTSIKGFVETLKD---------GAIEDKEIA 382
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++l++i+ e+erl rLindlL ls ie+
CP000679_1 383 ----------GKFLNIIEVEVERLVRLINDLLYLSEIENA 412
CP000679_1922: domain 1 of 1, from 246 to 312: score 75.4, E = 3.3e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ f++n+sHELR PLt+Irg++e +l+ + ++e+r
CP000679_1 246 MRKSFISNVSHELRSPLTSIRGFIEGILD---------RTIPDEKR- 282
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
yl+ +++e+ l++Lin+lL+lsr+e g
CP000679_1 283 ----------DFYLNLVREEVIKLNNLINQLLELSRLEWG 312
CP000679_2206: domain 1 of 1, from 238 to 304: score 74.6, E = 5.7e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a+ +fl n+sHEL+TPLt+Irgyae+L++ + l+ e++
CP000679_2 238 AQIRFLQNISHELKTPLTSIRGYAEVLKM---------GLLDVEKA- 274
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ ++i++ ++ l Lin++ dl++ie+
CP000679_2 275 ----------NHAANRIIEHVDKLKILINQIIDLTKIESV 304
CP000679_0393: domain 1 of 1, from 210 to 277: score 57.8, E = 3.2e-16
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ef+an H+L+TP+++I+gy+e + + +++++
CP000679_0 210 VRREFVANFVHDLKTPIAVIKGYCESTKY---------LDIDDKEK- 246
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++ + i+++ +++++Li+++ +ls+i ag
CP000679_0 247 ---------IKQNMDGIEKQCDYMQNLISNMVELSKISAG 277
Histogram of all scores:
score obs exp (one = represents 9 sequences)
----- --- ---
-15 1 0|=
-14 2 0|=
-13 31 2|*===
-12 136 61|======*=========
-11 312 317|===================================*
-10 466 589|==================================================== *
-9 498 608|======================================================== *
-8 418 455|=============================================== *
-7 317 287|===============================*====
-6 220 165|==================*======
-5 113 90|=========*===
-4 68 48|=====*==
-3 33 25|==*=
-2 21 13|=*=
-1 16 7|*=
0 5 3|*
1 6 1|*
2 5 0|=
3 2 0|=
> 4 9 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 192.7459
P(chi-square) = 2.373e-35
Total sequences searched: 2679
Whole sequence top hits:
tophits_s report:
Total hits: 5
Satisfying E cutoff: 5
Total memory: 21K
Domain top hits:
tophits_s report:
Total hits: 5
Satisfying E cutoff: 5
Total memory: 22K
