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

Query:       Gcd10p  [M=300]
Accession:   PF04189.12
Description: Gcd10p family
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
    ------- ------ -----    ------- ------ -----   ---- --  --------              -----------
     3e-102  340.2   4.2   3.7e-102  339.9   4.2    1.1  1  sp|P41814|TRM6_YEAST   tRNA (adenine(58)-N(1))-methyltransfer
     0.0024   15.5   0.1      0.007   14.0   0.0    1.7  2  sp|P46959|TRM61_YEAST  tRNA (adenine(58)-N(1))-methyltransfer


Domain annotation for each sequence (and alignments):
>> sp|P41814|TRM6_YEAST  tRNA (adenine(58)-N(1))-methyltransferase non-catalytic subunit TRM6 OS=Saccharomyces cerevisia
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  339.9   4.2  1.2e-105  3.7e-102       1     300 []       6     329 ..       6     329 .. 0.92

  Alignments for each domain:
  == domain 1  score: 339.9 bits;  conditional E-value: 1.2e-105
                Gcd10p   1 lIkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeks..............rlrvvpaaeleae 75 
                           +I  nqhv+++Lps+n+kiv++kpnt++sLgKfg+f++n+iig+p+gltfei+++ ee ++               ++r++++  + ++
  sp|P41814|TRM6_YEAST   6 TIDFNQHVIVRLPSKNYKIVELKPNTSVSLGKFGAFEVNDIIGYPFGLTFEIYYDGEEVSSDenrdskpknkipigKVRLLSQ--EIKD 92 
                           6999***************************************************9887777899999999888888888833..2222 PP

                Gcd10p  76 slaeeeaeeeeeeee...................ardnreiiddga.rQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSq 144
                           + +++++ ++    e                   +++n+++++ g+ +Q+Lt eeIe++K+e+ s +keiI+k+++sh++f++KT++Sq
  sp|P41814|TRM6_YEAST  93 VNNDKDDGQS----EpplsikeksvslelssidsSATNQNLVNMGSkAQELTVEEIEKMKQESLS-SKEIIDKIIKSHKSFHNKTVYSQ 176
                           2222222222....12345678899******************8888******************.*********************** PP

                Gcd10p 145 eKYlkrKkkKYlkrftvlpldvsllleyllekkdaqkilelreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgif.essake 232
                           eKY++rKk+K++k+ftv++l++s+ll++l++k d+q++l++++e++g+ll+lan++++g+yL++D+tgGLlv++l+eRm++  ++s+++
  sp|P41814|TRM6_YEAST 177 EKYVNRKKQKFAKYFTVEYLSSSNLLQFLIDKGDIQRVLDMSQESMGMLLNLANIQSEGNYLCMDETGGLLVYFLLERMFGGdNESKSK 265
                           ***************************************************************************************** PP

                Gcd10p 233 gtitliheneqpnlsllkyfnydaaepehplkkhlktlswlqllepeedetyeeepeevseeelaelk 300
                           g++++ihene++nl+llk+ ny    +e+++k+h++t+s+l+++ep + ++ +++++++++ee ++lk
  sp|P41814|TRM6_YEAST 266 GKVIVIHENEHANLDLLKFANY----SEKFIKEHVHTISLLDFFEPPTLQEIQSRFTPLPKEEARALK 329
                           **********************....9***********************************998876 PP

>> sp|P46959|TRM61_YEAST  tRNA (adenine(58)-N(1))-methyltransferase catalytic subunit TRM61 OS=Saccharomyces cerevisiae 
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !   14.0   0.0   2.3e-06     0.007       2      66 ..      13      75 ..      12      95 .. 0.84
   2 ?   -1.4   0.1      0.11   3.3e+02      30      76 ..     243     290 ..     235     340 .. 0.69

  Alignments for each domain:
  == domain 1  score: 14.0 bits;  conditional E-value: 2.3e-06
                 Gcd10p  2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrv 66
                           Ik+++  l+ + ++n+k v++  ++ +   ++gsfp ++iig+pyg++  i  k ++ ++  ++v
  sp|P46959|TRM61_YEAST 13 IKEGDLTLIWVSRDNIKPVRMHSEE-VFNTRYGSFPHKDIIGKPYGSQIAIRTKGSN-KFAFVHV 75
                           899999999999********99888.5568*****************9999988854.4554554 PP

  == domain 2  score: -1.4 bits;  conditional E-value: 0.11
                 Gcd10p  30 LgKfg..sfplnliigrpygltfeildkreeeeksrlrvvpaaeleaes 76 
                            L+K+g   +++ +i gr+y++  +++++ ++   +rlr ++ ++l+  +
  sp|P46959|TRM61_YEAST 243 LEKYGwtDVEMVEIQGRQYESRRQMVRSLNDA-LERLRDIKRHKLQGVE 290
                            66666223556688899999999998888664.4568888776655433 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (300 nodes)
Target sequences:                         6049  (2936363 residues searched)
Passed MSV filter:                       591  (0.0977021); expected 121.0 (0.02)
Passed bias filter:                      109  (0.0180195); expected 121.0 (0.02)
Passed Vit filter:                         6  (0.000991899); expected 6.0 (0.001)
Passed Fwd filter:                         2  (0.000330633); 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.07u 0.00s 00:00:00.07 Elapsed: 00:00:00.03
# Mc/sec: 29363.63
//
[ok]