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_2.hmm [HisKA_2]
Sequence database: AL123456.Gprot.fsa
per-sequence score cutoff: >= 32.2 [TC1]
per-domain score cutoff: >= 32.2 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HisKA_2
Accession: PF07568.2
Description: Histidine kinase
[HMM has been calibrated; E-values are empirical estimates]
Scores for complete sequences (score includes all domains):
Sequence Description Score E-value N
-------- ----------- ----- ------- ---
AL123456_3283 CDS_3597534-3596029 118.2 1.9e-33 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
AL123456_3283 1/1 300 382 .. 1 94 [] 118.2 1.9e-33
Alignments of top-scoring domains:
AL123456_3283: domain 1 of 1, from 300 to 382: score 118.2, E = 1.9e-33
*->EIHHRVKNNLQvISSLLrLQarkfkdeevikdEVleaLrEsqnRVlS
EIHHRVKNNLQ++++LLrLQar++++ e +eaL Es++RV S
AL123456_3 300 EIHHRVKNNLQTVAALLRLQARRTSNAEG-----REALIESVRRVSS 341
mAliHEeLYksedldtvdfseYlekLtenLfqsygvgsdgisLklDl<-*
+Al+H+ L++s d ++v+++e++++++ + + ++v++ ++ +
AL123456_3 342 IALVHDALSMSVD-EQVNLDEVIDRILPIMNDVASVDRP---IR--I 382
Histogram of all scores:
score obs exp (one = represents 20 sequences)
----- --- ---
-15 6 0|=
-14 21 0|==
-13 77 24|=*==
-12 262 262|=============*
-11 1169 729|====================================*======================
-10 889 944|============================================= *
-9 561 797|============================= *
-8 501 535|==========================*
-7 222 317|============ *
-6 157 177|========*
-5 50 95|=== *
-4 34 50|==*
-3 21 26|=*
-2 9 13|*
-1 5 7|*
0 1 3|*
1 4 1|*
2 0 1|*
3 0 0|
> 4 2 -|=
% Statistical details of theoretical EVD fit:
mu = -9.5227
lambda = 0.6548
chi-sq statistic = 509.7957
P(chi-square) = 0
Total sequences searched: 3999
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
