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: FM991728.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
-------- ----------- ----- ------- ---
FM991728_0235 CDS_244662-243517 43.5 2.1e-12 1
FM991728_1297 CDS_1337004-1335886 31.8 4.4e-09 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
FM991728_0235 1/1 170 234 .. 1 87 [] 43.5 2.1e-12
FM991728_1297 1/1 165 223 .. 1 87 [] 31.8 4.4e-09
Alignments of top-scoring domains:
FM991728_0235: domain 1 of 1, from 170 to 234: score 43.5, E = 2.1e-12
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++++ a+l+HE+R+P+ +I+++a++Ll+ +++e++
FM991728_0 170 SMGKISAHLAHEIRNPVGSISLLASVLLK----------HANEKT-- 204
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ ++ +++++ r +r+i L +s+ +
FM991728_0 205 ----------KPIVVELQKALWRVERIIKATLLFSKGIQA 234
FM991728_1297: domain 1 of 1, from 165 to 223: score 31.8, E = 4.4e-09
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ + fl n HEL+TP++a+ + +++L++
FM991728_1 165 RIDHFLKNTTHELNTPMSALVLSLKTLED------------------ 193
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+ ++i+ +++r++ L + l l + +
FM991728_1 194 ----------SQQHRRIKIAIQRMSFLYRSLSYLVMQDIE 223
Histogram of all scores:
score obs exp (one = represents 5 sequences)
----- --- ---
-14 5 0|=
-13 18 1|*===
-12 72 31|======*========
-11 136 164|============================ *
-10 239 304|================================================ *
-9 285 314|========================================================= *
-8 236 235|==============================================*=
-7 156 148|=============================*==
-6 97 85|================*===
-5 55 46|=========*=
-4 36 25|====*===
-3 18 13|==*=
-2 10 6|=*
-1 9 3|*=
0 7 1|*=
1 0 0|
2 1 0|=
3 0 0|
> 4 4 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 84.3414
P(chi-square) = 7.038e-14
Total sequences searched: 1528
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
