# 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/Tetrahymena_thermophila.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.9e-64 218.4 9.6 1.4e-62 212.3 4.5 3.2 2 tr|Q231G1|Q231G1_TETTS U6 small nuclear RNA (adenine-(43)-N( 2.9e-05 24.1 0.0 4e-05 23.7 0.0 1.1 1 tr|W7XCF3|W7XCF3_TETTS HemK methyltransferase family protein ------ inclusion threshold ------ 0.05 13.5 0.0 0.081 12.8 0.0 1.3 1 tr|I7MAE9|I7MAE9_TETTS UbiE/COQ5 family methyltransferase OS Domain annotation for each sequence (and alignments): >> tr|Q231G1|Q231G1_TETTS U6 small nuclear RNA (adenine-(43)-N(6))-methyltransferase OS=Tetrahymena thermophila (strain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 212.3 4.5 1.6e-66 1.4e-62 14 278 .. 71 317 .. 58 322 .. 0.91 2 ? 2.7 0.2 0.011 1e+02 271 297 .. 424 450 .. 343 451 .] 0.78 Alignments for each domain: == domain 1 score: 212.3 bits; conditional E-value: 1.6e-66 XXXXXXXXXXXXXXXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS----HHHHHHHHHHHHHHHT..--TC CS Methyltransf_10 14 kdkydfaalikvvPelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPpvPgradyihyvadllg..kads 98 k + df +l+ +P l+ +v +n +Gr++id+ + +a k ln+a+l+ ++++ +wd+P++ liP++P r++y+ ++ + l k ++ tr|Q231G1|Q231G1_TETTS 71 KFQWDFDELAIHYPPLNDYVFINERGRTTIDWTNTAAHKTLNQAILAVQFKIVYWDMPQDQLIPTIPSRLNYLEWIHNKLQslKINN 157 4467*************************************************************************9996225678 PP GCS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE--ST-SSTTTSTT--S.-EEEEEE- CS Methyltransf_10 99 disaqvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaivesnsklsdaiklrkqkekkaifdglieene.rydvtlCn 184 + s+q++a+diG+Ga+ iyP+lG+ +y W+f+++e++d s ++ak+i+++n++l+++++lrk + +if+ +i++++ +t+Cn tr|Q231G1|Q231G1_TETTS 158 EQSNQIKAIDIGTGASVIYPILGASKYNWSFIATELSDASYQNAKKIINNNPQLKNKVELRKS--QGNIFKDIIQKSDgDISFTMCN 242 8999*******************************************************9975..689*******99989******* PP ----SS------------...---.....---------TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEEEEESSGGGHHHHHHHHH CS Methyltransf_10 185 PPfhaslaeakagserklknlgkanktralpslnfgGqkaelyceGGevafikklieeskqlakqvlWfttlvskkenlaplkeelk 271 PPf+++ +e + + n +++e +ceGGe++f+k++ eesk +a+q+ ft+lv+k++n++++ + lk tr|Q231G1|Q231G1_TETTS 243 PPFFSEHEEKQI-----WVN--------------AQMKENEYICEGGELQFLKQMFEESKIFANQIKIFTSLVGKQTNFEQFLQYLK 310 ****98776443.....322..............23356899********************************************* PP HTT-SEE CS Methyltransf_10 272 klGaskv 278 ++ ++++ tr|Q231G1|Q231G1_TETTS 311 QQNIQNL 317 9998876 PP == domain 2 score: 2.7 bits; conditional E-value: 0.011 HHTT-SEEEEEEEEETTEEEEEEEEES CS Methyltransf_10 271 kklGaskvtvvelaqGqklsrfiaWsf 297 +k+ +++ ++qG++ ++++ W f tr|Q231G1|Q231G1_TETTS 424 NKQVSVQIFEGSIVQGNTTRWIVGWRF 450 3455567888889************98 PP >> tr|W7XCF3|W7XCF3_TETTS HemK methyltransferase family protein OS=Tetrahymena thermophila (strain SB210) OX=312017 GN= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 23.7 0.0 4.4e-09 4e-05 56 193 .. 68 208 .. 59 235 .. 0.78 Alignments for each domain: == domain 1 score: 23.7 bits; conditional E-value: 4.4e-09 HHHHHHHH----.---TTS-..---HHHHHHHHHHHHHHHT--TCGCS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHH CS Methyltransf_10 56 kallrefygvsewdiPdgfl..iPpvPgradyihyvadllgkadsdisaqvraldiGvGanCiyPllGvteyGWrfvgsevddvsla 140 +a++ +y + ew+ d l Ppv + y+ +l+ k ++++ ++ + l+iG+G+ i l +++ ++ ++ + ++ +++ tr|W7XCF3|W7XCF3_TETTS 68 RAMMPIQYVIGEWEFRDLLLkmKPPVLIPRNETSYLVELVNKLSKQK-QKCQFLEIGIGTGAISLSLLKENSQFTGIAIDKQKFCIE 153 6788889999999999876633699999999999*******999998.7788899998888887656555555666666666779** PP HHHHHHHHT.T-TTTEEEEE--ST-SSTTTSTT--S.-EEEEEE-----S.S--- CS Methyltransf_10 141 sakaivesnsklsdaiklrkqkekkaifdglieene.rydvtlCnPPfha.slae 193 ak++++ n + ++l + + k + l ++ + ++d+ + nPP+ s+ + tr|W7XCF3|W7XCF3_TETTS 154 LAKQNLNLNKIDPQRLQLIHLECLKFFQNNLNQNPDqQFDFIVSNPPYIPtSQVQ 208 *********6666889999888888888888877555**********97414444 PP >> tr|I7MAE9|I7MAE9_TETTS UbiE/COQ5 family methyltransferase OS=Tetrahymena thermophila (strain SB210) OX=312017 GN=TTH # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.8 0.0 9e-06 0.081 97 154 .. 46 104 .. 31 130 .. 0.77 Alignments for each domain: == domain 1 score: 12.8 bits; conditional E-value: 9e-06 TCGCS---EEEEES-TT.TTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TT CS Methyltransf_10 97 dsdisaqvraldiGvGa.nCiyPllGvteyGWrfvgsevddvslasakaivesnsklsd 154 +++ + +r+ld+GvG+ +y ++ + +g ++d + ak+i ++ s+++ tr|I7MAE9|I7MAE9_TETTS 46 KGNEKNFKRILDVGVGTgHPLYSIIDRIPKSVQVTGIDIDTNYIPAAKKIFSKHSNVEI 104 44455679********8357*******************************99887653 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (299 nodes) Target sequences: 26971 (16912921 residues searched) Passed MSV filter: 1867 (0.0692225); expected 539.4 (0.02) Passed bias filter: 292 (0.0108264); expected 539.4 (0.02) Passed Vit filter: 23 (0.000852768); expected 27.0 (0.001) Passed Fwd filter: 3 (0.000111231); expected 0.3 (1e-05) Initial search space (Z): 26971 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.31u 0.01s 00:00:00.32 Elapsed: 00:00:00.11 # Mc/sec: 45972.39 // [ok]