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: AP008957.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
-------- ----------- ----- ------- ---
AP008957_2189 CDS_2383636-2385138 128.1 4.3e-36 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
AP008957_2189 1/1 301 381 .. 1 94 [] 128.1 4.3e-36
Alignments of top-scoring domains:
AP008957_2189: domain 1 of 1, from 301 to 381: score 128.1, E = 4.3e-36
*->EIHHRVKNNLQvISSLLrLQarkfkdeevikdEVleaLrEsqnRVlS
EIHHRVKNNLQ++++LLrLQar++ + e+ + aL Es++RV S
AP008957_2 301 EIHHRVKNNLQTVAALLRLQARRTGNDEA-----RQALTESVRRVTS 342
mAliHEeLYksedldtvdfseYlekLtenLfqsygvgsdgisLklDl<-*
+Al+H+ L++s d + vd++e++++L+ L + +gvg +
AP008957_2 343 IALVHDTLSMSVD-EEVDLDEVIDRLVPMLSDVAGVGVP-------V 381
Histogram of all scores:
score obs exp (one = represents 29 sequences)
----- --- ---
-15 3 0|=
-14 19 0|=
-13 81 37|=*=
-12 341 396|============ *
-11 1658 1103|======================================*===================
-10 1393 1428|=================================================*
-9 976 1206|================================== *
-8 725 809|========================= *
-7 342 480|============ *
-6 285 267|=========*
-5 94 144|====*
-4 57 76|==*
-3 31 40|=*
-2 9 20|*
-1 12 10|*
0 3 5|*
1 2 2|*
2 0 1|*
3 1 0|=
> 4 2 -|=
% Statistical details of theoretical EVD fit:
mu = -9.5227
lambda = 0.6548
chi-sq statistic = 461.6715
P(chi-square) = 0
Total sequences searched: 6034
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