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: CP000437.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
-------- ----------- ----- ------- ---
CP000437_1405 CDS_1695374-1696801 57.0 3.1e-16 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000437_1405 1/1 249 314 .. 1 87 [] 57.0 3.1e-16
Alignments of top-scoring domains:
CP000437_1405: domain 1 of 1, from 249 to 314: score 57.0, E = 3.1e-16
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ +f +++HEL+TP++++++++e++l+ + +++
CP000437_1 249 REKRFSGDAAHELKTPIAGVKTLVEIALA----------SDDIDEI- 284
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++ le+i+++++r +++i++lL+lsri+ +
CP000437_1 285 ----------KQKLERIKSSTNRYSHIIDQLLTLSRIQPN 314
Histogram of all scores:
score obs exp (one = represents 5 sequences)
----- --- ---
-15 1 0|=
-14 10 0|==
-13 43 1|*========
-12 88 35|======*===========
-11 192 184|====================================*==
-10 282 341|========================================================= *
-9 286 353|==========================================================*
-8 224 264|============================================= *
-7 159 166|================================ *
-6 140 95|==================*=========
-5 56 52|==========*=
-4 27 28|=====*
-3 18 14|==*=
-2 4 7|=*
-1 7 4|*=
0 7 2|*=
1 6 1|*=
2 3 0|=
> 3 2 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 128.7594
P(chi-square) = 2.134e-23
Total sequences searched: 1555
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