# 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/Mus_musculus.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.5e-148 494.1 0.0 4.5e-148 493.7 0.0 1.0 1 sp|Q3TX08|TRM1_MOUSE tRNA (guanine(26)-N(2))-dimethyltransf 1.5e-66 225.4 1.7 5.1e-56 190.8 0.1 2.8 3 sp|A2RSY6|TRM1L_MOUSE TRMT1-like protein OS=Mus musculus OX= ------ inclusion threshold ------ 0.042 13.3 0.0 0.066 12.7 0.0 1.3 1 sp|Q9D0C4|TRM5_MOUSE tRNA (guanine(37)-N1)-methyltransferas 0.079 12.4 0.0 0.11 12.0 0.0 1.1 1 sp|Q80ZM3|ETKMT_MOUSE Electron transfer flavoprotein beta su Domain annotation for each sequence (and alignments): >> sp|Q3TX08|TRM1_MOUSE tRNA (guanine(26)-N(2))-dimethyltransferase OS=Mus musculus OX=10090 GN=Trmt1 PE=1 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 493.7 0.0 8.1e-152 4.5e-148 9 376 .] 62 499 .. 55 499 .. 0.94 Alignments for each domain: == domain 1 score: 493.7 bits; conditional E-value: 8.1e-152 ---SS--TT--SS--GGGHHHHHHHHHH........HHHHHHHH-S-.......................................... CS TRM 9 dkevtvskkeeVFYNPvqefnRDlsvlvi....r.qlellkkklgek.......................................... 50 +++++++++eVFYNPvqefnRDl+++vi r +l +++ +++ + sp|Q3TX08|TRM1_MOUSE 62 AAKIVFPSANEVFYNPVQEFNRDLTCAVItefaRiHLGAKGIQIKVPgekdsekiavdlsdqeeetagknenlapgdwprtaavgeice 150 34488***********************965554555556666655599**************************************** PP ..EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEEEE--SS--GG CS TRM 51 ..ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdviDlDPfGspap 134 ++vl++L+asGlR+iRfalEvpg ++v+aND+s +Avel+++Nv+ln v+++v++n++dA++l++ + +rfdviDlDP+Gspap sp|Q3TX08|TRM1_MOUSE 151 egLRVLEGLAASGLRSIRFALEVPGLQSVVANDASARAVELMHRNVELNGVAHLVQPNQADARMLMYqhqKAPERFDVIDLDPYGSPAP 239 *******************************************************************66566689************** PP GHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEE CS TRM 135 fldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvk 223 fld+avq+v+++glL+vT+TD+avLaG+s e+c++kYga+alks++cHE+alRi+l++++ +a++Y+++++plls+sadfY+rvfv+v sp|Q3TX08|TRM1_MOUSE 240 FLDAAVQAVSDGGLLCVTCTDMAVLAGNSGETCYSKYGAMALKSRACHEMALRIVLHSLDLHANCYQRYIVPLLSISADFYIRVFVRVF 328 ***************************************************************************************** PP E-HHHHHHHHTTEEEEEEETTT.........--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT. CS TRM 224 tgaakvkkvleklglvlhCsgC.........lseekvkglakfk........aecehcgaklhlaGPlwlgplkdkefveevleeaeke 295 tg+akvk +++k++lv++C+gC ++ g+ kf+ ecehcg++++l+GP+w++p++d +fv++vl+++ + sp|Q3TX08|TRM1_MOUSE 329 TGQAKVKSSASKQALVFQCVGCgafylqrlgKASGDPGGRIKFSaacgppvtPECEHCGQRHQLGGPMWAEPIHDLDFVGRVLDAVTTN 417 **************************98876444888899999999999999***********************************97 PP .-TTS..-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 296 keeel..kkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvevlr 376 ++ + r++++L++++eEl+d+plyy+l+ql+s++++++p+l +++sal+++gf+vs++Ha+++a+kT+ap ea+++++r sp|Q3TX08|TRM1_MOUSE 418 P-GRFhtSMRIQGVLSVVTEELPDVPLYYTLDQLSSTIHCNTPRLLQLRSALLHAGFRVSLSHACKNAVKTDAPPEALWDIMR 499 5.4444599*************************************************************************7 PP >> sp|A2RSY6|TRM1L_MOUSE TRMT1-like protein OS=Mus musculus OX=10090 GN=Trmt1l PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.9 0.0 0.33 1.8e+03 229 253 .. 170 198 .. 165 219 .. 0.63 2 ! 190.8 0.1 9.2e-60 5.1e-56 13 291 .. 236 542 .. 220 563 .. 0.87 3 ! 32.8 0.0 9e-12 5e-08 316 374 .. 625 682 .. 595 684 .. 0.83 Alignments for each domain: == domain 1 score: -1.9 bits; conditional E-value: 0.33 HHHHHTTEEEEEEETTT.........--EEEE-- CS TRM 229 vkkvleklglvlhCsgC.........lseekvkg 253 ++++ +klg +hC+ C + g sp|A2RSY6|TRM1L_MOUSE 170 GEQISSKLGAHYHCIICsatitrrtdM-----LG 198 5678899**********6655333220.....22 PP == domain 2 score: 190.8 bits; conditional E-value: 9.2e-60 S....--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHH CS TRM 13 t....vskkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvkl 96 +++k++ ++NP ++ nR + + l++l k+ ++++ lda+ a+G+ ++ +a+ + ++ +v +ND++e+ v li++N +l sp|A2RSY6|TRM1L_MOUSE 236 FpnysIPQKTDSYFNPKMKLNRQIIFCT---LAALAKE-RKPLECLDAFGATGIMGLQWAKHLGNAVKVTINDLNENSVTLIQKNCHL 319 144459999**************98888...6666655.589*******************************************999 PP TT--.....................GGEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSS CS TRM 97 neve.....................eivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGas 163 n+ + + +++ ++dAnvl+ + + fd+i lDPfG+++++ldsa + v++ g+++vT+TD ++L ++ sp|A2RSY6|TRM1L_MOUSE 320 NKLKvvvdseekeegdaleddgtlgD-IQVTRMDANVLM--HLRSFDFIHLDPFGTSVNYLDSAFRNVRNLGIVSVTSTDISSLYAKA 404 99988999999********9998855.888889******..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 sp|A2RSY6|TRM1L_MOUSE 405 QHVARRHYGCNIVRTEYYKELAARIVVAAVARAAARCNKGIEVLFAVALEHFVLVVVRVLRGPTSADETAKKIQYLIHCQWCeerifq 492 **********************************************************************************666543 PP ....--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHH CS TRM 246 ....lseekvkglakfk........aecehcgaklhlaGPlwlgplkdkefveevlee 291 + + + + + g++ GPlw l ++ f++++l e sp|A2RSY6|TRM1L_MOUSE 493 kdgnM-----VEE---NpyrqlpcnCHGSMPGKTAIELGPLWSSSLFNTGFLKRMLFE 542 32222.....222...1333444213334456666778****************9855 PP == domain 3 score: 32.8 bits; conditional E-value: 9e-12 S--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHH CS TRM 316 ldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvev 374 + p+yy++++ s+ +++p+lkk+ l ++gf+vsrtH+ p +i+T+apl ++ + sp|A2RSY6|TRM1L_MOUSE 625 AHPPFYYNIHRH-SIKGMNMPKLKKFLCCLSQAGFRVSRTHFDPMGIRTDAPLMQFKSI 682 667899999885.6789***********************************9888766 PP >> sp|Q9D0C4|TRM5_MOUSE tRNA (guanine(37)-N1)-methyltransferase OS=Mus musculus OX=10090 GN=Trmt5 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.7 0.0 1.2e-05 0.066 74 114 .. 313 353 .. 256 360 .. 0.85 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 1.2e-05 EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHH CS TRM 74 vkevlaNDisekAvelikeNvklneveeivvineddAnvll 114 +v+aND++++ + + +N kln+v++ v++ ++d + ++ sp|Q9D0C4|TRM5_MOUSE 313 NCTVFANDLNPESHKWLLHNCKLNKVDQKVKVFNMDGKDFI 353 5689***********************99999998987766 PP >> sp|Q80ZM3|ETKMT_MOUSE Electron transfer flavoprotein beta subunit lysine methyltransferase OS=Mus musculus OX=10090 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.0 0.0 1.9e-05 0.11 53 124 .. 112 177 .. 93 182 .. 0.75 Alignments for each domain: == domain 1 score: 12.0 bits; conditional E-value: 1.9e-05 EEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC--SS-EEEE CS TRM 53 vldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvllaekkkrfdvi 124 vld s+ G ai a ++g++++laNDi + A i N kln + v + +l++++ +fd i sp|Q80ZM3|ETKMT_MOUSE 112 VLDLGSGCGATAI--AAKMSGASKILANDIDPIAGMAITLNCKLNGLNPFPVLTK----NILNTQQGKFDLI 177 5555555555555..77789*************************9986333333....3555666666665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 22289 (11757923 residues searched) Passed MSV filter: 733 (0.0328862); expected 445.8 (0.02) Passed bias filter: 519 (0.023285); expected 445.8 (0.02) Passed Vit filter: 38 (0.00170488); expected 22.3 (0.001) Passed Fwd filter: 4 (0.000179461); expected 0.2 (1e-05) Initial search space (Z): 22289 [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.10 # Mc/sec: 44209.79 // [ok]