# 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/Methyltransf_4.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_4 [M=173] Accession: PF02390.16 Description: Putative 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-43 149.3 0.0 2.5e-43 149.0 0.0 1.1 1 tr|A2ESN8|A2ESN8_TRIVA tRNA (guanine-N(7)-)-methyltransferas 9.3e-42 143.9 0.0 1.1e-41 143.6 0.0 1.1 1 tr|A2DE67|A2DE67_TRIVA tRNA (guanine-N(7)-)-methyltransferas 2.4e-06 28.6 0.2 4.7e-06 27.6 0.2 1.4 1 tr|A2ELU0|A2ELU0_TRIVA Uncharacterized protein OS=Trichomona ------ inclusion threshold ------ 0.081 13.8 0.9 0.45 11.4 0.2 2.5 2 tr|A2DNW0|A2DNW0_TRIVA Uncharacterized protein OS=Trichomona 0.14 13.1 0.0 7.5 7.4 0.0 2.3 2 tr|A2FK19|A2FK19_TRIVA Methyltransferase, putative OS=Tricho 0.26 12.2 0.1 0.53 11.2 0.1 1.5 1 tr|A2FNS0|A2FNS0_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2ESN8|A2ESN8_TRIVA tRNA (guanine-N(7)-)-methyltransferase OS=Trichomonas vaginalis OX=5722 GN=TVAG_069970 PE=3 S # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 149.0 0.0 2.9e-47 2.5e-43 3 172 .. 61 232 .. 59 233 .. 0.95 Alignments for each domain: == domain 1 score: 149.0 bits; conditional E-value: 2.9e-47 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHC....TTSSEEEEES-CCCHHHHCS-TTSEEEEEEES----- CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkieal....elkNlrilrgdalkvllelledgslekifilFPDPwf 85 p++l+iGcG+G+++ al++++P+++++++Eir++v+ +v kki+a+ + N+++ + +++l ++l+ +++ekifilFPDP+f tr|A2ESN8|A2ESN8_TRIVA 61 PTMLDIGCGYGGMTFALSTEFPQNLIFALEIRSQVTTYVEKKIQAYhaqgLALNIQVQWANTMRTLMRYLRAHTIEKIFILFPDPHF 147 789*******************************************55544559********************************* PP SGGGGGGSTTSHHHHHHHHHHEEEEEEEEEEESHHHHHHHHHHHHHCSTTEEESSCCCSCCTCHCHHCCCSSTTCCCCCCTTSS-EE CS Methyltransf_4 86 KkrhhKrRllqpellkelarvLkkgGvlelaTDveeyaeemlkvleeeplferlskeeeleeeplsvlriatkyeqkvqskgkeiye 172 Kk++ K+R+++++ll+e+a+++k g+ ++l+TDv +y+e ++ le++plf+r+++ ++++++ + +t++ +kv+++g++++ tr|A2ESN8|A2ESN8_TRIVA 148 KKKKAKWRIISEQLLDEYAYIMKGGARFYLVTDVRQYFEYAVPILEAHPLFKRVTDI--DSDKCVNLAKTSTEESNKVTENGGQKFA 232 *****************************************************9987..567899*****************99885 PP >> tr|A2DE67|A2DE67_TRIVA tRNA (guanine-N(7)-)-methyltransferase OS=Trichomonas vaginalis OX=5722 GN=TVAG_166600 PE=3 S # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 143.6 0.0 1.3e-45 1.1e-41 3 172 .. 61 232 .. 59 233 .. 0.95 Alignments for each domain: == domain 1 score: 143.6 bits; conditional E-value: 1.3e-45 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTS....SEEEEES-CCCHHHHCS-TTSEEEEEEES----- CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelk....NlrilrgdalkvllelledgslekifilFPDPwf 85 p++l+iGcG G++ al++++P++++l++Eir++v +v kki+a +++ N+++ ++ +++l ++++ +++ekifilFPDP+f tr|A2DE67|A2DE67_TRIVA 61 PTILDIGCGFGGMSFALSPEFPNNLILALEIRSQVATYVEKKIQAFHAQglslNIQVQWTNTMRTLMRYIRAHTIEKIFILFPDPHF 147 789*****************************************995444445********************************** PP SGGGGGGSTTSHHHHHHHHHHEEEEEEEEEEESHHHHHHHHHHHHHCSTTEEESSCCCSCCTCHCHHCCCSSTTCCCCCCTTSS-EE CS Methyltransf_4 86 KkrhhKrRllqpellkelarvLkkgGvlelaTDveeyaeemlkvleeeplferlskeeeleeeplsvlriatkyeqkvqskgkeiye 172 Kk++ K+R+++++ll+e+a ++k g+ ++l+TDv +y+e ++ le++plf+r+++ ++++++ + +t++ +kv+++g++++ tr|A2DE67|A2DE67_TRIVA 148 KKKKSKWRIISQQLLDEYAFIMKGGARMYLVTDVRQYFEYAVPLLEAHPLFKRITDL--DSDKCVELSKTSTEESNKVSENGGQKFA 232 *****************************************************9987..56789******************99885 PP >> tr|A2ELU0|A2ELU0_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_151750 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 27.6 0.2 5.7e-10 4.7e-06 3 70 .. 377 442 .. 375 454 .. 0.90 Alignments for each domain: == domain 1 score: 27.6 bits; conditional E-value: 5.7e-10 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTSSEEEEES-CCCHHHHCS-T CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelkNlrilrgdalkvllelled 70 v+++++cG G + + +Akk k+++Gi+i+++ ++ + k++ ++++kN + + a++vl+ +l + tr|A2ELU0|A2ELU0_TRIVA 377 TVLIDVCCGTGVIGLSMAKKA--KKVIGIDIEKDAIEDAKKNATKNKIKNAEFIAAPAQDVLPGILDQ 442 689********9999999985..889************************************988765 PP >> tr|A2DNW0|A2DNW0_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_225810 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.4 0.2 5.4e-05 0.45 22 78 .. 319 380 .. 310 381 .. 0.88 2 ? -1.4 0.0 0.47 3.9e+03 52 81 .. 543 577 .. 535 578 .. 0.77 Alignments for each domain: == domain 1 score: 11.4 bits; conditional E-value: 5.4e-05 HSTTSEEEEEES-HHHHHHHHHHHHHCTTSSEEEEES-C.....CCHHHHCS-TTSEEEEEE CS Methyltransf_4 22 knPdknflGiEirvkvvakvlkkiealelkNlrilrgda.....lkvllelledgslekifi 78 + P++++l +E+ +a vl++i +++l+Nl+ ++ a ++ l +++ ++sle i++ tr|A2DNW0|A2DNW0_TRIVA 319 NIPNIQYLSLEKCDLDIADVLNSITKYRLNNLKAINLSAnnnidTNSLDSIILPESLESISV 380 56999*******************************99977777788888888888888765 PP == domain 2 score: -1.4 bits; conditional E-value: 0.47 SEEEEES-C.....CCHHHHCS-TTSEEEEEEES- CS Methyltransf_4 52 Nlrilrgda.....lkvllelledgslekifilFP 81 Nl++l d l v++++le++ + +f+ FP tr|A2DNW0|A2DNW0_TRIVA 543 NLQVLVADGlcpdsLPVFRDFLETAANSHVFLSFP 577 55555555444444999*****************9 PP >> tr|A2FK19|A2FK19_TRIVA Methyltransferase, putative OS=Trichomonas vaginalis OX=5722 GN=TVAG_346300 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 7.4 0.0 0.0009 7.5 4 48 .. 67 112 .. 64 138 .. 0.82 2 ? 3.2 0.0 0.018 1.5e+02 97 115 .. 158 176 .. 156 183 .. 0.87 Alignments for each domain: == domain 1 score: 7.4 bits; conditional E-value: 0.0009 EEEEET-TTSHHH.HHHHHHSTTSEEEEEES-HHHHHHHHHHHHHC CS Methyltransf_4 4 vvleiGcGkGefl.valAkknPdknflGiEirvkvvakvlkkieal 48 + +++G+G G+ v A + + n++Gi++ ++ ++ ++++i+ + tr|A2FK19|A2FK19_TRIVA 67 TSIDLGSGTGRDCfVISALAGKTGNVIGIDMTQEQIDVANEAIAYH 112 579*******86514455666678***************9999866 PP == domain 2 score: 3.2 bits; conditional E-value: 0.018 HHHHHHHHHHEEEEEEEEE CS Methyltransf_4 97 pellkelarvLkkgGvlel 115 +e+ ke+ rvLk+gG+l++ tr|A2FK19|A2FK19_TRIVA 158 KEVFKEIYRVLKEGGELHF 176 67889*************8 PP >> tr|A2FNS0|A2FNS0_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_137180 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.2 0.1 6.3e-05 0.53 2 56 .. 349 401 .. 348 411 .. 0.90 Alignments for each domain: == domain 1 score: 11.2 bits; conditional E-value: 6.3e-05 -EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTSSEEEE CS Methyltransf_4 2 apvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelkNlril 56 + v+++ +cG G + l+k k ++Gi+ +++ +a+ +++e++++kN + + tr|A2FNS0|A2FNS0_TRIVA 349 NTVLVDFCCGTGVNCLSLSKS--VKRCIGIDSSESNIATSQRNAEENQIKNAYFV 401 6789********999999876..699*************************8765 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (173 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 2677 (0.0533373); expected 1003.8 (0.02) Passed bias filter: 1173 (0.0233712); expected 1003.8 (0.02) Passed Vit filter: 101 (0.00201235); 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.28u 0.01s 00:00:00.29 Elapsed: 00:00:00.11 # Mc/sec: 26823.18 // [ok]