# 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/Danio_rerio.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.6e-142 473.7 0.0 8.4e-142 473.3 0.0 1.1 1 tr|F1QGD9|F1QGD9_DANRE tRNA methyltransferase 1 OS=Dani 4.9e-68 230.5 0.0 3e-55 188.4 0.0 2.3 2 tr|A0A2R8QS20|A0A2R8QS20_DANRE tRNA (guanine(26)-N(2))-dimethyl ------ inclusion threshold ------ 0.023 14.4 0.0 0.036 13.7 0.0 1.3 1 sp|B8A5G9|TRM5_DANRE tRNA (guanine(37)-N1)-methyltran 0.11 12.2 0.0 0.14 11.8 0.0 1.2 1 sp|Q0P466|TYW2_DANRE tRNA wybutosine-synthesizing pro 0.29 10.8 2.4 0.41 10.3 2.4 1.1 1 tr|F1R8M8|F1R8M8_DANRE Ribosomal L1 domain-containing 1 Domain annotation for each sequence (and alignments): >> tr|F1QGD9|F1QGD9_DANRE tRNA methyltransferase 1 OS=Danio rerio OX=7955 GN=trmt1 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 473.3 0.0 1.6e-145 8.4e-142 1 376 [] 89 542 .. 89 542 .. 0.93 Alignments for each domain: == domain 1 score: 473.3 bits; conditional E-value: 1.6e-145 TTEEEEE----SS--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-..................................... CS TRM 1 GeakikvedkevtvskkeeVFYNPvqefnRDlsvlvirqlellkkklgek..................................... 50 Ge+ +++ ++++ +++++eVFYNPvqefnRDl+++vi e + +l+++ tr|F1QGD9|F1QGD9_DANRE 89 GETVVREGKAAILFPNANEVFYNPVQEFNRDLTCAVIT--EFARDQLAQRgvrilvpgekdrvvvqlsedkngagdqeapeqpeqeh 173 6677777788899***********************75..33344444444557888****************************** PP .......................EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHH CS TRM 51 .......................ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvll 114 ++vl++L+asGlR++RfalEvpg k++ aND s kA +li +N++ n+v +i+++++ dA++++ tr|F1QGD9|F1QGD9_DANRE 174 pkeeeaaekeyitatvgerceqgLRVLEGLAASGLRSVRFALEVPGLKSITANDYSAKAAALIARNAQHNNVTHILQASQRDASMVM 260 *************************************************************************************** PP C...--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHH CS TRM 115 a...ekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaa 198 + +k+r+dviDlDP+Gsp+pfldsavq+v+e+glL++T+TD+av+aG+s e+c++kYg+v++ks++cHE+alRi+l++++++aa tr|F1QGD9|F1QGD9_DANRE 261 YearGRKERYDVIDLDPYGSPSPFLDSAVQAVSEGGLLCITCTDMAVMAGNSGETCYSKYGSVSIKSKYCHEMALRIILHSLDQRAA 347 *888899******************************************************************************** PP TTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT...........--EEEE--GGG--..........SB-TTT CS TRM 199 kYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgC...........lseekvkglakfk..........aecehc 264 Y++++eplls+s+dfY+rvfv+vktg++ vk++++k++lv++C+gC +++ +++k +cehc tr|F1QGD9|F1QGD9_DANRE 348 VYQRYIEPLLSVSVDFYIRVFVRVKTGQVVVKNSASKQALVYNCVGCgafhlqrmgkrINQ-----GKHMKysaatgppvgPNCEHC 429 ************************************************9998775444222.....222223366688899****** PP -SBEEEEEEEE-S-S--HHHHHHHHHHHCT...-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....T CS TRM 265 gaklhlaGPlwlgplkdkefveevleeaekek.eeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgf 350 g++++l+GP w++pl+d +fv++vl +++ + + +++kr++++L++++eEl d+plyy+l++l+s+++ s+pp+ +++sal+++g+ tr|F1QGD9|F1QGD9_DANRE 430 GQRHQLGGPAWAEPLHDINFVQKVLGAVSGNPtRFRTSKRIEGMLSMVTEELQDVPLYYTLDHLSSTVHRSTPPMMQFRSALLHAGY 516 *****************************97534455************************************************** PP T-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 351 evsrtHakpkaikTnapleavvevlr 376 +vs +Ha+++a+kT+ap+ ++++++r tr|F1QGD9|F1QGD9_DANRE 517 RVSFSHACKNAVKTDAPAGVLWDIMR 542 *************************7 PP >> tr|A0A2R8QS20|A0A2R8QS20_DANRE tRNA (guanine(26)-N(2))-dimethyltransferase OS=Danio rerio OX=7955 GN=trmt1l PE=1 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 188.4 0.0 5.9e-59 3e-55 13 306 .. 173 490 .. 163 519 .. 0.89 2 ! 39.7 0.0 8.9e-14 4.6e-10 288 375 .. 530 618 .. 516 619 .. 0.81 Alignments for each domain: == domain 1 score: 188.4 bits; conditional E-value: 5.9e-59 S--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHH CS TRM 13 tvskkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelik 91 t ++k++ ++ +++nR l+ l++l ++ ++++ lda+ a+G+ ++ +a+ + + +v +NDi+e v++ik tr|A0A2R8QS20|A0A2R8QS20_DANRE 173 TTPQKTDTYFSRKMKTNRQLVFCS---LAVLAEE-RNPLECLDAFGATGIMGLQWAKHLRNSVKVTINDINEACVKMIK 247 677888999999********9988...7887777.699***************************************** PP HHHHHTT--....................GGEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEE CS TRM 92 eNvklneve....................eivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLl 150 eN +ln++ v++ ++dAnv++ + + fd+i lDPfG++ ++ld+a + v++ g+++ tr|A0A2R8QS20|A0A2R8QS20_DANRE 248 ENCHLNHIRvegnrtgrqldgagdpespiASVEVVKMDANVIM--HLRPFDYIHLDPFGTAINYLDAAFRNVRNLGIVS 324 *******9999**************9887679999********..8899****************************** PP EEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHH CS TRM 151 vTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakv 229 vT+TDt +L +s + lr+Yg + +++e++ E a R++l +var+aa+ +k +e+ll+++ +++v v v+v +g++++ tr|A0A2R8QS20|A0A2R8QS20_DANRE 325 VTSTDTGSLYSKSLNVTLRHYGCQIVRTEYYRELAARMVLASVARAAARCNKGIEVLLAVALEHFVLVVVRVLRGPTQA 403 ******************************************************************************* PP HHHHTTEEEEEEETTT..........--EEEE--GGG--SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-T CS TRM 230 kkvleklglvlhCsgC..........lseekvkglakfkaecehcgaklhlaGPlwlgplkdkefveevleeaekekee 298 +++++kl +lhC C + ee+ + ++ + + g++ GPlw gpl ++ f+ ++l +a +++ e tr|A0A2R8QS20|A0A2R8QS20_DANRE 404 DESAKKLRQLLHCQWCeervflkpgsMVEENLYRQLPCQCHGSMPGKTAVELGPLWGGPLFNTGFLRRMLVAAVQHNME 482 ****************777775554433344444444445556667778889******************999996666 PP TS-HHHHH CS TRM 299 elkkrvkk 306 ++ vk tr|A0A2R8QS20|A0A2R8QS20_DANRE 483 DIQPLVKT 490 65555554 PP == domain 2 score: 39.7 bits; conditional E-value: 8.9e-14 HHHHHCT..-TTS-HHHHHHHHHHHHHH.S--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-S CS TRM 288 vleeaekekeeelkkrvkklLklikeEl.ldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTn 365 ++ k++++ + ++Lk++k ++ l++p +y + s+ +++p+l+k+ + l e+gf+vsrtH+ p++++T+ tr|A0A2R8QS20|A0A2R8QS20_DANRE 530 ADTSDQSGKRKTTGDDGGNVLKKVKPDTsLEHPAFYYSIHRHSIRGMNMPKLNKFLQYLTEAGFRVSRTHFDPTGVRTD 608 44444444455556667789999998876899988888888999*********************************** PP S-HHHHHHHH CS TRM 366 apleavvevl 375 a+l ++ evl tr|A0A2R8QS20|A0A2R8QS20_DANRE 609 ATLAQFKEVL 618 ***9999987 PP >> sp|B8A5G9|TRM5_DANRE tRNA (guanine(37)-N1)-methyltransferase OS=Danio rerio OX=7955 GN=trmt5 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.7 0.0 7.1e-06 0.036 53 114 .. 286 344 .. 229 352 .. 0.79 Alignments for each domain: == domain 1 score: 13.7 bits; conditional E-value: 7.1e-06 EEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHH CS TRM 53 vldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvll 114 v+d +++ G ai a +vlaND++++ +++N+kln+v++ +++++ d + ++ sp|B8A5G9|TRM5_DANRE 286 VVDVFAGVGPFAIPAARR---GCAVLANDLNPESFRWLQHNAKLNKVDQKITTSNLDGRDFI 344 566666666666666655...5579***********************98888877877666 PP >> sp|Q0P466|TYW2_DANRE tRNA wybutosine-synthesizing protein 2 homolog OS=Danio rerio OX=7955 GN=trmt12 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.8 0.0 2.8e-05 0.14 52 112 .. 220 279 .. 132 288 .. 0.88 Alignments for each domain: == domain 1 score: 11.8 bits; conditional E-value: 2.8e-05 EEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHH CS TRM 52 kvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnv 112 +v+d ++ G ++ +++ +++ +v a + ++ Av+ +++N ++n v+++ +++++d + sp|Q0P466|TYW2_DANRE 220 TVVDLYAGIGYFTLPYLVH-ANAAHVHACEWNPDAVAALQRNLEINGVSNRCTVHQGDNRQ 279 5777888888888889888.89999*************************99999999765 PP >> tr|F1R8M8|F1R8M8_DANRE Ribosomal L1 domain-containing 1 OS=Danio rerio OX=7955 GN=rsl1d1 PE=1 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.3 2.4 7.9e-05 0.41 222 371 .. 206 351 .. 200 355 .. 0.70 Alignments for each domain: == domain 1 score: 10.3 bits; conditional E-value: 7.9e-05 EEE-HHHHHHHHTTEEEEEEETTT--EEEE--GGG--SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT...-TTS..-HHHH CS TRM 222 vktgaakvkkvleklglvlhCsgClseekvkglakfkaecehcgaklhlaGPlwlgplkdkefveevleeaekek.eeel..kkrvk 305 v ++++ +++++e++ + +s k ++ k + + +h+++++++a P++ +l ++vee ++a+ +k +++ k + tr|F1R8M8|F1R8M8_DANRE 206 VAHSEMTADQIVENVISAVS---TIS-AKLASTGK-NIKIIHLKSQTSVALPIYTSDLSHLALVEEAQKKARLTKgAQKRkrKTENT 287 55555666666666544443...233.33444444.4899****************************9998877522221233334 PP HHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHH CS TRM 306 klLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleav 371 + k +++E+ + +++ + + +++p++k+ +++l + g + + a ++++k + +++ tr|F1R8M8|F1R8M8_DANRE 288 EAPKPVTKETE--EEKDSEEEIPQLVPIETPKKKTKREGLSKKGLKKVPKPAAGNSVKKSVKDKKL 351 44556666633..33447889999999999999999999999999999999999988877666655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 25687 (13798380 residues searched) Passed MSV filter: 906 (0.0352708); expected 513.7 (0.02) Passed bias filter: 576 (0.0224238); expected 513.7 (0.02) Passed Vit filter: 44 (0.00171293); expected 25.7 (0.001) Passed Fwd filter: 5 (0.000194651); expected 0.3 (1e-05) Initial search space (Z): 25687 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.37u 0.01s 00:00:00.38 Elapsed: 00:00:00.11 # Mc/sec: 47165.37 // [ok]