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: AE007871.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
-------- ----------- ----- ------- ---
AE007871_0001 CDS_1-2502 50.1 3.5e-15 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
AE007871_0001 1/1 468 534 .. 1 87 [] 50.1 3.5e-15
Alignments of top-scoring domains:
AE007871_0001: domain 1 of 1, from 468 to 534: score 50.1, E = 3.5e-15
*->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl
a + ++ +++HE+++ L aI+gyae++++ l +s+ +
AE007871_0 468 AVGTLAGGIAHEFNNILGAILGYAEMAQNML--------RRSSVT-- 504
peLRledpeyreyletilreaerllrLindlLdlsrieag<-*
r+ ++ i+++++r + +i+++L+lsr +
AE007871_0 505 ----------RRHIDQIISSGDRARLIIDQILTLSRKLER 534
Histogram of all scores:
score obs exp (one = represents 1 sequences)
----- --- ---
-13 1 0|=
-12 6 4|===*==
-11 23 23|======================*
-10 30 43|============================== *
-9 34 44|================================== *
-8 32 33|================================*
-7 27 21|====================*======
-6 20 12|===========*========
-5 9 6|=====*===
-4 5 3|==*==
-3 3 1|*==
-2 4 0|====
-1 1 0|=
0 0 0|
1 1 0|=
2 0 0|
> 3 2 -|==
% Statistical details of theoretical EVD fit:
mu = -8.9758
lambda = 0.6525
chi-sq statistic = 14.3142
P(chi-square) = 0.02632
Total sequences searched: 198
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
