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

Query:       WBS_methylT  [M=82]
Accession:   PF12589.7
Description: Methyltransferase involved in Williams-Beuren syndrome
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.6e-14   52.2   3.7    5.7e-14   51.1   3.7    1.6  1  sp|P25627|BUD23_YEAST  18S rRNA (guanine(1575)-N(7))-methyltr
  ------ inclusion threshold ------
       0.75    9.1   7.8        1.8    7.9   7.8    1.7  1  sp|P14063|RM31_YEAST   54S ribosomal protein L31, mitochondri


Domain annotation for each sequence (and alignments):
>> sp|P25627|BUD23_YEAST  18S rRNA (guanine(1575)-N(7))-methyltransferase OS=Saccharomyces cerevisiae (strain ATCC 20450
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !   51.1   3.7   1.9e-17   5.7e-14      25      81 ..     219     272 ..     198     273 .. 0.68

  Alignments for each domain:
  == domain 1  score: 51.1 bits;  conditional E-value: 1.9e-17
            WBS_methylT  25 eaeqvkyekrrrkrkkrkkkkkkavksreWilkKKerrrkkGkeVkrdSKYTGRkRk 81 
                            ++e+v+ +k+ r+r    k++k++++++ +il+KKe+++++G++V +dSK+TGRkR+
  sp|P25627|BUD23_YEAST 219 DEENVNLKKQLRQRL---KGGKDKESAKSFILRKKELMKRRGRKVAKDSKFTGRKRR 272
                            334444444443333...344444569*****************************8 PP

>> sp|P14063|RM31_YEAST  54S ribosomal protein L31, mitochondrial OS=Saccharomyces cerevisiae (strain ATCC 204508 / S288
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?    7.9   7.8   0.00058       1.8       8      79 ..      52     126 ..      22     129 .. 0.67

  Alignments for each domain:
  == domain 1  score: 7.9 bits;  conditional E-value: 0.00058
           WBS_methylT   8 qqqlpkglgedgeeeseeaeqvkyekrrrkrkkr..kkkkkkavksreWilkKKerrrkkGk.eVkrdSKYTGRk 79 
                           q    kgl+ ++++es+++ + + ++ r  +k +   k+++ + k++ W + KK    +kG  +V++ +K + Rk
  sp|P14063|RM31_YEAST  52 QRCQDKGLTYQEAMESKKKYKPRSKSLRLLNKPSvfPKENQMSSKDKYWTFDKKAVGYRKGIhKVPKWTKISIRK 126
                           4444455555554444422222222222223333343344444559***************99*******99998 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (82 nodes)
Target sequences:                         6049  (2936363 residues searched)
Passed MSV filter:                       695  (0.114895); expected 121.0 (0.02)
Passed bias filter:                      209  (0.0345512); expected 121.0 (0.02)
Passed Vit filter:                        21  (0.00347165); expected 6.0 (0.001)
Passed Fwd filter:                         3  (0.00049595); expected 0.1 (1e-05)
Initial search space (Z):               6049  [actual number of targets]
Domain search space  (domZ):               2  [number of targets reported over threshold]
# CPU time: 0.05u 0.00s 00:00:00.05 Elapsed: 00:00:00.02
# Mc/sec: 12039.09
//
[ok]