# 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/Methyltransf_15.hmm.txt
# target sequence database:        proteomes/Drosophila_melanogaster.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query:       Methyltransf_15  [M=165]
Accession:   PF09445.9
Description: RNA cap guanine-N2 methyltransferase
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
    ------- ------ -----    ------- ------ -----   ---- --  --------               -----------
    1.7e-44  151.1   0.0    3.6e-44  150.0   0.0    1.6  1  tr|Q9VEA4|Q9VEA4_DROME  Trimethylguanosine synthase 1 OS=Dros
    2.3e-05   23.8   0.0    4.5e-05   22.8   0.0    1.5  1  tr|Q9VVV7|Q9VVV7_DROME  Uncharacterized protein OS=Drosophila
  ------ inclusion threshold ------
      0.023   14.0   0.0      0.048   13.0   0.0    1.5  1  sp|Q8IRE4|TRM5_DROME    tRNA (guanine(37)-N1)-methyltransfera


Domain annotation for each sequence (and alignments):
>> tr|Q9VEA4|Q9VEA4_DROME  Trimethylguanosine synthase 1 OS=Drosophila melanogaster OX=7227 GN=Tgs1 PE=4 SV=3
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  150.0   0.0   7.9e-48   3.6e-44       2     164 ..     300     455 ..     299     456 .. 0.95

  Alignments for each domain:
  == domain 1  score: 150.0 bits;  conditional E-value: 7.9e-48
         Methyltransf_15   2 evvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpPwGGpsYkk 88 
                              v++D+f+G+GGn+iqfan+  +via+di+ ++l++a+hNa +Ygv++k++++++d+ +  a++kl+     +vflSpPwGGp+Y+k
  tr|Q9VEA4|Q9VEA4_DROME 300 DVIVDAFCGCGGNAIQFANTCGRVIAVDIDAEKLAMAKHNAGIYGVAHKIEFIHADFLQFAASTKLRPN---VVFLSPPWGGPDYQK 383
                             689**********************************************************99998877...*************** PP

         Methyltransf_15  89 kekfdlekkleplklekllkesvklskniilfLPrnsdlkqlaaltaevlenaskvkvvklkenglvkallaylge 164
                             + +fd+e+ l p+  ++l++ s +l + + +fLPrn+++kq+ al     +  ++++v +   + ++ al+ay+g+
  tr|Q9VEA4|Q9VEA4_DROME 384 QATFDIETGLLPVGASQLMQLSRSLASDVAFFLPRNANMKQVVAL----SGVGQQCEVEHNYLDTRMVALTAYYGK 455
                             ******************************************998....577788999999999999999999986 PP

>> tr|Q9VVV7|Q9VVV7_DROME  Uncharacterized protein OS=Drosophila melanogaster OX=7227 GN=cg3808 PE=1 SV=3
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !   22.8   0.0   9.8e-09   4.5e-05       3      65 ..     418     479 ..     416     526 .. 0.84

  Alignments for each domain:
  == domain 1  score: 22.8 bits;  conditional E-value: 9.8e-09
         Methyltransf_15   3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellakl 65 
                             ++lD+++G+G  + +fa+  +kv++++i ++ +k a+ Nae+ g+++  ++ +g+  + ++++
  tr|Q9VVV7|Q9VVV7_DROME 418 TMLDICCGTGTITLAFAKHCKKVMGVEIVPDAIKDAEFNAEANGIKN-AKFFTGNADDFIKSM 479
                             79*****************************************9964.456666666655544 PP

>> sp|Q8IRE4|TRM5_DROME  tRNA (guanine(37)-N1)-methyltransferase OS=Drosophila melanogaster OX=7227 GN=CG32281 PE=2 SV=2
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?   13.0   0.0     1e-05     0.048       3      60 ..     252     309 ..     250     341 .. 0.86

  Alignments for each domain:
  == domain 1  score: 13.0 bits;  conditional E-value: 1e-05
       Methyltransf_15   3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfe 60 
                           v+ D+f+G G  +i  a+  ++v a d n++  +  qhNa+      ++++   d  +
  sp|Q8IRE4|TRM5_DROME 252 VLYDVFAGVGPFSIPAAKKRCHVLANDLNPESFRWLQHNAKRNKCLPNIKMSNKDGRQ 309
                           567**************************************99999999887777665 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (165 nodes)
Target sequences:                        13786  (7394964 residues searched)
Passed MSV filter:                       375  (0.0272015); expected 275.7 (0.02)
Passed bias filter:                      300  (0.0217612); expected 275.7 (0.02)
Passed Vit filter:                        32  (0.0023212); 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):               3  [number of targets reported over threshold]
# CPU time: 0.10u 0.00s 00:00:00.10 Elapsed: 00:00:00.04
# Mc/sec: 30504.23
//
[ok]