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_3.hmm [HisKA_3]
Sequence database: CP001892.Gprot.fsa
per-sequence score cutoff: >= 30.1 [TC1]
per-domain score cutoff: >= 30.1 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HisKA_3
Accession: PF07730.2
Description: Histidine kinase
[HMM has been calibrated; E-values are empirical estimates]
Scores for complete sequences (score includes all domains):
Sequence Description Score E-value N
-------- ----------- ----- ------- ---
CP001892_1018 CDS_1208982-1211537 84.5 5e-25 2
CP001892_1088 CDS_1301652-1302818 47.9 3.1e-14 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001892_1018 1/2 229 296 .. 1 70 [] 69.7 1.1e-20
CP001892_1088 1/1 191 260 .. 1 70 [] 47.9 3.1e-14
Alignments of top-scoring domains:
CP001892_1018: domain 1 of 2, from 229 to 296: score 69.7, E = 1.1e-20
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ERaRIAR++HD v+++Ls i++q+++ r+++ + d + ar ++d+i
CP001892_1 229 ERARIARDMHDIVAHTLSIIIIQSDGGRYAA--THDVQLARSTMDTI 273
relarealaevRrllgdLRpaal<-*
r+ a++a+ ++ +llg++ + +
CP001892_1 274 RHEAEHAMHDMKQLLGVFGGSSH 296
CP001892_1088: domain 1 of 1, from 191 to 260: score 47.9, E = 3.1e-14
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpe.earealde
ER+RIARE+HD vg+ L+ +++ + + + ++ ++ ++ ++++ +
CP001892_1 191 ERTRIAREIHDNVGHLLTRAIMLTNADQVIA--TSTGDtLHADQFGQ 235
irelarealaevRrllgdLR.paal<-*
+ ++++ea++ +R +++L++++++
CP001892_1 236 LGHTLDEAMTMIRKSVHGLKdEGTD 260
Histogram of all scores:
score obs exp (one = represents 5 sequences)
----- --- ---
-14 9 0|==
-13 16 0|====
-12 35 17|===*===
-11 67 144|============== *
-10 200 334|======================================== *
-9 286 379|==========================================================*
-8 246 292|================================================== *
-7 183 183|====================================*
-6 159 104|====================*===========
-5 152 55|==========*====================
-4 87 29|=====*============
-3 51 14|==*========
-2 37 7|=*======
-1 10 3|*=
0 7 1|*=
1 11 1|*==
2 5 0|=
3 1 0|=
4 3 0|=
5 2 0|=
6 0 0|
7 0 0|
> 8 5 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 651.2999
P(chi-square) = 0
Total sequences searched: 1572
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
