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: CP001138.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
-------- ----------- ----- ------- ---
CP001138_1285 CDS_1338587-1336746 74.3 1.5e-21 1
CP001138_2453 CDS_2552406-2554106 72.9 3.7e-21 1
CP001138_3746 CDS_3915066-3913564 65.1 7.8e-19 1
Parsed for domains:
Sequence Domain seq-f seq-t hmm-f hmm-t score E-value
-------- ------- ----- ----- ----- ----- ----- -------
CP001138_1285 1/1 405 473 .. 1 70 [] 74.3 1.5e-21
CP001138_2453 1/1 361 428 .. 1 70 [] 72.9 3.7e-21
CP001138_3746 1/1 304 370 .. 1 70 [] 65.1 7.8e-19
Alignments of top-scoring domains:
CP001138_1285: domain 1 of 1, from 405 to 473: score 74.3, E = 1.5e-21
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ERa IARELHDs++QsLs +k+q++ ++++ d pe re l +i
CP001138_1 405 ERATIARELHDSIAQSLSCMKMQVSCLQMQG--DVLPESSRELLSQI 449
relarealaevRrllgdLR.paal<-*
r ++ + a++R+ll+++R ++
CP001138_1 450 RNELNCSWAQLRELLTTFRlQLTE 473
CP001138_2453: domain 1 of 1, from 361 to 428: score 72.9, E = 3.7e-21
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
ERa IARELHDs++Q Ls +++ql l++r++ ++d++ a+ ++
CP001138_2 361 ERATIARELHDSLAQVLSYLRIQLTLLKRAI--PEDNAGAQSIMADF 405
relarealaevRrllgdLRpaal<-*
+++++a++++R+ll+++R +
CP001138_2 406 SRALNDAYRQLRELLTTFRLTLQ 428
CP001138_3746: domain 1 of 1, from 304 to 370: score 65.1, E = 7.8e-19
*->ERaRIARELHDsvgQsLsaiklqlelarrlldsdkdpeearealdei
R+ ARELHD +gQ+++ai+ q++ ++rl+ +d++ +++ + i
CP001138_3 304 VRRDVARELHDDIGQTITAIRTQAGIVQRLA---ADNGGVKQSGQLI 347
relarealaevRrllgdLRpaal<-*
++l+ ++ vRrllg+LRp +l
CP001138_3 348 EQLSLGVYDAVRRLLGRLRPRQL 370
Histogram of all scores:
score obs exp (one = represents 14 sequences)
----- --- ---
-15 9 0|=
-14 97 0|=======
-13 164 0|============
-12 216 50|===*============
-11 298 418|====================== *
-10 754 970|====================================================== *
-9 808 1101|==========================================================*
-8 633 849|============================================== *
-7 530 533|======================================*
-6 396 302|=====================*=======
-5 314 162|===========*===========
-4 139 84|=====*====
-3 75 43|===*==
-2 55 22|=*==
-1 29 11|*==
0 21 5|*=
1 8 2|*
2 4 1|*
3 4 0|=
4 1 0|=
5 3 0|=
> 6 4 -|=
% Statistical details of theoretical EVD fit:
mu = -8.7900
lambda = 0.6787
chi-sq statistic = 1096.3506
P(chi-square) = 0
Total sequences searched: 4562
Whole sequence top hits:
tophits_s report:
Total hits: 3
Satisfying E cutoff: 3
Total memory: 20K
Domain top hits:
tophits_s report:
Total hits: 3
Satisfying E cutoff: 3
Total memory: 21K
