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: CP000025.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
-------- ----------- ----- ------- ---
CP000025_0937 CDS_907678-906389 47.2 2.2e-13 1
CP000025_0857 CDS_821817-822836 38.6 6e-11 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000025_0937 1/1 224 291 .. 1 87 [] 47.2 2.2e-13
CP000025_0857 1/1 131 195 .. 1 87 [] 38.6 6e-11
Alignments of top-scoring domains:
CP000025_0937: domain 1 of 1, from 224 to 291: score 47.2, E = 2.2e-13
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ ++++sHEL++PL++I + +e+L++ ++l ++++
CP000025_0 224 QLASVISAISHELKNPLSVIDLSLEMLKD---------EKLEDKKLK 261
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+e+le+i r++ l++L l + ++
CP000025_0 262 ----------KELLEKISRQSVKLNALTHKLNFVFNLNSE 291
CP000025_0857: domain 1 of 1, from 131 to 195: score 38.6, E = 6e-11
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++ +a+l+HE+R+P+ +I+++a++L + +se ++
CP000025_0 131 SMGSVAAHLAHEIRNPIGSISLLASTLFA-----------RSELKN- 165
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ ++ i++++ r +r++n L +++ +
CP000025_0 166 ----------KHIVLEIQKAIARVERIVNSTLLFTKGVHI 195
Histogram of all scores:
score obs exp (one = represents 6 sequences)
----- --- ---
-15 5 0|=
-14 9 0|==
-13 41 1|*======
-12 101 42|======*==========
-11 212 217|====================================*
-10 291 404|================================================= *
-9 347 417|==========================================================*
-8 297 312|================================================== *
-7 210 196|================================*==
-6 130 113|==================*===
-5 101 62|==========*======
-4 36 33|=====*
-3 27 17|==*==
-2 11 9|=*
-1 2 4|*
0 4 2|*
1 4 1|*
2 1 0|=
3 3 0|=
4 0 0|
5 0 0|
> 6 6 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 160.5010
P(chi-square) = 2.554e-29
Total sequences searched: 1838
Whole sequence top hits:
tophits_s report:
Total hits: 2
Satisfying E cutoff: 2
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 2
Satisfying E cutoff: 2
Total memory: 21K
