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: BA000032.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
-------- ----------- ----- ------- ---
BA000032_1196 CDS_1270699-1268963 72.2 2.3e-21 1
BA000032_0965 CDS_1007737-1009239 53.0 1.1e-15 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
BA000032_1196 1/1 369 441 .. 1 70 [] 72.2 2.3e-21
BA000032_0965 1/1 303 368 .. 1 70 [] 53.0 1.1e-15
Alignments of top-scoring domains:
BA000032_1196: domain 1 of 1, from 369 to 441: score 72.2, E = 2.3e-21
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpe.....eare
ERa IARELHDs++Q+Ls +k+q+ l++r + + p+++ +a +
BA000032_1 369 ERATIARELHDSLAQALSYLKIQVALLKRSV--KNLPDekaiaQANQ 413
aldeirelarealaevRrllgdLRpaal<-*
++e+ + ++ a++++R+ll+++R ++
BA000032_1 414 VIAELDTGLSAAYTQLRELLTTFRLTIK 441
BA000032_0965: domain 1 of 1, from 303 to 368: score 53.0, E = 1.1e-15
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
R+ IARELHD++gQ+++ai+ q+ ++r+ ++ + + +d+i
BA000032_0 303 VRRDIARELHDEIGQNITAIRTQANIIKRID--NAE--MSAHCADTI 345
relarealaevRrllgdLRpaal<-*
++l+ ++ + rll+ LRp l
BA000032_0 346 EGLSLNVYDTTKRLLSKLRPKML 368
Histogram of all scores:
score obs exp (one = represents 7 sequences)
----- --- ---
-15 1 0|=
-14 47 0|=======
-13 65 0|==========
-12 105 19|==*============
-11 143 160|===================== *
-10 327 372|=============================================== *
-9 360 423|==================================================== *
-8 231 326|================================= *
-7 169 205|========================= *
-6 136 116|================*===
-5 77 62|========*==
-4 32 32|====*
-3 26 16|==*=
-2 14 8|=*
-1 8 4|*=
0 4 2|*
> 1 7 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 439.5069
P(chi-square) = 0
Total sequences searched: 1752
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
