# 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]