# 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/Petromyzon_marinus.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.9e-25 87.7 0.0 6.2e-25 87.4 0.0 1.0 1 tr|S4RSM5|S4RSM5_PETMA rRNA adenine N(6)-methyltransferase O 3.3e-18 65.3 0.0 4.1e-18 65.0 0.0 1.0 1 tr|S4RUG9|S4RUG9_PETMA rRNA adenine N(6)-methyltransferase O 3.7e-11 42.2 0.0 7.9e-09 34.6 0.0 2.2 2 tr|S4R4I3|S4R4I3_PETMA Phosphoethanolamine methyltransferase 8.9e-08 31.1 0.0 7.1e-07 28.2 0.0 2.1 2 tr|S4RMG0|S4RMG0_PETMA Phosphoethanolamine methyltransferase 0.00073 18.3 0.0 0.0012 17.6 0.0 1.2 1 tr|S4RWE3|S4RWE3_PETMA Acetylserotonin O-methyltransferase O ------ inclusion threshold ------ 0.083 11.5 0.0 0.12 11.0 0.0 1.2 1 tr|S4R939|S4R939_PETMA Methyltransferase like 9 OS=Petromyzo 0.14 10.8 0.0 0.2 10.3 0.0 1.1 1 tr|S4RA41|S4RA41_PETMA Acetylserotonin O-methyltransferase O Domain annotation for each sequence (and alignments): >> tr|S4RSM5|S4RSM5_PETMA rRNA adenine N(6)-methyltransferase OS=Petromyzon marinus OX=7757 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 87.4 0.0 3.2e-28 6.2e-25 12 263 .. 6 304 .. 2 306 .. 0.90 Alignments for each domain: == domain 1 score: 87.4 bits; conditional E-value: 3.2e-28 EE--HHHHHHHHHHH-GG..TTSEEEEES-TTTCCHHHHHCC-.SEEEEE-S-HHHHHHHHCHTT.TCCCEEEE-S-CCGHHHHHHH CS RrnaAD 12 flknkkvineIldklnle..ekltvleigpGkGilTeeLakra.kqvvvlElderLakllqkkls.kdeklkvvhqDflkfetpeev 94 f++++ v+ +I ++l+++ ++ v + +G+GilT+ L++ + vvvlE d+ ++ +lq+ + ++++++ h+D+ k +t + tr|S4RSM5|S4RSM5_PETMA 6 FITEPRVAHSITEALQINegDQAVVFDCSAGPGILTRALLNSGaNRVVVLERDKLFLPPLQELEEsAEGRMNLFHCDVFKLDTHDGL 92 999***************6667789**************66655******************5558999****************99 PP CT..................TTCCCCEEEEEEE--TTTHHHHHHHHHH.......GGGGCEEEEEEEEEHHHHHHHH--TTSTT-SH CS RrnaAD 95 te..................kdsihqeflvvanlPyeistqlvkqlle.......esrfglvkmllwlqkefarrllarpgskkrsr 156 + +++ +++ v+ lP + +l+ rfg+v+ l + + +++ ++pg k++s+ tr|S4RSM5|S4RSM5_PETMA 93 LRppfmsskqlqvqlglkerSWNDEIPIRLVGILPLRNPRTILWRLIYaafqgisVFRFGRVELNLFISERELQKIRSSPGCKNYSS 179 99****************99999999**********999999999998899999999****************************** PP HHHHHHHHEEEEEEEEC-GG................CECS--SSSEEEEEEEE-SS..-S.--S-HHHHHHHHHHHHCTTTSBCCCH CS RrnaAD 157 lsvlreaftdvklvakveks................ifsPppkvdsalvelerkde..dl.pqkdlkkyesvvrkllnrkrktlsts 224 ls+l +++ d+++v++ + + + +++ l++l+++ e ++ + +++++ ++++++l r+rk+l+ tr|S4RSM5|S4RSM5_PETMA 180 LSILWQIMSDITVVYQESLTsvvtcsnkalktdffsNHVEA-GEKLCLMQLNPRREitSVlNYHSMDRFMAMINHCLSRRRKKLKDI 265 **************9987779***********977543344.7899********98223455667899*****************99 PP CCCTSTCHHHHHCTT..T--TTTCCCCSGHHHHHHHHHHH CS RrnaAD 225 lksllpggelqalsst.gindnalvkklsaeqtldifkel 263 l+ l + g +++ + g+ +a+ + +s+ ++l++f+++ tr|S4RSM5|S4RSM5_PETMA 266 LN-LWNPGSGELIMAQlGLSPEATSGTVSPRDYLRLFEII 304 88.55555556666777********************998 PP >> tr|S4RUG9|S4RUG9_PETMA rRNA adenine N(6)-methyltransferase OS=Petromyzon marinus OX=7757 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 65.0 0.0 2.1e-21 4.1e-18 104 242 .. 5 151 .. 2 169 .. 0.90 Alignments for each domain: == domain 1 score: 65.0 bits; conditional E-value: 2.1e-21 EEEEEE--TTTHHHHHHHHHH.......GGGGCEEEEEEEEEHHHHHHHH.--TTSTT-SHHHHHHHHHEEEEEEEEC-GGCECS-- CS RrnaAD 104 flvvanlPyeistqlvkqlle.......esrfglvkmllwlqkefarrll.arpgskkrsrlsvlreaftdvklvakveksifsPpp 182 + +++nlP +st+lv +le g+ + +l +qke+a+rll a+ s++r+rlsv+++ + +vkl +++ s+f P p tr|S4RUG9|S4RUG9_PETMA 5 IHIIGNLPFSVSTALVIDWLEcmanrsgPFSYGRTRLTLTFQKEVAERLLvASVHSSQRCRLSVMAQSLSSVKLRFIIPGSAFIPVP 91 689********************999884445689*************7548999******************************** PP SSSEEEEEEEE-SS-S--S-HHHHHHHHHHHHCTTTSBCCCHCCCTSTCHHHHHCTT.T- CS RrnaAD 183 kvdsalvelerkdedlpqkdlkkyesvvrkllnrkrktlstslksllpggelqalsstgi 242 kvd +v le+ ++ q+ + e+vvr l+ +rk + l+p++e+ +ls++ + tr|S4RUG9|S4RUG9_PETMA 92 KVDVGVVTLEPLVTPRIQQPFPLVEKVVRSLFQFRRKYCLRGVELLFPQAERLSLSQRML 151 ***********988657999*******************************999987643 PP >> tr|S4R4I3|S4R4I3_PETMA Phosphoethanolamine methyltransferase OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 34.6 0.0 4.1e-12 7.9e-09 30 95 .. 50 115 .. 36 138 .. 0.89 2 ! 5.3 0.0 0.0034 6.6 19 104 .. 272 367 .. 266 371 .. 0.73 Alignments for each domain: == domain 1 score: 34.6 bits; conditional E-value: 4.1e-12 TTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S-CCGHHHHHHHC CS RrnaAD 30 ekltvleigpGkGilTeeLakrakqvvvlElderLakllqkklskdeklkvvhqDflkfetpeevt 95 ++ +vle+g+G G +T +Lak+a++v+++ + ++ ++ s++ ++k +++D+ + ++pe+ tr|S4R4I3|S4R4I3_PETMA 50 DDAKVLELGAGIGRFTGHLAKQARHVTAVDFMQTFLDKNRELNSDHVNIKFLQADVTQLQFPENSF 115 6789*************************************888******************9754 PP == domain 2 score: 5.3 bits; conditional E-value: 0.0034 HHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-.SEEEEE-S-HHHHHHHHCHTTTCCCEEE..E-S-CCGHHHHHHHCT...... CS RrnaAD 19 ineIldklnleekltvleigpGkGilTeeLakra.kqvvvlElderLakllqkklskdeklkv..vhqDflkfetpeevte...... 96 +e ++ l+l+ ++vl++g G G +ak + v + l ++++++ ++ ++d+ l v D+ k ++pee + tr|S4R4I3|S4R4I3_PETMA 272 TKEFVSLLDLQAGQKVLDVGCGIGGGNFYMAKAFkVDVLGMDLSSNMVEIAMERATEDSTLLVqfEISDATKRDFPEESFDvvysrd 358 67788899****************999999999944566678888888888888887666655424568888888876666666655 PP .TTCCCCEE CS RrnaAD 97 .kdsihqef 104 ++i ++f tr|S4R4I3|S4R4I3_PETMA 359 tILHIADKF 367 555555555 PP >> tr|S4RMG0|S4RMG0_PETMA Phosphoethanolamine methyltransferase OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 28.2 0.0 3.7e-10 7.1e-07 30 95 .. 47 112 .. 36 134 .. 0.86 2 ? -0.0 0.0 0.15 2.9e+02 21 42 .. 272 293 .. 266 295 .. 0.83 Alignments for each domain: == domain 1 score: 28.2 bits; conditional E-value: 3.7e-10 TTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S-CCGHHHHHHHC CS RrnaAD 30 ekltvleigpGkGilTeeLakrakqvvvlElderLakllqkklskdeklkvvhqDflkfetpeevt 95 + + +le+g+G G +T +Lak+a++v+++ e++++ ++ ++l+ +++D+ + ++p + tr|S4RMG0|S4RMG0_PETMA 47 DGTDILELGAGIGRFTGHLAKIANHVTAVDFMENFIAKNKEDNGHRGNLTFLQADVTQLDFPANSF 112 5678*************************************8889*****************9765 PP == domain 2 score: -0.0 bits; conditional E-value: 0.15 HHHHHH-GGTTSEEEEES-TTT CS RrnaAD 21 eIldklnleekltvleigpGkG 42 e l l l++ ++vl++g G G tr|S4RMG0|S4RMG0_PETMA 272 EFLAMLGLRQGQKVLDVGCGIG 293 5667788999**********98 PP >> tr|S4RWE3|S4RWE3_PETMA Acetylserotonin O-methyltransferase OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.6 0.0 6.2e-07 0.0012 22 92 .. 177 248 .. 175 275 .. 0.87 Alignments for each domain: == domain 1 score: 17.6 bits; conditional E-value: 6.2e-07 HHHHH-GGTTSEEEEES-TTTCCHHHHHCC-..SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S-CCGHHHHH CS RrnaAD 22 IldklnleekltvleigpGkGilTeeLakra..kqvvvlElderLakllqkklskdeklkvvhqDflkfetpe 92 Il++ +l+ + v ++g G+G+l +++a + v+++ l ++ +++ qk++ +++++k+ Df k ++pe tr|S4RWE3|S4RWE3_PETMA 177 ILNAFDLSSFKDVCDVGGGTGVLAQQFARMYpdARVTLFDL-PEVLEVAQKHFLSSGNVKLHAGDFFKGKIPE 248 8999999999*********************9666666665.666799*******************998887 PP >> tr|S4R939|S4R939_PETMA Methyltransferase like 9 OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.0 0.0 6.1e-05 0.12 32 71 .. 148 187 .. 130 208 .. 0.90 Alignments for each domain: == domain 1 score: 11.0 bits; conditional E-value: 6.1e-05 SEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCH CS RrnaAD 32 ltvleigpGkGilTeeLakrakqvvvlElderLakllqkk 71 + +l++g+G G T+++ +kqv v E + + +l+k+ tr|S4R939|S4R939_PETMA 148 DSLLDLGAGDGEVTAIMGAHFKQVFVTEVSPTMQWRLRKR 187 789***************************9999999984 PP >> tr|S4RA41|S4RA41_PETMA Acetylserotonin O-methyltransferase OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.3 0.0 0.0001 0.2 22 122 .. 181 273 .. 178 280 .. 0.81 Alignments for each domain: == domain 1 score: 10.3 bits; conditional E-value: 0.0001 HHHHH-GGTTSEEEEES-TTTCCHHHHHCC-..SEEEEE-S-HHHHHHHHCHTTTCCCEEEE-S-CCGHHHHHHHCTTTCCCCEEEE CS RrnaAD 22 IldklnleekltvleigpGkGilTeeLakra..kqvvvlElderLakllqkklskdeklkvvhqDflkfetpeevtekdsihqeflv 106 +l + +l++ + v+++g +G+l +La++ v+v+ + + +++ +k++ ++ +++ Df k e+p+ + + tr|S4RA41|S4RA41_PETMA 181 VLTAFDLSPFKDVYDLGGCTGALALQLAEIHpdSRVTVF-DLPPVVEVARKHFLASGTVQLQGGDFFKDEIPQ---------GDLYI 257 5778899999**********************5555555.557789**********************99998.........45566 PP EEE--TTTHHHHHHHH CS RrnaAD 107 vanlPyeistqlvkql 122 a + + s++ +ql tr|S4RA41|S4RA41_PETMA 258 LARILHDWSKEQCMQL 273 6666666655555555 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (265 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 344 (0.0254081); expected 270.8 (0.02) Passed bias filter: 273 (0.020164); expected 270.8 (0.02) Passed Vit filter: 32 (0.00236354); expected 13.5 (0.001) Passed Fwd filter: 7 (0.000517025); expected 0.1 (1e-05) Initial search space (Z): 13539 [actual number of targets] Domain search space (domZ): 7 [number of targets reported over threshold] # CPU time: 0.12u 0.01s 00:00:00.13 Elapsed: 00:00:00.04 # Mc/sec: 33448.91 // [ok]