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: FM162591.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
-------- ----------- ----- ------- ---
FM162591_0765 CDS_873439-871895 55.1 6.6e-16 1
FM162591_2548 CDS_2919442-2920671 48.2 7e-14 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
FM162591_0765 1/1 316 380 .. 2 70 .] 55.1 6.6e-16
FM162591_2548 1/1 213 281 .. 1 70 [] 48.2 7e-14
Alignments of top-scoring domains:
FM162591_0765: domain 1 of 1, from 316 to 380: score 55.1, E = 6.6e-16
*->RaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdeir
R+ IARELHD++gQ+++ai +q+ l++r + +++ + +++i
FM162591_0 316 RKSIARELHDEIGQNITAIQIQATLVKRTA--TSP--TTEMIAEQIN 358
elarealaevRrllgdLRpaal<-*
+la++ + +R+ll++LRp++l
FM162591_0 359 HLAQQIHQSTRQLLRQLRPPVL 380
FM162591_2548: domain 1 of 1, from 213 to 281: score 48.2, E = 7e-14
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ERa I+RELHDs++QsLs +k+q+ + d + + l +i
FM162591_2 213 ERASISRELHDSIAQSLSFLKIQVNCLYICS--DNLSKQSLASLNMI 257
relarealaevRrllgdLR.paal<-*
++ ++ a++++R+ll ++R+++
FM162591_2 258 KDELNIAYRQLRELLAAFRnEQSQ 281
Histogram of all scores:
score obs exp (one = represents 14 sequences)
----- --- ---
-15 2 0|=
-14 80 0|======
-13 177 0|=============
-12 261 49|===*===============
-11 354 402|========================== *
-10 797 933|========================================================= *
-9 815 1060|==========================================================*
-8 618 817|============================================= *
-7 472 513|================================== *
-6 289 290|====================*
-5 231 155|===========*=====
-4 126 81|=====*===
-3 72 41|==*===
-2 34 21|=*=
-1 30 10|*==
0 15 5|*=
1 13 2|*
2 0 1|*
3 1 0|=
4 0 0|
> 5 3 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 1190.0350
P(chi-square) = 0
Total sequences searched: 4395
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