# 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/Dictyostelium_discoideum.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
    ------- ------ -----    ------- ------ -----   ---- --  --------              -----------
    1.5e-59  201.3  11.5    6.1e-37  127.1   0.6    3.0  3  sp|Q554C9|MET16_DICDI  U6 small nuclear RNA (adenine-(43)-N(6


Domain annotation for each sequence (and alignments):
>> sp|Q554C9|MET16_DICDI  U6 small nuclear RNA (adenine-(43)-N(6))-methyltransferase OS=Dictyostelium discoideum OX=4468
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  127.1   0.6   4.7e-41   6.1e-37       3     189 ..      38     231 ..      36     239 .. 0.86
   2 !   81.4   0.3     4e-27     5e-23     175     299 .]     239     355 ..     234     355 .. 0.85
   3 ?   -2.4   0.5      0.13   1.7e+03     196     226 ..     515     540 ..     476     562 .. 0.58

  Alignments for each domain:
  == domain 1  score: 127.1 bits;  conditional E-value: 4.7e-41
                            XXXXXXXXXXXXXXX.XXXXXXXXXXXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS----HHHHHHHHHHH CS
        Methyltransf_10   3 lksglhprnrhkdky.dfaalikvvPelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPpvPgradyihyv 89 
                            +++++h  n ++ ++ +f  l++++P + +++i + +  + id++d  a k l++ ll +++g+++ ++Pd++l+P++ +r++y++++
  sp|Q554C9|MET16_DICDI  38 SNDSMHINNFYRYNPpNFKLLASKYPTFDKYIINKTEKIYNIDWKDSNATKELTRVLLDHDFGLRI-ELPDNYLCPTLTLRINYLYWI 124
                            5789********9988**************************************************.********************* PP

                            HH......HHT--TCGCS---EEEEES-TTTTHHHHHHHHH.H--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE--ST-SSTTT CS
        Methyltransf_10  90 ad......llgkadsdisaqvraldiGvGanCiyPllGvte.yGWrfvgsevddvslasakaivesnsklsdaiklrkqkekkaifdg 170
                             d      ++++++++ ++ ++++diG+G++Ci+PllG+++ + W+f+g ++dd+ l++a++++  n  l+++i l k ++++ i+  
  sp|Q554C9|MET16_DICDI 125 SDqlknlkIILNDNDNDNKIIKGIDIGTGTSCIFPLLGAKLfNNWSFIGIDIDDKVLEYAQNNITIN-SLNSKITLFKNEKNSDILLK 211
                            *93333222344555566699******************98368***********************.9******9998888888777 PP

                            STT--S.-EEEEEE-----S CS
        Methyltransf_10 171 lieene.rydvtlCnPPfha 189
                            l + +e   ++   n    +
  sp|Q554C9|MET16_DICDI 212 LLNYKEgSNTFNSSNDDHQD 231
                            77766646666666655555 PP

  == domain 2  score: 81.4 bits;  conditional E-value: 4e-27
                            -S-E..EEEEE-----SS------------...---.....---------TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEEEEESSG CS
        Methyltransf_10 175 nery..dvtlCnPPfhaslaeakagserklknlgkanktralpslnfgGqkaelyceGGevafikklieeskqlakqvlWfttlvskk 260
                            +e+y  d++lCnPPf+ +l+e ++ +             + +p++n +G  +e++++GGe +f+k++i+es ql+ +++++tt++++k
  sp|Q554C9|MET16_DICDI 239 DEEYfaDFCLCNPPFFKDLNENNNNK-------------NNNPKSNCTGSVNEMVTDGGEFEFVKRIIKESFQLKCKIRFYTTMIGRK 313
                            4444459*********9998753221.............224789******************************************* PP

                            GGHHHHHHHHHHT...T-SEEEEEEEEETTEEEEEEEEES-- CS
        Methyltransf_10 261 enlaplkeelkkl...GaskvtvvelaqGqklsrfiaWsfld 299
                             nl+pl + l k+     ++++++el+qG++ +++++W fl+
  sp|Q554C9|MET16_DICDI 314 VNLNPLINILIKQyylPKNQIQTTELVQGNTSRWVLSWYFLN 355
                            *******9997752225789*******************987 PP

  == domain 3  score: -2.4 bits;  conditional E-value: 0.13
                            -------...---.....---------TTTT CS
        Methyltransf_10 196 agserklknlgkanktralpslnfgGqkael 226
                            + ++ +     + n+++ +++ n+  ++++ 
  sp|Q554C9|MET16_DICDI 515 NNNNNN-----NNNNNNNNNNNNNNNKNNNN 540
                            111111.....11222222222222222333 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (299 nodes)
Target sequences:                        12739  (6856946 residues searched)
Passed MSV filter:                       933  (0.0732397); expected 254.8 (0.02)
Passed bias filter:                      197  (0.0154643); expected 254.8 (0.02)
Passed Vit filter:                        15  (0.00117749); expected 12.7 (0.001)
Passed Fwd filter:                         1  (7.84991e-05); expected 0.1 (1e-05)
Initial search space (Z):              12739  [actual number of targets]
Domain search space  (domZ):               1  [number of targets reported over threshold]
# CPU time: 0.14u 0.00s 00:00:00.14 Elapsed: 00:00:00.06
# Mc/sec: 34170.45
//
[ok]