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_1.hmm [HisKA]
Sequence database: CP001505.Gprot.fsa
per-sequence score cutoff: >= 26.8 [TC1]
per-domain score cutoff: >= 26.8 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HisKA
Accession: PF00512.12
Description: His Kinase A (phosphoacceptor) domain
[HMM has been calibrated; E-values are empirical estimates]
Scores for complete sequences (score includes all domains):
Sequence Description Score E-value N
-------- ----------- ----- ------- ---
CP001505_0056 CDS_76271-73584 72.1 1.1e-21 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001505_0056 1/1 447 510 .. 1 87 [] 72.1 1.1e-21
Alignments of top-scoring domains:
CP001505_0056: domain 1 of 1, from 447 to 510: score 72.1, E = 1.1e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
aks fla +sH+LRTPLt+I+g+a+l+++ d+ ++e
CP001505_0 447 AKSNFLARISHDLRTPLTSIIGFADLIQA--------DGREDAER-- 483
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++i+r+a+++l ++ndl d++r ++
CP001505_0 484 -------------GRVIRRSANHMLGMVNDLIDYARGREA 510
Histogram of all scores:
score obs exp (one = represents 1 sequences)
----- --- ---
-12 10 2|=*========
-11 10 12|========== *
-10 15 22|=============== *
-9 14 23|============== *
-8 15 17|=============== *
-7 19 10|=========*=========
-6 12 6|=====*======
-5 2 3|==*
-4 1 1|*
-3 1 0|=
> -2 3 -|===
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 17.9065
P(chi-square) = 0.003066
Total sequences searched: 102
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