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_2.hmm [HisKA_2]
Sequence database: CP000518.Gprot.fsa
per-sequence score cutoff: >= 32.2 [TC1]
per-domain score cutoff: >= 32.2 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HisKA_2
Accession: PF07568.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
-------- ----------- ----- ------- ---
CP000518_1393 CDS_1501440-1502930 117.6 3.8e-33 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000518_1393 1/1 295 379 .. 1 94 [] 117.6 3.8e-33
Alignments of top-scoring domains:
CP000518_1393: domain 1 of 1, from 295 to 379: score 117.6, E = 3.8e-33
*->EIHHRVKNNLQvISSLLrLQarkfkdeevikdEVleaLrEsqnRVlS
EIHHRVKNNLQ++++LLrLQar++ + e +eaL+Es++RV S
CP000518_1 295 EIHHRVKNNLQTVAALLRLQARRTNNAEG-----REALMESVRRVSS 336
mAliHEeLYksedldtvdfseYlekLtenLfqsygvgsdgisLklDl<-*
+A++H+ L++s d + v+++e++++++ + + ++vg i++k
CP000518_1 337 IAQVHDALSMSVD-EEVNLDEVVDRILPIMNDVASVGPP-IRIK--R 379
Histogram of all scores:
score obs exp (one = represents 28 sequences)
----- --- ---
-15 3 0|=
-14 20 0|=
-13 83 34|=*=
-12 325 358|============*
-11 1606 998|===================================*======================
-10 1252 1292|============================================= *
-9 824 1091|============================== *
-8 657 732|======================== *
-7 278 434|========== *
-6 223 242|========*
-5 74 130|=== *
-4 56 69|==*
-3 28 36|=*
-2 12 18|*
-1 14 9|*
0 2 5|*
1 1 2|*
> 2 2 -|=
% Statistical details of theoretical EVD fit:
mu = -9.5227
lambda = 0.6548
chi-sq statistic = 608.0228
P(chi-square) = 0
Total sequences searched: 5460
Whole sequence top hits:
tophits_s report:
Total hits: 1
Satisfying E cutoff: 1
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 1
Satisfying E cutoff: 1
Total memory: 21K