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: CP000644.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
-------- ----------- ----- ------- ---
CP000644_2576 CDS_2971287-2972999 78.9 5.6e-23 1
CP000644_1404 CDS_1607189-1605732 36.0 2.5e-10 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000644_2576 1/1 367 435 .. 1 70 [] 78.9 5.6e-23
CP000644_1404 1/1 306 372 .. 2 70 .] 36.0 2.5e-10
Alignments of top-scoring domains:
CP000644_2576: domain 1 of 1, from 367 to 435: score 78.9, E = 5.6e-23
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpe.earealde
ERa IARELHDs++Q+Ls +k+q+ l++r + +k + +a++a+++
CP000644_2 367 ERAVIARELHDSLAQALSYLKIQSTLLKRSF--AKGDNaKAEAAMQQ 411
irelarealaevRrllgdLRpaal<-*
i e ++ a++++R+llg++R ++
CP000644_2 412 IDEGLSSAYTQLRELLGTFRLSIG 435
CP000644_1404: domain 1 of 1, from 306 to 372: score 36.0, E = 2.5e-10
*->RaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdeir
R+ ARELHD +gQ+++ai+ q++ +rrl+ d++ + + + i+
CP000644_1 306 RKEVARELHDDIGQTITAIRTQAGIVRRLA---PDNAGVGQSSALIE 349
elarealaevRrllgdLR.paal<-*
+l+ ++ vR ll++L +
CP000644_1 350 TLSLGIYDAVRGLLRELElERHG 372
Histogram of all scores:
score obs exp (one = represents 14 sequences)
----- --- ---
-15 1 0|=
-14 16 0|==
-13 60 0|=====
-12 120 45|===*=====
-11 203 374|=============== *
-10 562 868|========================================= *
-9 784 986|======================================================== *
-8 607 760|============================================ *
-7 513 478|==================================*==
-6 372 270|===================*=======
-5 382 145|==========*=================
-4 193 75|=====*========
-3 98 38|==*====
-2 63 19|=*===
-1 45 10|*===
0 27 5|*=
1 16 2|*=
2 10 1|*
3 3 0|=
4 3 0|=
5 2 0|=
6 1 0|=
> 7 5 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 1384.5400
P(chi-square) = 0
Total sequences searched: 4086
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