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

Query:       Methyltransf_10  [M=299]
Accession:   PF05971.11
Description: Protein of unknown function (DUF890)
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
    ------- ------ -----    ------- ------ -----   ---- --  --------               -----------
    2.4e-76  256.6   0.0    2.8e-76  256.3   0.0    1.0  1  sp|Q7K3B9|MET16_DROME   U6 small nuclear RNA (adenine-(43)-N(
  ------ inclusion threshold ------
       0.03   13.3   0.0      0.041   12.9   0.0    1.2  1  tr|Q9VQJ8|Q9VQJ8_DROME  Protein-lysine N-methyltransferase CG


Domain annotation for each sequence (and alignments):
>> sp|Q7K3B9|MET16_DROME  U6 small nuclear RNA (adenine-(43)-N(6))-methyltransferase OS=Drosophila melanogaster OX=7227 
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  256.3   0.0   4.1e-80   2.8e-76       5     298 ..      11     296 ..       7     297 .. 0.94

  Alignments for each domain:
  == domain 1  score: 256.3 bits;  conditional E-value: 4.1e-80
                            XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS----HHHHHHHHHHHHHH CS
        Methyltransf_10   5 sglhprnrhkdkydfaalikvvPelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPpvPgradyihyvadl 92 
                             g+hprn  + ++d+++++ ++ +++q  +l+++G+vs++f++ ++++ l+k ll+e+y++   d   g l+P++ +r++yi +++dl
  sp|Q7K3B9|MET16_DROME  11 VGMHPRNVLRTQPDYTKMAIKYKDFRQQCQLELNGKVSVNFRNEKTLRELTKMLLKEYYDLDV-DFAPGSLVPTLALRLNYILWLEDL 97 
                            58************************************************************9.************************ PP

                            HT--TCGCS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE--ST-SSTTTSTT--S...- CS
        Methyltransf_10  93 lgkadsdisaqvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaivesnsklsdaiklrkqkekkaifdglieene...r 177
                            +   + +   ++r++diG G++Ciy llG++++GW+ ++ e  ++ +++ak +v+ n +++++i++  q+++++if++  e+++   +
  sp|Q7K3B9|MET16_DROME  98 MEPLNLQ---NIRGIDIGCGSSCIYSLLGAKKNGWHMLALESKPQNIEYAKENVKRN-HMESLIEVYAQPDNTNIFKSYFEQDQqqlQ 181
                            9887765...59*********************************************.9********************998876678 PP

                            EEEEEE-----SS--.----------...---.....---------TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEEEEESSGGGHH CS
        Methyltransf_10 178 ydvtlCnPPfhasla.eakagserklknlgkanktralpslnfgGqkaelyceGGevafikklieeskqlakqvlWfttlvskkenla 264
                            y+++lCnPPf++s + +  +g+ r++   +++      p+   +G ++el c GGev+f++++i+es + +++v+ ftt+++ k n++
  sp|Q7K3B9|MET16_DROME 182 YQFCLCNPPFFDSNLpNPLGGNTRNP--ERRP-----APNNARTGSQEELTCVGGEVQFVQRIIDESLENKERVRIFTTMLGVKANVP 262
                            ***********987626678899999..3444.....5777899******************************************** PP

                            HHHHHHHHTT-SEEEEEEEEETTEEEEEEEEES- CS
        Methyltransf_10 265 plkeelkklGaskvtvvelaqGqklsrfiaWsfl 298
                             + + lk+l + +v ++e+ qG++ ++++aWsf 
  sp|Q7K3B9|MET16_DROME 263 RILDYLKELQVANVSTTEFHQGHTTRWAVAWSFH 296
                            ********************************96 PP

>> tr|Q9VQJ8|Q9VQJ8_DROME  Protein-lysine N-methyltransferase CG9643 OS=Drosophila melanogaster OX=7227 GN=Dmel\CG9643 P
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?   12.9   0.0   5.9e-06     0.041      92     181 ..      50     139 ..      43     173 .. 0.74

  Alignments for each domain:
  == domain 1  score: 12.9 bits;  conditional E-value: 5.9e-06
                             HHT--TCGCS---EEEEES-TTTTHHHHHHHHHH--E..EEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE..--ST-SSTTTSTT- CS
         Methyltransf_10  92 llgkadsdisaqvraldiGvGanCiyPllGvteyGWr..fvgsevddvslasakaivesnsklsdaiklrk..qkekkaifdgliee 174
                             ll +++ d+ +  r+ld+G G      l+G  + G++   +g + ++++++ a++i e n kls + k+    q++++     +++ 
  tr|Q9VQJ8|Q9VQJ8_DROME  50 LLNEEKIDK-EASRVLDLGCGNGMF--LVGLANEGFTgdLTGVDYSPKAVELAQNIAEDN-KLSITYKVADltQPQNELGQFDVVHD 132
                             444444444.689********7654..89999999854489*****************99.99998888751145544443445666 PP

                             -S-EEEE CS
         Methyltransf_10 175 nerydvt 181
                             +  yd++
  tr|Q9VQJ8|Q9VQJ8_DROME 133 KGTYDAV 139
                             6677754 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (299 nodes)
Target sequences:                        13786  (7394964 residues searched)
Passed MSV filter:                       407  (0.0295227); expected 275.7 (0.02)
Passed bias filter:                      289  (0.0209633); expected 275.7 (0.02)
Passed Vit filter:                        11  (0.000797911); expected 13.8 (0.001)
Passed Fwd filter:                         2  (0.000145075); 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.15u 0.01s 00:00:00.16 Elapsed: 00:00:00.05
# Mc/sec: 44221.88
//
[ok]