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: CP001849.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
-------- ----------- ----- ------- ---
CP001849_0674 CDS_816767-818713 71.7 1.7e-20 1
CP001849_0557 CDS_652924-651749 62.9 5.5e-18 1
CP001849_0144 CDS_177004-175667 57.1 2.3e-16 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001849_0674 1/1 327 396 .. 1 87 [] 71.7 1.7e-20
CP001849_0557 1/1 163 227 .. 1 87 [] 62.9 5.5e-18
CP001849_0144 1/1 216 282 .. 1 87 [] 57.1 2.3e-16
Alignments of top-scoring domains:
CP001849_0674: domain 1 of 1, from 327 to 396: score 71.7, E = 1.7e-20
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ +f++++sHELRTPL+aI gyael + + +++ eq
CP001849_0 327 KMKRFVSDASHELRTPLAAIHGYAELYTM--------WRNMPGEQ-- 363
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++ + + i+ ++er+ +L++dlL l+r ++g
CP001849_0 364 ERA-------DDSIAHIKASSERMTELVEDLLSLARFDEG 396
CP001849_0557: domain 1 of 1, from 163 to 227: score 62.9, E = 5.5e-18
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++e++an+sHELRTP++a+++++e L + + ++
CP001849_0 163 MRRELIANVSHELRTPVSALQAMLENLAD---------GVVDPTP-- 198
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
e il++++rl +Li lLdlsr eag
CP001849_0 199 -----------SNMEAILNQTHRLADLITFLLDLSRVEAG 227
CP001849_0144: domain 1 of 1, from 216 to 282: score 57.1, E = 2.3e-16
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a+ ++a++sH++RTPL+ ++gya L++ + + ee+
CP001849_0 216 ARANWIAGVSHDIRTPLSMVMGYAGQLEN--------SLHILEED-- 252
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
r+ i++++ r+ +Lindl s++e +
CP001849_0 253 ----------RKKATMIIKQSVRMKNLINDLNLASKLEYN 282
Histogram of all scores:
score obs exp (one = represents 4 sequences)
----- --- ---
-15 1 0|=
-14 8 0|==
-13 21 0|======
-12 60 28|======*========
-11 121 149|=============================== *
-10 186 276|=============================================== *
-9 222 285|======================================================== *
-8 216 213|=====================================================*
-7 186 134|=================================*=============
-6 94 77|===================*====
-5 60 42|==========*====
-4 33 22|=====*===
-3 21 12|==*===
-2 12 6|=*=
-1 7 3|*=
0 4 1|*
1 1 0|=
2 0 0|
3 1 0|=
4 0 0|
5 0 0|
> 6 4 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 129.6540
P(chi-square) = 5.496e-23
Total sequences searched: 1258
Whole sequence top hits:
tophits_s report:
Total hits: 3
Satisfying E cutoff: 3
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 3
Satisfying E cutoff: 3
Total memory: 22K
