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:          CR628338.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 

--------      -----------                               -----    ------- ---

CR628338_0099 CDS_97498-95894                            45.9    4.1e-14   1



Parsed for domains:

Sequence      Domain  seq-f seq-t    hmm-f hmm-t      score  E-value

--------      ------- ----- -----    ----- -----      -----  -------

CR628338_0099   1/1     315   381 ..     1    87 []    45.9  4.1e-14



Alignments of top-scoring domains:

CR628338_0099: domain 1 of 1, from 315 to 381: score 45.9, E = 4.1e-14

                   *->akseflanlsHELRTPLtaIrgyaelLleeldkvlagdtelseeqrl

                      +  ++ + ++HE+RTPL++I+++++ +++          ++ +  + 

  CR628338_0   315    SHKQLCRFVAHEIRTPLSTIQMATDSIKR---------RNAEDALL- 351  



                   peLRledpeyreyletilreaerllrLindlLdlsrieag<-*

                              + +++i++++  ++rL++ +L +s++ ++   

  CR628338_0   352 ----------NKQINSIQEDIADMNRLVSTFLIYSKMHSN    381  





Histogram of all scores:

score    obs    exp  (one = represents 1 sequences)

-----    ---    ---

  -14      2      0|==                                                         

  -13      3      0|===                                                        

  -12     10      3|==*=======                                                 

  -11     18     16|===============*==                                         

  -10     23     30|=======================      *                             

   -9     26     31|==========================    *                            

   -8     20     23|====================  *                                    

   -7     17     14|=============*===                                          

   -6      8      8|=======*                                                   

   -5      7      4|===*===                                                    

   -4      1      2|=*                                                         

>  -3      4      -|====                                                       





% Statistical details of theoretical EVD fit:

              mu =    -8.9758

          lambda =     0.6525

chi-sq statistic =     3.8899

  P(chi-square)  =     0.5654



Total sequences searched: 142



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