# 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/Trichomonas_vaginalis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.5e-107 360.8 0.0 2.9e-107 360.6 0.0 1.0 1 tr|A2FQY9|A2FQY9_TRIVA tRNA (guanine(26)-N(2))-dimethyltrans 1.4e-104 351.7 0.0 1.8e-104 351.4 0.0 1.0 1 tr|A2FDH4|A2FDH4_TRIVA tRNA (guanine(26)-N(2))-dimethyltrans 1e-54 187.6 0.4 1.2e-54 187.4 0.4 1.0 1 tr|A2FGC5|A2FGC5_TRIVA Uncharacterized protein OS=Trichomona 1.2e-23 85.4 0.0 1.5e-23 85.0 0.0 1.1 1 tr|A2FGC4|A2FGC4_TRIVA Uncharacterized protein OS=Trichomona 2e-08 35.3 0.0 2.9e-08 34.8 0.0 1.1 1 sp|A2E5K9|TRM5_TRIVA tRNA (guanine(37)-N1)-methyltransfera ------ inclusion threshold ------ 0.038 14.6 0.1 0.044 14.4 0.1 1.1 1 tr|A2D9B5|A2D9B5_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2FQY9|A2FQY9_TRIVA tRNA (guanine(26)-N(2))-dimethyltransferase OS=Trichomonas vaginalis OX=5722 GN=TVAG_439220 P # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 360.6 0.0 3.4e-111 2.9e-107 15 376 .] 23 396 .. 4 396 .. 0.94 Alignments for each domain: == domain 1 score: 360.6 bits; conditional E-value: 3.4e-111 -TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--G CS TRM 15 skkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklnevee 101 +k+eVFYNPvq+fnRD+++lv l++++++++++ ++++a++a+GlR+iR+a Ev+g +++ aNDi ++A +i +N++ln+v++ tr|A2FQY9|A2FQY9_TRIVA 23 IDKSEVFYNPVQNFNRDITCLV---LQAYNNRVKKPTNLFEAFAATGLRSIRYAREVTGFNSITANDIDKNAGSIIDRNIQLNKVDN 106 4679******************...9999999999**************************************************** PP GEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHH CS TRM 102 ivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRi 188 iv+++ +dAn +++++k ++++iDlDP+ +++ fld+a++++++++lL+vT+TD +L+G+++e + Y+ + l++e+cHE+++R tr|A2FQY9|A2FQY9_TRIVA 107 IVKSSLGDANKVMNSHKGEYQIIDLDPYSTVSIFLDAAIHAAADGALLCVTSTDGRTLCGQQQEVSYGWYNVMVLNTEFCHEIGIRT 193 *************************************************************************************** PP HHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT...........--EEEE--GGG--.....S CS TRM 189 llkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgC...........lseekvkglakfk.....a 259 ll+++++ aa+Y+k++epl+sl+a+fY+rvfv v+ + +++k++++ +++l+ +s++++ +++ + tr|A2FQY9|A2FQY9_TRIVA 194 LLTLIKQIAARYKKTIEPLISLAANFYFRVFVIVHGKSNETKHLASDNAYLLFSPDShlfwlqplgrdMSKGENITVKPAQvtipsF 280 ****************************************************9777799**99999886665555555555667668 PP B-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH. CS TRM 260 ecehcgaklhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalk 346 + ++ +++l+l+GP++ gp+++kef++++++ + k +++ r++++L+ + E ld+p+yyd+s l+ ++++s+p+ + +vs l+ tr|A2FQY9|A2FQY9_TRIVA 281 KDPYTESQLQLGGPIYTGPIHNKEFIQSLIDLLPTMKFFKTTARIEAVLNSCIRE-LDSPFYYDISVLTGIIHSSTPSREILVSILE 366 999********************************988899**************.******************************* PP ...TT-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 347 elgfevsrtHakpkaikTnapleavvevlr 376 +lg+++s+tH+kp +kT+ap e ++++lr tr|A2FQY9|A2FQY9_TRIVA 367 KLGYKTSGTHCKPGMLKTDAPGELIWDILR 396 ****************************97 PP >> tr|A2FDH4|A2FDH4_TRIVA tRNA (guanine(26)-N(2))-dimethyltransferase OS=Trichomonas vaginalis OX=5722 GN=TVAG_353820 P # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 351.4 0.0 2.2e-108 1.8e-104 16 376 .] 30 402 .. 17 402 .. 0.93 Alignments for each domain: == domain 1 score: 351.4 bits; conditional E-value: 2.2e-108 TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GG CS TRM 16 kkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveei 102 ++ +VFYNPvq+fnRDl+++v ++++++ + +i++++a++asGlR+iR+a+E++g ke++ ND+ + Ave+ik+N+ +n+v++i tr|A2FDH4|A2FDH4_TRIVA 30 AEGDVFYNPVQQFNRDLTCAV---IQAYRNSVRDDISIFEAFAASGLRSIRYAKEITGIKEIISNDLDPGAVEIIKRNIAINKVSNI 113 5679*****************...67777777888**************************************************** PP EEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHH CS TRM 103 vvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRil 189 v+++++dA+ + e++++++v+DlDP+ ++apfl++av+++++++lL++T+TD +L+G +p+ ++ Yg v+l+ e+ HE+++R l tr|A2FDH4|A2FDH4_TRIVA 114 VTASQGDARDKMLEEENKYQVLDLDPYSTAAPFLEAAVKAATNGALLCITSTDGRTLCGVQPDVAYSWYGCVPLNVEFEHEFGIRCL 200 *************************************************************************************** PP HHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT............--EEEE--GGG--...... CS TRM 190 lkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgC............lseekvkglakfk...... 258 l+++++ a++ ++ +epl++l+a+fY+r+fv+++++++++kk ++ + ++++ + + k tr|A2FDH4|A2FDH4_TRIVA 201 LTTLMNVASRLKRSIEPLICLAANFYIRIFVRIHDKPIEAKKCASYTSIMTYSTDSksfwiqplgtctM-----------Kstnaiv 276 *********************************************999999977776676655544331...........2233444 PP .........SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS- CS TRM 259 .........aecehcgaklhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklikeElldaplyydlsqlasvlklsvp 336 ++ ++ + l+l+GP++ gpl++kef+++vl+ + k + r++++L+ +++E +dap+yyd+ +l+ ++k s p tr|A2FDH4|A2FDH4_TRIVA 277 rppnvefefTRDPYTNGPLKLGGPFYNGPLHNKEFIDKVLQLLPTMKYLGTIPRIEGVLNCCRNE-IDAPFYYDIGALTTIIKASCP 362 4667777789999********************************9888889*************.********************* PP --HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 337 plkkvvsalkelgfevsrtHakpkaikTnapleavvevlr 376 p + +++++ g+++s+tH+ p ikT+ap e ++++ r tr|A2FDH4|A2FDH4_TRIVA 363 PRALFYTQIEKQGYHCSLTHCDPGKIKTDAPPELLWDICR 402 *************************************975 PP >> tr|A2FGC5|A2FGC5_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_309480 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 187.4 0.4 1.4e-58 1.2e-54 15 165 .. 29 176 .. 13 179 .] 0.94 Alignments for each domain: == domain 1 score: 187.4 bits; conditional E-value: 1.4e-58 -TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--G CS TRM 15 skkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklnevee 101 ++++e+FYNPvq+fnRDl+++v ++++++ ++ +i++++a++asGlR+iR+a+Ev+gvk+++aND+ + Ave+ik+N+ +n+v++ tr|A2FGC5|A2FGC5_TRIVA 29 YTENEIFYNPVQQFNRDLTCAV---IQAYRNSIKDDISIFEAFAASGLRSIRYAKEVSGVKKIIANDLDQAAVEIIKRNILINKVSN 112 5689******************...77888888899*************************************************** PP GEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHH CS TRM 102 ivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspe 165 +v+++++dA+ + e++++f+v+DlDP+ ++apfl+ +v+++++++lL++T+TD +L+G +p+ tr|A2FGC5|A2FGC5_TRIVA 113 LVTASQGDAREKMLENENQFQVLDLDPYSTAAPFLEGVVKAATDGALLCITSTDGRTLCGVQPD 176 ************************************************************9887 PP >> tr|A2FGC4|A2FGC4_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_309470 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 85.0 0.0 1.8e-27 1.5e-23 259 376 .] 37 153 .. 11 153 .. 0.91 Alignments for each domain: == domain 1 score: 85.0 bits; conditional E-value: 1.8e-27 SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH CS TRM 259 aecehcgaklhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsal 345 ++ ++ g l+++GP++ gpl++kef+++vle + k ++ ++ r++++L+ +++E ++a +yyd+ +l+s++k s pp + +++ tr|A2FGC4|A2FGC4_TRIVA 37 TKDPYTGGALKIGGPFYNGPLHSKEFIDQVLELLPKMHNLQTIPRIEGVLNCCRNE-IEALFYYDIGALTSIIKASCPPRALIYTQI 122 778899999****************************9******************.****************************** PP ....TT-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 346 kelgfevsrtHakpkaikTnapleavvevlr 376 ++ +f+vs+tH+ ikT+ap e ++++ r tr|A2FGC4|A2FGC4_TRIVA 123 ERQNFHVSLTHCDAGKIKTDAPSELIWDICR 153 ****************************975 PP >> sp|A2E5K9|TRM5_TRIVA tRNA (guanine(37)-N1)-methyltransferase OS=Trichomonas vaginalis OX=5722 GN=TVAG_038970 PE=3 SV # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 34.8 0.0 3.4e-12 2.9e-08 18 144 .. 184 299 .. 156 303 .. 0.81 Alignments for each domain: == domain 1 score: 34.8 bits; conditional E-value: 3.4e-12 --SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE CS TRM 18 eeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvin 106 ++V++N + + D +v++ +++ v da+ + G a+R a+ +v aND+++ + +keN k+n v+e v+ sp|A2E5K9|TRM5_TRIVA 184 SKVYWNSRLQYEHDSVVAT-----FKEN-----SLVCDAMCGIGPFAVRAAKK---GCRVRANDLNPDSYYWLKENCKINGVSENVECF 259 5688888888888866666.....2333.....3699**************98...67899**************************** PP -S-HHHHHC..--SS-EEEEEE--SS--GGGHHHHHHHEE CS TRM 107 eddAnvlla..ekkkrfdviDlDPfGspapfldsavqsvk 144 ++dA+ +++ + d+i + G++++fld+ + +k sp|A2E5K9|TRM5_TRIVA 260 NMDAREFIRkqFDNGGCDYIVMNLPGTAVEFLDAIGEGAK 299 ********9544566799*****99********9766555 PP >> tr|A2D9B5|A2D9B5_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_183780 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.4 0.1 5.3e-06 0.044 248 313 .. 47 115 .. 21 126 .. 0.74 Alignments for each domain: == domain 1 score: 14.4 bits; conditional E-value: 5.3e-06 E.....EEE--GGG--SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHH CS TRM 248 e.....ekvkglakfkaecehcgaklhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklike 313 + e + ++k+++e +h + + GP w ++ +++ef+++v+ ++e+ +e++++r+++ + +k+ tr|A2D9B5|A2D9B5_TRIVA 47 KqleyfEPNGAMQKTNQEFCHSIQICSKPGPKWQQNYNNNEFLDRVIGAMEE--QENNARRIREAAQNLKN 115 23333333444556567777888889999**********************9..68888899887766665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 2582 (0.0514445); expected 1003.8 (0.02) Passed bias filter: 1318 (0.0262602); expected 1003.8 (0.02) Passed Vit filter: 89 (0.00177326); expected 50.2 (0.001) Passed Fwd filter: 6 (0.000119546); expected 0.5 (1e-05) Initial search space (Z): 50190 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.46u 0.02s 00:00:00.48 Elapsed: 00:00:00.14 # Mc/sec: 45805.39 // [ok]