# 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/RrnaAD.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: RrnaAD [M=265] Accession: PF00398.19 Description: Ribosomal RNA adenine dimethylase 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.2e-58 198.0 0.0 5.2e-58 197.7 0.0 1.1 1 tr|A2EVN6|A2EVN6_TRIVA rRNA adenine N(6)-methyltransferase O 0.00017 22.2 0.0 0.00022 21.9 0.0 1.1 1 tr|A2FEZ1|A2FEZ1_TRIVA S-adenosylmethionine-dependent methyl 0.0058 17.2 0.0 0.0066 17.0 0.0 1.2 1 tr|A2F2A3|A2F2A3_TRIVA Menaquinone biosynthesis methyltransf 0.0086 16.6 0.1 0.032 14.8 0.0 1.8 2 tr|A2FAI4|A2FAI4_TRIVA Uncharacterized protein OS=Trichomona ------ inclusion threshold ------ 0.018 15.6 0.0 0.025 15.1 0.0 1.2 1 tr|A2DNX4|A2DNX4_TRIVA Arginine methyltransferase, putative 0.042 14.4 0.0 0.067 13.7 0.0 1.3 1 tr|A2EXH2|A2EXH2_TRIVA Uncharacterized protein OS=Trichomona 0.047 14.2 0.0 0.064 13.8 0.0 1.2 1 tr|A2FCQ2|A2FCQ2_TRIVA Arginine methyltransferase, putative 0.077 13.5 0.0 0.12 12.9 0.0 1.2 1 tr|A2H3A0|A2H3A0_TRIVA Uncharacterized protein OS=Trichomona 0.12 12.9 0.0 0.18 12.3 0.0 1.2 1 tr|A2G2M2|A2G2M2_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2EVN6|A2EVN6_TRIVA rRNA adenine N(6)-methyltransferase OS=Trichomonas vaginalis OX=5722 GN=TVAG_299080 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 197.7 0.0 9.4e-62 5.2e-58 3 254 .. 16 249 .. 14 262 .. 0.89 Alignments for each domain: == domain 1 score: 197.7 bits; conditional E-value: 9.4e-62 CCXXGCG--EE--HHHHHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-.SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S-CCGH CS RrnaAD 3 klrksygqnflknkkvineIldklnleekltvleigpGkGilTeeLakra.kqvvvlElderLakllqkklskdeklkvvhqDflkf 88 k+ s+gqn+l++k v+++I+d+ ++++ +++leigpG G++Tee++ r +v ++E+d+r++ +l+kk+ +++l ++++D+l tr|A2EVN6|A2EVN6_TRIVA 16 KMDHSLGQNILRSKVVVKNIVDAGEPRPGDKILEIGPGNGNMTEEMLSREgIEVIAIEKDQRMCVELKKKFPRHPNLRIINADVLSV 102 67789********************************************99************************************ PP HHHHHHCTTTCCCCEEEEEEE--TTTHHHHHHHHHHGGGGCEEEEEEEEEHHHHHHHH--TTSTT-SHHHHHHHHHEEEEEEEEC-G CS RrnaAD 89 etpeevtekdsihqeflvvanlPyeistqlvkqlleesrfglvkmllwlqkefarrllarpgskkrsrlsvlreaftdvklvakvek 175 ++pe+ l+++n+Py+is+++v++ll ++f ++ +l++qkef +r++arpg+ ++rl ++++++++vklv++v + tr|A2EVN6|A2EVN6_TRIVA 103 DLPEFD----------LCISNIPYNISSAIVFKLLARPTF--RRTVLMVQKEFGERIVARPGKDGWGRLAINTQLYASVKLVMNVSR 177 ****77..........9******************98888..9******************************************** PP GCECS--SSSEEEEEEEE-SS-S--S-HHHHHHHHHHHHCTTTSBCCCHCCCTSTCHHHHHCTT.T--TTTCCCCSGHH CS RrnaAD 176 sifsPppkvdsalvelerkdedlpqkdlkkyesvvrkllnrkrktlstslksllpggelqalsstgindnalvkklsae 254 ++f Pppkvds +v +++++ + p+ d ++++ +++ +++r ++tl++ +++ ++ +al ++ n ++ + ++ + tr|A2EVN6|A2EVN6_TRIVA 178 KNFVPPPKVDSIVVSIQPREVQ-PDIDFNEWDGLIKLCFTRPNRTLQATFHK----KKIRALLNE--NREKYLSSINQQ 249 *********************9.8***********************98653....333333332..333333333333 PP >> tr|A2FEZ1|A2FEZ1_TRIVA S-adenosylmethionine-dependent methyltransferase, putative OS=Trichomonas vaginalis OX=5722 G # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 21.9 0.0 3.9e-08 0.00022 19 93 .. 29 105 .. 22 133 .. 0.88 Alignments for each domain: == domain 1 score: 21.9 bits; conditional E-value: 3.9e-08 HHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTT..TCCCEEEE-S-CCGHHHHHH CS RrnaAD 19 ineIldklnleekltvleigpGkGilTeeLakrakqvvvlElderLakllqkkls..kdeklkvvhqDflkfetpee 93 ++e+ +++n++++ +l+ g G+G+ L + ak l + +++ ++k+l +++++++++q++ ++ + + tr|A2FEZ1|A2FEZ1_TRIVA 29 VAEVRKQANITPESRILDFGCGTGLNGIYLINDAKTIGFLDPSSGMIEQVKKELAeiNKTNVEIYNQEISQTTMEPF 105 667889**********************************************999999***********99877664 PP >> tr|A2F2A3|A2F2A3_TRIVA Menaquinone biosynthesis methyltransferase, putative OS=Trichomonas vaginalis OX=5722 GN=TVAG # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.0 0.0 1.2e-06 0.0066 32 96 .. 60 125 .. 46 190 .. 0.87 Alignments for each domain: == domain 1 score: 17.0 bits; conditional E-value: 1.2e-06 SEEEEES-TTTCCHHHH.HCC-SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S-CCGHHHHHHHCT CS RrnaAD 32 ltvleigpGkGilTeeL.akrakqvvvlElderLakllqkklskdeklkvvhqDflkfetpeevte 96 t l +g G ++T eL ++ +++vv + e+ + ++kk ++++kl++ D+ k+++p++ + tr|A2F2A3|A2F2A3_TRIVA 60 GTALNLGCGNSNMTSELlLNGFDKVVGIDFSEVVIGQMRKKYQLEQKLEWETGDITKMKFPNNHFD 125 57799**********99899******************************************7655 PP >> tr|A2FAI4|A2FAI4_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_415010 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 14.8 0.0 5.7e-06 0.032 19 83 .. 25 91 .. 14 112 .. 0.88 2 ? -0.4 0.0 0.24 1.4e+03 160 196 .. 109 143 .. 98 167 .. 0.61 Alignments for each domain: == domain 1 score: 14.8 bits; conditional E-value: 5.7e-06 HHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTT..TCCCEEEE-S CS RrnaAD 19 ineIldklnleekltvleigpGkGilTeeLakrakqvvvlElderLakllqkkls..kdeklkvvhq 83 ++ I +++nl++++tvl+ g kG+l +L + ak+ l d++ + + ++ l+ +++ +++++ tr|A2FAI4|A2FAI4_TRIVA 25 VDIIKSQSNLNKSMTVLDFGCEKGFLGLKLIDDAKKIAFLDPDPNAISIAEESLKmmGKTNYEIINK 91 5667889*************************************99999999888888888888876 PP == domain 2 score: -0.4 bits; conditional E-value: 0.24 HHHHHEEEEEEEEC-GGCECS--SSSEEEEEEEE-SS CS RrnaAD 160 lreaftdvklvakveksifsPppkvdsalvelerkde 196 ++++ dvk ++ ks+ P+ +++++e+ k++ tr|A2FAI4|A2FAI4_TRIVA 109 SLHYVPDVKSALTKIKSLLNPD--GKLVIIEIYHKHF 143 4444445544444444444444..5556666655555 PP >> tr|A2DNX4|A2DNX4_TRIVA Arginine methyltransferase, putative OS=Trichomonas vaginalis OX=5722 GN=TVAG_225950 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 15.1 0.0 4.4e-06 0.025 7 93 .. 20 113 .. 14 142 .. 0.75 Alignments for each domain: == domain 1 score: 15.1 bits; conditional E-value: 4.4e-06 GCG--EE--HHHHHHHHHHH-GG....TTSEEEEES-TTTCCHHHHHCC-.SEEEEE-S-HHHHHHHHCHTT....TCCCEEEE-S- CS RrnaAD 7 sygqnflknkkvineIldklnle....ekltvleigpGkGilTeeLakra.kqvvvlElderLakllqkkls....kdeklkvvhqD 84 s+ +l+++ v + +++ + +++ vleig+G+Gi++ ak k+v + E ++L+++++kk ++ k++++ q+ tr|A2DNX4|A2DNX4_TRIVA 20 SFHDIILRDTAVFSQFEKAIRYNstlfKDKIVLEIGSGTGIFSMMAAKSGaKHVYAWE--PSLLNIYSKKTIednnLQDKITILSQN 104 5555566666666655554444445559999***************999857777776..899999998876577788889999998 PP CCGHHHHHH CS RrnaAD 85 flkfetpee 93 + ++ e+ tr|A2DNX4|A2DNX4_TRIVA 105 LEEIKLEEK 113 877776665 PP >> tr|A2EXH2|A2EXH2_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_253130 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.7 0.0 1.2e-05 0.067 32 89 .. 77 135 .. 56 163 .. 0.84 Alignments for each domain: == domain 1 score: 13.7 bits; conditional E-value: 1.2e-05 SEEEEES-TTTCCHHHH.HCC-SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S-CCGHH CS RrnaAD 32 ltvleigpGkGilTeeL.akrakqvvvlElderLakllqkklskdeklkvvhqDflkfe 89 ++ +ig+G G +e + k++k+ ++E ++++ + ++kl+++ lk + + + +++ tr|A2EXH2|A2EXH2_TRIVA 77 QKCADIGAGIGRVSELIlSKYFKEIDLVEPVQKFVDVAKEKLKNKVILKTYTCGAQDWK 135 6778**********96659***********************98888888888887776 PP >> tr|A2FCQ2|A2FCQ2_TRIVA Arginine methyltransferase, putative OS=Trichomonas vaginalis OX=5722 GN=TVAG_048280 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.8 0.0 1.1e-05 0.064 7 93 .. 20 113 .. 14 141 .. 0.75 Alignments for each domain: == domain 1 score: 13.8 bits; conditional E-value: 1.1e-05 GCG--EE--HHHHHHHHHHH-GG....TTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTT....TCCCEEEE-S-C CS RrnaAD 7 sygqnflknkkvineIldklnle....ekltvleigpGkGilTeeLakrakqvvvlElderLakllqkkls....kdeklkvvhqDf 85 s+ +l+++ v + +++ + +++ vleig+G+Gi++ ak ++v ++L+++++kk ++ k++++ q+ tr|A2FCQ2|A2FCQ2_TRIVA 20 SFHDIILRDTAVFSQFEKAIRYNstlfKDKIVLEIGSGTGIFSMMAAKSGAKHVYA-WKPSLLNIYSKKTIednnLQDKITILSQNL 105 5555566666666655554444445559999*****************99888755.567888888777655777888899999988 PP CGHHHHHH CS RrnaAD 86 lkfetpee 93 + ++ e+ tr|A2FCQ2|A2FCQ2_TRIVA 106 EEIKLEEK 113 77766665 PP >> tr|A2H3A0|A2H3A0_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_207000 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.9 0.0 2.2e-05 0.12 55 126 .. 108 186 .. 103 203 .. 0.89 Alignments for each domain: == domain 1 score: 12.9 bits; conditional E-value: 2.2e-05 EEEE-S-HHHHHHHHCHTTTCCCEEE.E-S-CCGHHHHHHHCT........TTCCCCEEEEEEE--TTTHHHHHHHHHHGG CS RrnaAD 55 vvvlElderLakllqkklskdeklkv.vhqDflkfetpeevte........kdsihqeflvvanlPyeistqlvkqllees 126 ++ +E++ +L+ +++ +k++ +k+ ++ ++k + + k +++ +vv ++P ++++q+++ +++++ tr|A2H3A0|A2H3A0_TRIVA 108 LTNIESRIALVDFYNL--KKANSIKIlFNKQIYKDRILDVTYAplkkildpKNPPNNGTIVVFHIPEQVTEQMIFSIFSKY 186 789***********99..89999**9889999999999977779999******************************9955 PP >> tr|A2G2M2|A2G2M2_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_467070 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.3 0.0 3.3e-05 0.18 55 126 .. 170 248 .. 165 265 .. 0.89 Alignments for each domain: == domain 1 score: 12.3 bits; conditional E-value: 3.3e-05 EEEE-S-HHHHHHHHCHTTTCCCEEE.E-S-CCGHHHHHHHCT........TTCCCCEEEEEEE--TTTHHHHHHHHHHGG CS RrnaAD 55 vvvlElderLakllqkklskdeklkv.vhqDflkfetpeevte........kdsihqeflvvanlPyeistqlvkqllees 126 ++ +E++ +L+ +++ +k++ +k+ ++ ++k + + k +++ +vv ++P ++++q+++ +++++ tr|A2G2M2|A2G2M2_TRIVA 170 LTNIESRIALVDFYNL--KKANSIKIlFNKQIYKDRILDVTYAplkkildpKNPPNNGTIVVFHIPEQVTEQMIFSIFSKY 248 789***********99..89999**9889999999999977779999******************************9955 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (265 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 1918 (0.0382148); expected 1003.8 (0.02) Passed bias filter: 719 (0.0143256); expected 1003.8 (0.02) Passed Vit filter: 48 (0.000956366); expected 50.2 (0.001) Passed Fwd filter: 9 (0.000179319); expected 0.5 (1e-05) Initial search space (Z): 50190 [actual number of targets] Domain search space (domZ): 9 [number of targets reported over threshold] # CPU time: 0.34u 0.01s 00:00:00.35 Elapsed: 00:00:00.13 # Mc/sec: 34766.37 // [ok]