# hmmsearch :: search profile(s) against a sequence database
# HMMER 3.1b2 (February 2015); http://hmmer.org/
# Copyright (C) 2015 Howard Hughes Medical Institute.
# Freely distributed under the GNU General Public License (GPLv3).
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query HMM file:                  pfam_list/TRM13.hmm.txt
# target sequence database:        proteomes/Drosophila_melanogaster.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query:       TRM13  [M=262]
Accession:   PF05206.13
Description: Methyltransferase TRM13
Scores for complete sequences (score includes all domains):
   --- full sequence ---   --- best 1 domain ---    -#dom-
    E-value  score  bias    E-value  score  bias    exp  N  Sequence               Description
    ------- ------ -----    ------- ------ -----   ---- --  --------               -----------
    4.4e-91  304.7   1.9    8.6e-91  303.7   1.9    1.5  1  tr|Q8SX24|Q8SX24_DROME  RE67941p OS=Drosophila melanogaster O
     0.0029   16.8   3.3        1.6    7.8   0.0    3.0  3  sp|Q9VJ34|GSTCD_DROME   Glutathione S-transferase C-terminal 


Domain annotation for each sequence (and alignments):
>> tr|Q8SX24|Q8SX24_DROME  RE67941p OS=Drosophila melanogaster OX=7227 GN=Dmel\CG18048 PE=2 SV=1
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  303.7   1.9   1.2e-94   8.6e-91       2     261 ..     184     428 ..     183     429 .. 0.97

  Alignments for each domain:
  == domain 1  score: 303.7 bits;  conditional E-value: 1.2e-94
                   TRM13   2 qqsslignleelkllekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnls 88 
                             q++sl+g+le+ ++l +++ y+EfGaGkg+l+ +++++l+e+k ++s++vLidr+s r+k+D+kl+++e   +++R+++dI dL+ls
  tr|Q8SX24|Q8SX24_DROME 184 QTTSLLGILEHDHQLMDHTSYIEFGAGKGQLAYFLATVLQEQKLSHSQVVLIDRMSLRHKKDNKLANRE---VVQRIRADIADLKLS 267
                             899******************************************************************...8************** PP

                   TRM13  89 aleelekkkkvvavsKHLCGaATDLtLrcllnsekaskkakleglliAlCChhvCswkeyvnkeyleelgitkeefqiltkmvSWav 175
                             al+el+k++++va sKHLCGaATDLtLrc+l   +as     +++liAlCChh+Csw++yv++++l+e gi  +ef iltkmvSWav
  tr|Q8SX24|Q8SX24_DROME 268 ALPELKKTQRTVAFSKHLCGAATDLTLRCILGDGNAS----SDYVLIALCCHHRCSWRSYVGRKFLQEAGIGPREFVILTKMVSWAV 350
                             *********************************9986....59******************************************** PP

                   TRM13 176 sgkrkseaeeedekeekekeekkeeeeseeekksklsseereeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLla 261
                             +g+          ++e++k++++++ + +e+++++l+++ere+iG+++kr++d+gRl++l+++g++ael++Yv+++v+lEnv+Lla
  tr|Q8SX24|Q8SX24_DROME 351 CGTG--------LSRERRKAMESADFQLTETNTQRLTRQEREQIGQQCKRVLDYGRLEHLRSHGYQAELKFYVPRDVTLENVVLLA 428
                             ****........56788899999999**********************************************************98 PP

>> sp|Q9VJ34|GSTCD_DROME  Glutathione S-transferase C-terminal domain-containing protein homolog OS=Drosophila melanogas
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?   -0.4   0.0     0.076   5.2e+02      16      42 ..     384     410 ..     368     433 .. 0.70
   2 !    5.7   2.3    0.0011       7.3     100     122 ..     454     476 ..     449     489 .. 0.84
   3 !    7.8   0.0   0.00023       1.6     218     260 ..     536     578 ..     511     580 .. 0.89

  Alignments for each domain:
  == domain 1  score: -0.4 bits;  conditional E-value: 0.076
                  TRM13  16 lekssayvEfGaGkgelsryvnqalee 42 
                             ++   +v+f++G g+l+  ++  l +
  sp|Q9VJ34|GSTCD_DROME 384 AQPGDRIVDFCSGTGHLAILLALKLPN 410
                            3677899**********9988876644 PP

  == domain 2  score: 5.7 bits;  conditional E-value: 0.0011
                  TRM13 100 vavsKHLCGaATDLtLrcllnse 122
                            ++ s H CG+ATD+ L+ + + +
  sp|Q9VJ34|GSTCD_DROME 454 IGASLHACGTATDIVLQQCRRAK 476
                            57799***********9887665 PP

  == domain 3  score: 7.8 bits;  conditional E-value: 0.00023
                  TRM13 218 eiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLl 260
                              Gl +  ++d+ Rll  +e g+++ l+   +++ ++ n+lL+
  sp|Q9VJ34|GSTCD_DROME 536 LQGLHCMSVVDTDRLLQAEEAGYQVILTRLKPEQCTPKNHLLV 578
                            4699************************************997 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (262 nodes)
Target sequences:                        13786  (7394964 residues searched)
Passed MSV filter:                      1357  (0.0984332); expected 275.7 (0.02)
Passed bias filter:                      511  (0.0370666); expected 275.7 (0.02)
Passed Vit filter:                        50  (0.00362687); expected 13.8 (0.001)
Passed Fwd filter:                         3  (0.000217612); expected 0.1 (1e-05)
Initial search space (Z):              13786  [actual number of targets]
Domain search space  (domZ):               2  [number of targets reported over threshold]
# CPU time: 0.21u 0.01s 00:00:00.22 Elapsed: 00:00:00.07
# Mc/sec: 27678.29
//
[ok]