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: CP000057.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
-------- ----------- ----- ------- ---
CP000057_0276 CDS_303809-301968 74.8 3.2e-21 1
CP000057_1555 CDS_1620401-1621678 70.8 4.6e-20 1
CP000057_1760 CDS_1870740-1869385 46.3 3.8e-13 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000057_0276 1/1 136 201 .. 1 87 [] 74.8 3.2e-21
CP000057_1555 1/1 194 260 .. 1 87 [] 70.8 4.6e-20
CP000057_1760 1/1 234 300 .. 1 87 [] 46.3 3.8e-13
Alignments of top-scoring domains:
CP000057_0276: domain 1 of 1, from 136 to 201: score 74.8, E = 3.2e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
ks +la++sHE+RTPL++I g++++Ll+ +el++ q
CP000057_0 136 DKSTLLATISHEFRTPLNGIVGLSQILLD---------DELDDLQ-- 171
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
r yl+ti +a +l + +d+ dl++i+a
CP000057_0 172 ----------RNYLKTINISAVSLGYIFSDIIDLEKIDAS 201
CP000057_1555: domain 1 of 1, from 194 to 260: score 70.8, E = 4.6e-20
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++++fl+n+ HELRTPLt+++gy+e+L ++ +++ q
CP000057_1 194 SRQKFLSNINHELRTPLTVLQGYLEILADN--------NIQNPLQ-- 230
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++++ ++++++r+++L++++ l++ie+
CP000057_1 231 ----------KKAIFAMQEQSQRMEHLLQQFNFLAKIETT 260
CP000057_1760: domain 1 of 1, from 234 to 300: score 46.3, E = 3.8e-13
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ ++f+++++HELR PL+a+r +e+++ + + r
CP000057_1 234 RERRFTSDAAHELRSPLAALRIQTEVAQL---------AGDDVALR- 270
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++l + + ++r ++Li++lL+lsr+++
CP000057_1 271 ----------EQALLHLTQGIDRASQLIEQLLTLSRLDNL 300
Histogram of all scores:
score obs exp (one = represents 6 sequences)
----- --- ---
-14 6 0|=
-13 21 1|*===
-12 93 40|======*=========
-11 185 212|=============================== *
-10 334 393|======================================================== *
-9 332 406|======================================================== *
-8 294 304|================================================= *
-7 204 191|===============================*==
-6 143 110|==================*=====
-5 73 60|=========*===
-4 46 32|=====*==
-3 27 17|==*==
-2 16 8|=*=
-1 6 4|*
0 2 2|*
1 2 1|*
2 2 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 = 122.6741
P(chi-square) = 1.449e-21
Total sequences searched: 1792
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
