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: AP009178.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
-------- ----------- ----- ------- ---
AP009178_0514 CDS_503908-502079 114.7 8.8e-33 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
AP009178_0514 1/1 425 509 .. 1 94 [] 114.7 8.8e-33
Alignments of top-scoring domains:
AP009178_0514: domain 1 of 1, from 425 to 509: score 114.7, E = 8.8e-33
*->EIHHRVKNNLQvISSLLrLQarkfkdeevikdEVleaLrEsqnRVlS
E HHRVKNNLQ+I+SLLrLQ++k ++++ + +L E+++R+++
AP009178_0 425 ELHHRVKNNLQIIISLLRLQSDKLSNQKM-----QSILTEAESRIRA 466
mAliHEeLYksedldtvdfseYlekLtenLfqsygvgsdgisLklDl<-*
+ +HE+LYks+++ +d eY++ L+++L ++++ +++ i+ +
AP009178_0 467 LSSVHEMLYKSDNIAYIDAQEYIQTLASELQHLFSQKDIAIH----V 509
Histogram of all scores:
score obs exp (one = represents 7 sequences)
----- --- ---
-14 3 0|=
-13 13 11|=*
-12 61 121|========= *
-11 387 337|================================================*=======
-10 373 436|====================================================== *
-9 349 368|================================================== *
-8 245 247|===================================*
-7 158 146|====================*==
-6 119 81|===========*=====
-5 54 44|======*=
-4 34 23|===*=
-3 23 12|=*==
-2 10 6|*=
-1 5 3|*
0 2 1|*
1 3 0|=
2 1 0|=
3 1 0|=
4 0 0|
> 5 2 -|=
% Statistical details of theoretical EVD fit:
mu = -9.5227
lambda = 0.6548
chi-sq statistic = 84.0461
P(chi-square) = 2.423e-13
Total sequences searched: 1857
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
