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: CP000689.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
-------- ----------- ----- ------- ---
CP000689_0003 CDS_2918-4261 64.5 2.6e-19 1
CP000689_0059 CDS_75372-79709 47.2 2e-14 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000689_0003 1/1 223 289 .. 1 87 [] 64.5 2.6e-19
CP000689_0059 1/1 867 932 .. 1 87 [] 47.2 2e-14
Alignments of top-scoring domains:
CP000689_0003: domain 1 of 1, from 223 to 289: score 64.5, E = 2.6e-19
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++s+f a+l+HELRTP++ +rg +e++l ++s+++
CP000689_0 223 QLSQFSADLAHELRTPISNLRGELEVALS---------RPRSDDT-- 258
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
yre l ++l+e r++r i+ lL l+r+e
CP000689_0 259 ---------YREVLGSALEECTRISRVIQSLLFLARAEMV 289
CP000689_0059: domain 1 of 1, from 867 to 932: score 47.2, E = 2e-14
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
aks+fla++sH+LR+PL +I ++ elL e t+ +
CP000689_0 867 AKSRFLAAVSHDLRQPLQTIALLQELLTE---------TTSGKAA-- 902
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++ i++ + + ++n lLd+ +ieag
CP000689_0 903 ----------LGLIQQIKQTTSVMAGMLNTLLDINQIEAG 932
Histogram of all scores:
score obs exp (one = represents 1 sequences)
----- --- ---
-13 1 0|=
-12 6 3|==*===
-11 10 20|========== *
-10 39 37|====================================*==
-9 29 38|============================= *
-8 22 29|====================== *
-7 22 18|=================*====
-6 17 10|=========*=======
-5 10 5|====*=====
-4 1 3|= *
-3 8 1|*=======
-2 3 0|===
> -1 3 -|===
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 17.2262
P(chi-square) = 0.008487
Total sequences searched: 171
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
