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_3.hmm [HisKA_3]
Sequence database: CP000103.Gprot.fsa
per-sequence score cutoff: >= 30.1 [TC1]
per-domain score cutoff: >= 30.1 [TC2]
per-sequence Eval cutoff: [none]
per-domain Eval cutoff: [none]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query HMM: HisKA_3
Accession: PF07730.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
-------- ----------- ----- ------- ---
CP000103_1328 CDS_1524099-1525043 51.5 4.6e-15 1
CP000103_1230 CDS_1422090-1420012 42.5 2.1e-12 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP000103_1328 1/1 164 236 .. 1 70 [] 51.5 4.6e-15
CP000103_1230 1/1 454 521 .. 1 70 [] 42.5 2.1e-12
Alignments of top-scoring domains:
CP000103_1328: domain 1 of 1, from 164 to 236: score 51.5, E = 4.6e-15
*->ERaRIARELHDsvgQsLsaiklqlelarrllds.dkdpe..eareal
ER RIARE+HD +g +L+aik+ l+++ r++d +++ +e+ + + +
CP000103_1 164 ERIRIAREIHDDLGGNLTAIKIDLDWLVRRIDAgASSAEnaVLLAKV 210
deirelarealaevRrllgdLRpaal<-*
+ + +l+++++ + r+ +dLRp++
CP000103_1 211 RIVSDLVDRTIHSIQRISRDLRPGIM 236
CP000103_1230: domain 1 of 1, from 454 to 521: score 42.5, E = 2.1e-12
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
R+ IARE+HD++gQ L+a+ l ++l+r + d e++++ +++
CP000103_1 454 DRKYIAREIHDELGQLLAALHLEISLLRSAK--DHGSERLETIRHNM 498
relarealaevRrllgdLRpaal<-*
+el+++a + v +++ +LRpa l
CP000103_1 499 AELVEKAGQSVHNVAEHLRPASL 521
Histogram of all scores:
score obs exp (one = represents 9 sequences)
----- --- ---
-14 4 0|=
-13 49 0|======
-12 86 30|===*======
-11 149 252|================= *
-10 398 586|============================================= *
-9 509 665|========================================================= *
-8 405 513|============================================= *
-7 364 322|===================================*=====
-6 268 182|====================*=========
-5 224 97|==========*==============
-4 129 51|=====*=========
-3 80 26|==*======
-2 44 13|=*===
-1 20 6|*==
0 11 3|*=
1 4 1|*
2 5 0|=
> 3 8 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 792.5217
P(chi-square) = 0
Total sequences searched: 2757
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
