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: AL450380.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
-------- ----------- ----- ------- ---
AL450380_0775 CDS_916683-918371 75.1 4e-21 1
AL450380_2446 CDS_2916453-2915128 66.9 8.8e-19 1
AL450380_0175 CDS_241217-242776 57.8 3.2e-16 1
AL450380_2128 CDS_2524869-2526209 49.7 6.5e-14 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
AL450380_0775 1/1 295 362 .. 1 87 [] 75.1 4e-21
AL450380_2446 1/1 155 221 .. 1 87 [] 66.9 8.8e-19
AL450380_0175 1/1 239 311 .. 1 87 [] 57.8 3.2e-16
AL450380_2128 1/1 230 296 .. 1 87 [] 49.7 6.5e-14
Alignments of top-scoring domains:
AL450380_0775: domain 1 of 1, from 295 to 362: score 75.1, E = 4e-21
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+++f++++sHELRTPLt++r++a+l+ + ++++l++
AL450380_0 295 LQRRFTSDVSHELRTPLTTVRMAADLIYD-------HSSDLDPTL-- 332
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
r+ +e + +e++r + L+ndlL++sr +ag
AL450380_0 333 ----------RRSTELMVSELDRFETLLNDLLEISRHDAG 362
AL450380_2446: domain 1 of 1, from 155 to 221: score 66.9, E = 8.8e-19
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a+++f+an+sHEL+TP+ a+ ++ae+Ll+ ++ ++e++
AL450380_2 155 ARRDFVANVSHELKTPVGAMALLAEALLA---------SADDAETV- 191
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
++ e++l ea+rl ++ +l +lsr++
AL450380_2 192 ----------SRFAEKVLIEANRLGYMVAELIELSRLQGA 221
AL450380_0175: domain 1 of 1, from 239 to 311: score 57.8, E = 3.2e-16
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
++ +++ ++ HELRTPLt++r+++elL + + +++ + l+++++
AL450380_0 239 RQARLVTDAGHELRTPLTSLRTNVELLIASM---APEAPRLPDQEM- 281
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+++ +l ++e l+ L+ dl dl+r +ag
AL450380_0 282 ----------ADLRADVLAQIEELSTLVGDLVDLTRDDAG 311
AL450380_2128: domain 1 of 1, from 230 to 296: score 49.7, E = 6.5e-14
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
+ ++f+a sHELRTPLta+r+++e+L + +l ++qr
AL450380_2 230 SARDFAAVSSHELRTPLTAMRTNLEVLAT---------LDLADDQR- 266
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
+e l ++r r+++ ++ l l++ e
AL450380_2 267 ----------KEVLGDVIRTQSRIEATLSALERLAQGELS 296
Histogram of all scores:
score obs exp (one = represents 6 sequences)
----- --- ---
-14 3 0|=
-13 24 1|*===
-12 74 36|=====*=======
-11 174 190|============================= *
-10 278 352|=============================================== *
-9 315 364|===================================================== *
-8 271 272|=============================================*
-7 192 171|============================*===
-6 114 98|================*==
-5 70 54|========*===
-4 43 29|====*===
-3 20 15|==*=
-2 14 8|=*=
-1 1 4|*
0 3 2|*
1 1 1|*
2 1 0|=
3 1 0|=
4 1 0|=
5 0 0|
> 6 5 -|=
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 83.6424
P(chi-square) = 9.664e-14
Total sequences searched: 2720
Whole sequence top hits:
tophits_s report:
Total hits: 4
Satisfying E cutoff: 4
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 4
Satisfying E cutoff: 4
Total memory: 22K
