# 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/Pan_troglodytes.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.6e-152 506.0 0.0 1e-151 505.8 0.0 1.0 1 tr|H2R3Q5|H2R3Q5_PANTR tRNA (guanine(26)-N(2))-dimethyl 3.7e-75 253.8 0.1 1.6e-68 232.0 0.1 2.0 2 tr|A0A2I3S380|A0A2I3S380_PANTR tRNA (guanine(26)-N(2))-dimethyl 5.4e-66 223.7 3.4 1.1e-55 189.7 0.4 3.1 3 tr|G2HHV3|G2HHV3_PANTR N2,N2-dimethylguanosine tRNA met ------ inclusion threshold ------ 0.099 12.1 0.0 0.15 11.5 0.0 1.3 1 tr|H2Q8F0|H2Q8F0_PANTR tRNA (guanine(37)-N1)-methyltran Domain annotation for each sequence (and alignments): >> tr|H2R3Q5|H2R3Q5_PANTR tRNA (guanine(26)-N(2))-dimethyltransferase OS=Pan troglodytes OX=9598 GN=TRMT1 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 505.8 0.0 1.7e-155 1e-151 12 376 .] 57 500 .. 45 500 .. 0.94 Alignments for each domain: == domain 1 score: 505.8 bits; conditional E-value: 1.7e-155 SS.......--TT--SS--GGGHHHHHHHHHH........HHHHHHHH-S-.................................... CS TRM 12 vt.......vskkeeVFYNPvqefnRDlsvlvi....r.qlellkkklgek.................................... 50 vt +++++eVFYNPvqefnRDl+++vi r ql +++ +++ + tr|H2R3Q5|H2R3Q5_PANTR 57 VTegaakiaFPSANEVFYNPVQEFNRDLTCAVItefaRiQLGAKGIQIKVPgekdtqkvvvdlseqeeekvelkesenlasgdqprt 143 224455556***********************977766666777777766699********************************** PP ..........EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEE CS TRM 51 ..........ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdvi 124 ++vl++L+asGlR+iRfalEvpg ++v+aND+s++Av+li++Nv+ln+v+++v+++++dA++l++ ++++rfdvi tr|H2R3Q5|H2R3Q5_PANTR 144 aavgeiceegLHVLEGLAASGLRSIRFALEVPGLRSVVANDASTRAVDLIRRNVQLNDVAHLVQPSQADARMLMYqhqRVSERFDVI 230 ***************************************************************************8877889***** PP EE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEE CS TRM 125 DlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllsls 211 DlDP+Gspa+fld+avq+v+e+glL+vT+TD+avLaG+s e+c++kYga+alks++cHE+alRi+l++++ +a++Y+++v+plls+s tr|H2R3Q5|H2R3Q5_PANTR 231 DLDPYGSPATFLDAAVQAVSEGGLLCVTCTDMAVLAGNSGETCYSKYGAMALKSRACHEMALRIVLHSLDLRANCYQRFVVPLLSIS 317 *************************************************************************************** PP ETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT.........--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S-- CS TRM 212 adfYvrvfvkvktgaakvkkvleklglvlhCsgC.........lseekvkglakfk........aecehcgaklhlaGPlwlgplkd 281 adfYvrvfv+v tg+akvk++++k++lv++C+gC +++ +g+ kf+ ecehcg++++l+GP+w++p++d tr|H2R3Q5|H2R3Q5_PANTR 318 ADFYVRVFVRVFTGQAKVKASASKQALVFQCVGCgafhlqrlgKASGVPSGRVKFSaacgppvtPECEHCGQRHQLGGPMWAEPIHD 404 *************************************9888764448889999999******************************* PP HHHHHHHHHHHCT..-TTS..-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS CS TRM 282 kefveevleeaekekeeel..kkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTna 366 +fv++vle+++++ ++ ++r++++L++i+eEl+d+plyy+l+ql+s++++++p+l +++sal++++f+vs++Ha+++a+kT+a tr|H2R3Q5|H2R3Q5_PANTR 405 LDFVGRVLEAVSANP-GRFhtSERIRGVLSVITEELPDVPLYYTLDQLSSTIHCNTPSLLQLRSALLHADFRVSLSHACKNAVKTDA 490 *************75.44445****************************************************************** PP -HHHHHHHHH CS TRM 367 pleavvevlr 376 p++a+++++r tr|H2R3Q5|H2R3Q5_PANTR 491 PASALWDIMR 500 *********7 PP >> tr|A0A2I3S380|A0A2I3S380_PANTR tRNA (guanine(26)-N(2))-dimethyltransferase OS=Pan troglodytes OX=9598 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 19.8 0.0 8.3e-08 0.00048 50 75 .. 96 144 .. 60 148 .. 0.75 2 ! 232.0 0.1 2.8e-72 1.6e-68 105 365 .. 164 416 .. 156 423 .. 0.92 Alignments for each domain: == domain 1 score: 19.8 bits; conditional E-value: 8.3e-08 -.......................EEEEESS-TTSHHHHHHHHH-S-EE CS TRM 50 k.......................ikvldaLsasGlRaiRfalEvpgvk 75 k ++vl++L+asGl +iRfalEvpg++ tr|A0A2I3S380|A0A2I3S380_PANTR 96 KgsenlasgdqprtatvweiceegLHVLEGLAASGLLSIRFALEVPGAQ 144 266666666667777777777777*********************9875 PP == domain 2 score: 232.0 bits; conditional E-value: 2.8e-72 EE-S-HHHHHC...--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STT CS TRM 105 ineddAnvlla...ekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalksel 180 +++ d +l++ ++++rfdv lDP+Gspapfld+av +v+e+ lL+vT TDtavLa +s e c++k+ga+alks++ tr|A0A2I3S380|A0A2I3S380_PANTR 164 PSQVDGWMLMYqhqRVSERFDVTHLDPYGSPAPFLDAAVLAVSEGELLCVTRTDTAVLARNSGEMCYSKFGAMALKSRA 242 56789999**98777889************************************************************* PP HHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT--EEEE--GGG--S CS TRM 181 cHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgClseekvkglakfka 259 cHE+alRi+l++++ +a++Y+++v+plls+sadfYvrvfv+v tg+ kvk+++ +g +++ g+ k + tr|A0A2I3S380|A0A2I3S380_PANTR 243 CHEMALRIVLHSLDLRANCYQRFVVPLLSISADFYVRVFVRVFTGQPKVKASAMGCGSFHLQRLGKASGVPRGRVKSSV 321 ****************************************************999998766344455589999999885 PP B-.TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS-- CS TRM 260 ec.ehcgaklhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklikeElldaplyydlsqlasvlklsvpp 337 c ehcg++++l+ +w++p++d +fv+ vle+++++ + ++ ++ ++L++++eE l +l+s++ +++p+ tr|A0A2I3S380|A0A2I3S380_PANTR 322 ACrEHCGQRHQLGCSMWAEPIHDLDFVGCVLEAVSAN-PGRF--HTSEVLSVMTEE---------LPELSSTILCHTPS 388 5538*******************************94.6666..46799*******.........889*********** PP -HHHHHHH....TT-EE-TTSTTEEE-S CS TRM 338 lkkvvsalkelgfevsrtHakpkaikTn 365 l ++sal+++ f+vs++Ha+++a+kT tr|A0A2I3S380|A0A2I3S380_PANTR 389 LLPLRSALLHAHFRVSLSHACKNAVKTY 416 **************************95 PP >> tr|G2HHV3|G2HHV3_PANTR N2,N2-dimethylguanosine tRNA methyltransferase-like protein OS=Pan troglodytes OX=9598 GN=TRM # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.3 0.1 0.21 1.2e+03 230 246 .. 174 200 .. 168 223 .. 0.65 2 ! 189.7 0.4 1.9e-59 1.1e-55 12 292 .. 242 548 .. 226 571 .. 0.86 3 ! 32.9 0.0 8.7e-12 5e-08 313 374 .. 629 687 .. 590 689 .. 0.79 Alignments for each domain: == domain 1 score: -1.3 bits; conditional E-value: 0.21 HHHHTTEEEEEEETTT..........- CS TRM 230 kkvleklglvlhCsgC..........l 246 +++++k+g +hC+ C l tr|G2HHV3|G2HHV3_PANTR 174 EQITSKMGAHYHCIICsatitrrtdmL 200 78899***********66543333320 PP == domain 2 score: 189.7 bits; conditional E-value: 1.9e-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|G2HHV3|G2HHV3_PANTR 242 YSIPQKTDSYFNPKMKLNRQLIFCT---LAALAEE-RKPLECLDAFGATGIMGLQWAKHLGNAVKVTINDLNENSVTLIQENCHLNK 324 3499999**************9888...5666655.589*******************************************99999 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|G2HHV3|G2HHV3_PANTR 325 LKvvvdskekeksddileegeknlGNIKVTKMDANVLM--HLRSFDFIHLDPFGTSVNYLDSAFRNIRNLGIVSVTSTDISSLYAKA 409 8888999999999*999999987636999999******..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|G2HHV3|G2HHV3_PANTR 410 QHVARRHYGCNIVRTEYYKELAARIVVAAVARAAARCNKGIEVLFAVALEHFVLVVVRVLRGPTSADETAKKIQYLIHCQWCeerif 496 **********************************************************************************66654 PP .....--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHH CS TRM 246 .....lseekvkglakfk........aecehcgaklhlaGPlwlgplkdkefveevleea 292 + + + + + g++ GPlw l ++ f++++l e tr|G2HHV3|G2HHV3_PANTR 497 qkdgnM-----VEE---NpyrqlpcnCHGSMPGKTAIELGPLWSSSLFNTGFLKRMLFES 548 332222.....222...1333444213334456666778****************98654 PP == domain 3 score: 32.9 bits; conditional E-value: 8.7e-12 HHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHH CS TRM 313 eElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvev 374 + p+yy++++ s+ +++p+lkk+ l ++gf+vsrtH+ p +++T+apl ++ + tr|G2HHV3|G2HHV3_PANTR 629 T--EHPPFYYNIHR-HSIKGMNMPKLKKFLCYLSQAGFRVSRTHFDPMGVRTDAPLMQFKSI 687 2..36788888887.56889***********************************9888766 PP >> tr|H2Q8F0|H2Q8F0_PANTR tRNA (guanine(37)-N1)-methyltransferase OS=Pan troglodytes OX=9598 GN=TRMT5 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 2.6e-05 0.15 55 115 .. 303 360 .. 240 369 .. 0.80 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 2.6e-05 ESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC CS TRM 55 daLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla 115 d +++ G ai a+ +v+aND++++ + + N kln+v++ v++ + d + +l+ tr|H2Q8F0|H2Q8F0_PANTR 303 DVFAGVGPFAIPVAKK---NCTVFANDLNPESHKWLLYNCKLNKVDQKVKVFNLDGKDFLQ 360 5556666666666555...568************************988887778887775 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 740 (0.0321627); expected 460.2 (0.02) Passed bias filter: 516 (0.022427); expected 460.2 (0.02) Passed Vit filter: 44 (0.00191238); expected 23.0 (0.001) Passed Fwd filter: 4 (0.000173853); expected 0.2 (1e-05) Initial search space (Z): 23008 [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.12 # Mc/sec: 37419.07 // [ok]