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: CP001928.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
-------- ----------- ----- ------- ---
CP001928_1116 CDS_1213356-1216373 94.7 7.9e-27 1
CP001928_1831 CDS_2011429-2013189 60.6 3.8e-17 1
CP001928_1542 CDS_1648685-1649323 30.1 1.7e-08 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001928_1116 1/1 780 845 .. 1 87 [] 94.7 7.9e-27
CP001928_1831 1/1 361 428 .. 1 87 [] 60.6 3.8e-17
CP001928_1542 1/1 10 76 .. 1 87 [] 30.1 1.7e-08
Alignments of top-scoring domains:
CP001928_1116: domain 1 of 1, from 780 to 845: score 94.7, E = 7.9e-27
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
aks flan+sHE+RTPL++++g+++lLl+ t+l+ q
CP001928_1 780 AKSTFLANMSHEIRTPLNGVIGMTSLLLN---------TPLDKSQ-- 815
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++y ++i+ +++ ll++ind+Ld+s+ieag
CP001928_1 816 ----------EKYANRIDLSGKVLLEIINDILDFSKIEAG 845
CP001928_1831: domain 1 of 1, from 361 to 428: score 60.6, E = 3.8e-17
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ +f+an+sHEL+TP+t I+g+ae+L ++ + l+ e+
CP001928_1 361 MRKDFVANASHELKTPITIIQGFAEALHDN--------PGLPIETS- 398
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
e++++i r + Li dlL+l+ ie+
CP001928_1 399 ----------SEITSKIVRNCIKMTMLIKDLLTLTDIEHI 428
CP001928_1542: domain 1 of 1, from 10 to 76: score 30.1, E = 1.7e-08
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a +e ++ l H LR+PL +I+g+a+lL++ +l + +l
CP001928_1 10 AFQEAASLLIHRLRNPLGGIKGFASLLKR----------DLEGQPNL 46
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++ i + a+ l+r+ Ld++ +
CP001928_1 47 ----------QKMAAMIVEGADNLDRIMTAALDFAEPLNL 76
Histogram of all scores:
score obs exp (one = represents 6 sequences)
----- --- ---
-14 8 0|==
-13 31 1|*=====
-12 98 44|=======*=========
-11 202 229|================================== *
-10 337 425|========================================================= *
-9 339 439|========================================================= *
-8 331 328|======================================================*=
-7 222 207|==================================*==
-6 134 119|===================*===
-5 93 65|==========*=====
-4 54 35|=====*===
-3 38 18|==*====
-2 20 9|=*==
-1 8 5|*=
0 6 2|*
1 3 1|*
2 6 0|=
> 3 4 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 167.9002
P(chi-square) = 3.19e-30
Total sequences searched: 1934
Whole sequence top hits:
tophits_s report:
Total hits: 3
Satisfying E cutoff: 3
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 3
Satisfying E cutoff: 3
Total memory: 22K
