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: CP001349.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
-------- ----------- ----- ------- ---
CP001349_6949 CDS_7357447-7358808 44.0 2.1e-12 1
CP001349_7209 CDS_7620007-7621440 40.7 1.9e-11 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001349_6949 1/1 247 314 .. 1 70 [] 44.0 2.1e-12
CP001349_7209 1/1 264 335 .. 1 70 [] 40.7 1.9e-11
Alignments of top-scoring domains:
CP001349_6949: domain 1 of 1, from 247 to 314: score 44.0, E = 2.1e-12
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ER+ +AR+LHD+ gQ L+a +e++++ + ++ + +r+ + +i
CP001349_6 247 ERRALARDLHDEFGQCLTATGALAETIAAEA--ADGRPDLRADALAI 291
relarealaevRrllgdLRpaal<-*
+ a + ++ +R +l +L p++
CP001349_6 292 GRNAARMMETLRGALARLEPPEI 314
CP001349_7209: domain 1 of 1, from 264 to 335: score 40.7, E = 1.9e-11
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpe....earea
ER+R A ELHD+ g L+a+ + ++r++ + +e+ +++++++
CP001349_7 264 ERRRTALELHDEFGPCLFALEANSASIARIA--GTLDEpargRLQQR 308
ldeirelarealaevRrllgdLRpaal<-*
d i++++++ ++ R+ll++LRpa+l
CP001349_7 309 TDDIAAIVKRVQTINRALLSRLRPAGL 335
Histogram of all scores:
score obs exp (one = represents 20 sequences)
----- --- ---
-15 2 0|=
-14 28 0|==
-13 104 1|*=====
-12 165 82|====*====
-11 284 674|=============== *
-10 780 1564|======================================= *
-9 1145 1776|==========================================================*
-8 981 1369|================================================== *
-7 1032 860|==========================================*=========
-6 820 487|========================*================
-5 803 261|=============*===========================
-4 418 136|======*==============
-3 303 70|===*============
-2 174 35|=*=======
-1 121 18|*======
0 73 9|*===
1 46 4|*==
2 25 2|*=
3 21 1|*=
4 8 0|=
5 7 0|=
6 3 0|=
7 3 0|=
8 4 0|=
9 0 0|
10 0 0|
11 1 0|=
12 0 0|
13 1 0|=
14 1 0|=
15 0 0|
16 0 0|
> 17 3 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 5324.5835
P(chi-square) = 0
Total sequences searched: 7356
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
