# 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/TRM.hmm.txt # target sequence database: proteomes/Macaca_mulatta.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TRM [M=376] Accession: PF02005.15 Description: N2,N2-dimethylguanosine tRNA methyltransferase 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.3e-144 480.9 0.0 4.2e-144 480.6 0.0 1.0 1 tr|F6RHQ5|F6RHQ5_MACMU tRNA (guanine(26)-N(2))-dimethyl 5.3e-66 223.6 2.9 1.6e-55 189.1 0.5 3.1 3 tr|F6RCB6|F6RCB6_MACMU TRMT1-like protein isoform 1 OS= ------ inclusion threshold ------ 0.022 14.2 0.0 0.034 13.6 0.0 1.3 1 sp|F7GSQ4|TRM5_MACMU tRNA (guanine(37)-N1)-methyltran 0.023 14.1 0.0 0.035 13.5 0.0 1.2 1 tr|A0A1I7Q472|A0A1I7Q472_MACMU tRNA (guanine(37)-N1)-methyltran Domain annotation for each sequence (and alignments): >> tr|F6RHQ5|F6RHQ5_MACMU tRNA (guanine(26)-N(2))-dimethyltransferase OS=Macaca mulatta OX=9544 GN=TRMT1 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 480.6 0.0 7.9e-148 4.2e-144 7 376 .] 53 488 .. 42 488 .. 0.91 Alignments for each domain: == domain 1 score: 480.6 bits; conditional E-value: 7.9e-148 E----.SS--TT--SS--GGGHHHHHHHHHH........HHHHHHHH-S-..................................... CS TRM 7 vedke.vtvskkeeVFYNPvqefnRDlsvlvi....r.qlellkkklgek..................................... 50 +e+++ + +++++eVFYNPvqefnRDl+++vi r ql +++ +++ + tr|F6RHQ5|F6RHQ5_MACMU 53 TEGAAkIAFPSANEVFYNPVQEFNRDLTCAVItefaRiQLGAKGIQIKVPgekdtqkvvvdlseqeeekvelkenedlasgdqprta 139 33333266***********************977766666777777766699*********************************** PP .........EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEEE CS TRM 51 .........ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdviD 125 ++vl++L+asGlR+iRfalEvpg ++v+aND s +Av+li++Nv+ln+v+++v+++++dA++l++ ++++rfdviD tr|F6RHQ5|F6RHQ5_MACMU 140 avgeiceegLHVLEGLAASGLRSIRFALEVPGLRSVVANDTSARAVDLIRRNVQLNDVAHLVQPSQADARMLMYqhqRVSERFDVID 226 **************************************************************************8877889****** PP E--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEE CS TRM 126 lDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsa 212 lDP+Gspapfld+avq+v+e+glL+vT+TD+avLaG+s e+c++kYga+alks++cHE+alRi+l++++ +a++Y+++v+plls+sa tr|F6RHQ5|F6RHQ5_MACMU 227 LDPYGSPAPFLDAAVQAVSEGGLLCVTCTDMAVLAGNSGETCYSKYGAMALKSRACHEMALRIVLHSLDLRANCYQRFVVPLLSISA 313 *************************************************************************************** PP TTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT.......--EEEE--GGG--.....SB-TTT-SBEEEEEEEE-S-S--HHHHHH CS TRM 213 dfYvrvfvkvktgaakvkkvleklglvlhCsgC.......lseekvkglakfk.....aecehcgaklhlaGPlwlgplkdkefvee 287 dfYvrvfv+v tg+akvk++++k++lv++C+gC l +k +g+ + + +l+GP+w++p++d +fv++ tr|F6RHQ5|F6RHQ5_MACMU 314 DFYVRVFVRVFTGQAKVKASASKQALVFQCVGCgayhlqrL--GKASGVPGSRvsevmVPAHLSPVLPQLGGPVWAEPIHDLDFVGR 398 *********************************88775331..455555444444442333334456799***************** PP HHHHHCT..-TTS..-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHH CS TRM 288 vleeaekekeeel..kkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavv 372 vle+++++ ++ ++r++++L++i+eEl+d+plyy+l+ql+s++++++p+l +++sal++++f+vs++Ha+++a+kT+ap++a++ tr|F6RHQ5|F6RHQ5_MACMU 399 VLEAVSANP-GRFhtSERIRGVLSVITEELPDVPLYYTLDQLSSTIHCNTPSLLQLRSALLHADFRVSLSHACKNAVKTDAPASALW 484 *******75.44445************************************************************************ PP HHHH CS TRM 373 evlr 376 +++r tr|F6RHQ5|F6RHQ5_MACMU 485 DIMR 488 ***7 PP >> tr|F6RCB6|F6RCB6_MACMU TRMT1-like protein isoform 1 OS=Macaca mulatta OX=9544 GN=TRMT1L PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.9 0.0 0.32 1.7e+03 230 246 .. 176 201 .. 171 223 .. 0.64 2 ! 189.1 0.5 3e-59 1.6e-55 12 294 .. 244 552 .. 228 591 .. 0.84 3 ! 33.0 0.0 8.1e-12 4.3e-08 312 374 .. 630 689 .. 593 691 .. 0.78 Alignments for each domain: == domain 1 score: -1.9 bits; conditional E-value: 0.32 HHHHTTEEEEEEETTT.........- CS TRM 230 kkvleklglvlhCsgC.........l 246 +++ +k+g +hC+ C + tr|F6RCB6|F6RCB6_MACMU 176 EQISSKMGAHYHCIICsatitrrtdM 201 678899**********6554333220 PP == domain 2 score: 189.1 bits; conditional E-value: 3e-59 SS--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT CS TRM 12 vtvskkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklne 98 +++k++ ++NP ++ nR l + l++l ++ ++++ lda+ a+G+ ++ +a+ + ++ +v +ND++e+ v li+eN ln+ tr|F6RCB6|F6RCB6_MACMU 244 YSIPQKTDSYFNPKMKLNRQLIFCT---LAALAEE-RKPLECLDAFGATGIMGLQWAKHLGNAVKVTINDLNENSVTLIQENCYLNK 326 3499999**************9888...5666655.589*******************************************88888 PP --......................GGEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSS CS TRM 99 ve......................eivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGas 163 + +++ ++dAnvl+ + + fd+i lDPfG+++++ldsa + +++ g+++vT+TD ++L ++ tr|F6RCB6|F6RCB6_MACMU 327 LKvvvdskekeksddileegekniGNIKVTKMDANVLM--HLRSFDFIHLDPFGTSVNYLDSAFRNIRNLGIVSVTSTDISSLYAKA 411 77778888888999999999886635899999******..7899******************************************* PP HHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT..... CS TRM 164 pekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgC..... 245 ++ r+Yg + +++e++ E a Ri++ avar+aa+ +k +e+l++++ +++v v v+v +g++ ++++++k+ +++hC C tr|F6RCB6|F6RCB6_MACMU 412 QHVARRHYGCNIVRTEYYKELAARIVVAAVARAAARCNKGIEVLFAVALEHFVLVVVRVLRGPTSADETAKKIQYLIHCQWCeerif 498 **********************************************************************************66654 PP .....--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT CS TRM 246 .....lseekvkglakfk........aecehcgaklhlaGPlwlgplkdkefveevleeaek 294 + + + + + g++ GPlw l ++ f++++l e + tr|F6RCB6|F6RCB6_MACMU 499 qkdgnM-----VEE---NpyrqlpcnCHGSMPGKTAIELGPLWSSSLFNTGFLKRMLFESLH 552 332222.....222...1333444213334456666778****************9866544 PP == domain 3 score: 33.0 bits; conditional E-value: 8.1e-12 HHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHH CS TRM 312 keElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvev 374 + E + p+yy++++ s+ +++p+lkk+ l ++gf+vsrtH+ p +++T+apl ++ + tr|F6RCB6|F6RCB6_MACMU 630 STE--HPPFYYNIHR-HSIKGMNMPKLKKFLCYLSQAGFRVSRTHFDPMGVRTDAPLMQFKSI 689 333..6788888887.56889***********************************9888766 PP >> sp|F7GSQ4|TRM5_MACMU tRNA (guanine(37)-N1)-methyltransferase OS=Macaca mulatta OX=9544 GN=TRMT5 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.6 0.0 6.3e-06 0.034 54 115 .. 303 361 .. 241 370 .. 0.80 Alignments for each domain: == domain 1 score: 13.6 bits; conditional E-value: 6.3e-06 EESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC CS TRM 54 ldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla 115 +d +++ G ai a+ +v+aND++++ + + +N kln+v++ v+i + d + +l+ sp|F7GSQ4|TRM5_MACMU 303 FDVFAGVGPFAIPVAKK---NCTVFANDLNPESHKWLLHNCKLNKVDQKVKIFNLDGKDFLQ 361 56666666666666655...568************************999998889888776 PP >> tr|A0A1I7Q472|A0A1I7Q472_MACMU tRNA (guanine(37)-N1)-methyltransferase OS=Macaca mulatta OX=9544 GN=TRMT5 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.0 6.7e-06 0.035 55 115 .. 283 340 .. 232 349 .. 0.81 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 6.7e-06 ESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC CS TRM 55 daLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla 115 d +++ G ai a+ +v+aND++++ + + +N kln+v++ v+i + d + +l+ tr|A0A1I7Q472|A0A1I7Q472_MACMU 283 DVFAGVGPFAIPVAKK---NCTVFANDLNPESHKWLLHNCKLNKVDQKVKIFNLDGKDFLQ 340 5556666656665555...568************************999998889888776 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 21211 (10888574 residues searched) Passed MSV filter: 682 (0.0321531); expected 424.2 (0.02) Passed bias filter: 466 (0.0219697); expected 424.2 (0.02) Passed Vit filter: 46 (0.00216869); expected 21.2 (0.001) Passed Fwd filter: 5 (0.000235727); expected 0.2 (1e-05) Initial search space (Z): 21211 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.32u 0.01s 00:00:00.33 Elapsed: 00:00:00.09 # Mc/sec: 45490.04 // [ok]