# 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/Pan_troglodytes.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.4e-66 222.7 0.0 6.6e-66 222.4 0.0 1.1 1 tr|H2QQY8|H2QQY8_PANTR rRNA adenine N(6)-methyltransfer 2.3e-42 145.3 0.0 2.8e-42 145.0 0.0 1.1 1 tr|H2QTY1|H2QTY1_PANTR rRNA adenine N(6)-methyltransfer 9.2e-20 71.2 0.0 1.8e-19 70.2 0.0 1.4 1 tr|H2Q1G2|H2Q1G2_PANTR rRNA adenine N(6)-methyltransfer 0.00094 18.7 0.0 0.0021 17.6 0.0 1.5 1 tr|H2QQ96|H2QQ96_PANTR Protein arginine methyltransfera ------ inclusion threshold ------ 0.019 14.4 0.0 0.032 13.7 0.0 1.3 1 tr|H2QHP4|H2QHP4_PANTR tRNA methyltransferase 61B OS=Pa 0.085 12.3 0.0 0.13 11.6 0.0 1.2 1 tr|A0A2I3RRF6|A0A2I3RRF6_PANTR Ribosomal RNA adenine dimethylas Domain annotation for each sequence (and alignments): >> tr|H2QQY8|H2QQY8_PANTR rRNA adenine N(6)-methyltransferase OS=Pan troglodytes OX=9598 GN=DIMT1 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 222.4 0.0 1.7e-69 6.6e-66 4 237 .. 30 254 .. 27 274 .. 0.96 Alignments for each domain: == domain 1 score: 222.4 bits; conditional E-value: 1.7e-69 CXXGCG--EE--HHHHHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTT...TCCCEEEE-S-CCG CS RrnaAD 4 lrksygqnflknkkvineIldklnleekltvleigpGkGilTeeLakrakqvvvlElderLakllqkkls...kdeklkvvhqDflk 87 +++ gq++lkn+ +in+I+dk+ l++++ vle+gpG+G++T++L++ ak+vv+ Eld+rL+++l k+++ ++kl+v+ D+lk tr|H2QQY8|H2QQY8_PANTR 30 FNTGIGQHILKNPLIINSIIDKAALRPTDVVLEVGPGTGNMTVKLLEKAKKVVACELDPRLVAELHKRVQgtpVASKLQVLVGDVLK 116 78999*****************************************************************9999************* PP HHHHHHHCTTTCCCCEEEEEEE--TTTHHHHHHHHHHGGGGCEEEEEEEEEHHHHHHHH--TTSTT-SHHHHHHHHHEEEEEEEEC- CS RrnaAD 88 fetpeevtekdsihqeflvvanlPyeistqlvkqlleesrfglvkmllwlqkefarrllarpgskkrsrlsvlreaftdvklvakve 174 +++p + +vanlPy+is+++v++ll ++r+ ++ ++l++q+efa rl+a+pg+k+++rls+++++++ v +++kv+ tr|H2QQY8|H2QQY8_PANTR 117 TDLPFFD----------TCVANLPYQISSPFVFKLL-LHRPFFRCAILMFQREFALRLVAKPGDKLYCRLSINTQLLARVDHLMKVG 192 *****88..........6******************.788889******************************************** PP GGCECS--SSSEEEEEEEE-SS-S--S-HHHHHHHHHHHHCTTTSBCCCHCCCTSTCHHHHHC CS RrnaAD 175 ksifsPppkvdsalvelerkdedlpqkdlkkyesvvrkllnrkrktlstslksllpggelqal 237 k++f+Pppkv+s++v++e+k+++ p + ++++ +vr+ ++rk+ktls+++ks ++ l+++ tr|H2QQY8|H2QQY8_PANTR 193 KNNFRPPPKVESSVVRIEPKNPP-PPINFQEWDGLVRITFVRKNKTLSAAFKSSAVQQLLEKN 254 **********************8.89**************************98888877766 PP >> tr|H2QTY1|H2QTY1_PANTR rRNA adenine N(6)-methyltransferase OS=Pan troglodytes OX=9598 GN=TFB1M PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 145.0 0.0 7.3e-46 2.8e-42 3 244 .. 28 282 .. 26 304 .. 0.89 Alignments for each domain: == domain 1 score: 145.0 bits; conditional E-value: 7.3e-46 CCXXGCG--EE--HHHHHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-.SEEEEE-S-HHHHHHHHCHTT..TCCCEEEE-S-CC CS RrnaAD 3 klrksygqnflknkkvineIldklnleekltvleigpGkGilTeeLakra.kqvvvlElderLakllqkkls..kdeklkvvhqDfl 86 ++ k+++qnfl + + ++I++k+ + v+e+gpG+G T+++++ ++ v+E+d+r++ lq ls +kl +vh D+l tr|H2QTY1|H2QTY1_PANTR 28 QAAKQLSQNFLLDLRLTDKIVRKAGNLTNAYVYEVGPGPGGITRSILNADvAELLVVEKDTRFIPGLQM-LSdaAPGKLRIVHGDVL 113 56789*******************99999***************9966655******************.5558889********** PP GHHHHHHHCT.....TTCCCCEEEEEEE--TTTHHHHHHHHHH.......GGGGCEEEEEEEEEHHHHHHHH--TTSTT-SHHHHHH CS RrnaAD 87 kfetpeevte.....kdsihqeflvvanlPyeistqlvkqlle.......esrfglvkmllwlqkefarrllarpgskkrsrlsvlr 161 +f++ + e + + +++nlP +st+l+ ++le g+ +m+l +qke+a+rl+a+ gsk+rsrlsv++ tr|H2QTY1|H2QTY1_PANTR 114 TFKVEKAFSEslkrpWEDDPPNVHIIGNLPFSVSTPLIIKWLEniscrdgPFVYGRTQMTLTFQKEVAERLAANTGSKQRSRLSVMA 200 ******998888887555667899*********************999885566789****************************** PP HHHEEEEEEEEC-GGCECS--SSSEEEEEEEE-SS-S--S-HHHHHHHHHHHHCTTTSBCCCHCCCTSTCHHHHHCTT.T--T CS RrnaAD 162 eaftdvklvakveksifsPppkvdsalvelerkdedlpqkdlkkyesvvrkllnrkrktlstslksllpggelqalsstgind 244 +++++v+++ +++ ++f P+p+vd +v ++ ++ ++ k e+vv+ ++ +rk l l+p++++ + ++ + tr|H2QTY1|H2QTY1_PANTR 201 QYLCNVRHIFTIPGQAFVPKPEVDVGVVHFTPLIQPKIEQPFKLVEKVVQNVFQFRRKYCHRGLRMLFPEAQRLESTG-RLLE 282 *******************************97665347788999**************************9876644.4554 PP >> tr|H2Q1G2|H2Q1G2_PANTR rRNA adenine N(6)-methyltransferase OS=Pan troglodytes OX=9598 GN=TFB2M PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 70.2 0.0 4.6e-23 1.8e-19 30 263 .. 92 374 .. 69 376 .. 0.87 Alignments for each domain: == domain 1 score: 70.2 bits; conditional E-value: 4.6e-23 TTSEEEEES-TTTCCHHHHHCC-SEEEEE-S-HHHHHHHHCHTT.TCCCEEEE-S-CCGHHHHHHHCT................... CS RrnaAD 30 ekltvleigpGkGilTeeLakrakqvvvlElderLakllqkkls.kdeklkvvhqDflkfetpeevte................... 96 + +le pG+GilT+ L++ +vv+lE+d+ ++ +l+ + +d+kl v+h+Df k + + + tr|H2Q1G2|H2Q1G2_PANTR 92 PPHLLLECNPGPGILTQALLEAGAKVVALESDKTFIPHLESLGKnLDGKLRVIHCDFFKLDPRSGG-Vikppamssrglfknlgiea 177 566789**********************************97669*****************9943.2469***99*****999999 PP .TTCCCCEEEEEEE--TTTHHHHHHHHHH.......GGGGCEEEEEEEEEHHHHHHHH--TTS.TT-SHHHHHHHHHEEEEEEEEC- CS RrnaAD 97 .kdsihqeflvvanlPyeistqlvkqlle.......esrfglvkmllwlqkefarrllarpgs.kkrsrlsvlreaftdvklvakve 174 +++v++ P + +l ++fg+ ++ + + ++ ++l+a+pg+ ++ lsv+ ++++++k++ + tr|H2Q1G2|H2Q1G2_PANTR 178 vPWTADIPLKVIGMFPSRGEKRALWKLAYdlysctsIYKFGRIEVNMFIGEKEFQKLMADPGNpDLYHVLSVIWQLACEIKVLHMEP 264 98888999********965555444444458899999***********999999*******98468*****************9999 PP GG..................CECS--SSSEEEEEEEE-SS-S.--S...-HHHHHHHHHHHHCTTTSBCCCHCCCTSTCHHHHHCTT CS RrnaAD 175 ks..................ifsPppkvdsalvelerkdedl.pqk...dlkkyesvvrkllnrkrktlstslksllpggelqalss 239 s + + + ++ l++++++++ + +++ +++ + ++++++ r++ t+ l sl+p ++++ l + tr|H2Q1G2|H2Q1G2_PANTR 265 WSsfdiytrkgplenpkrreL-LDQLQQKLYLIQMTPRQNLFtKNLtpmNYNIFFHLLKHCFGRRSATVIDHLRSLTPLDARDILMQ 350 999**************9953.7788899***********99776533378888999**************************9987 PP .T--TTTCCCCSGHHHHHHHHHHH CS RrnaAD 240 tgindnalvkklsaeqtldifkel 263 g +++ +v++ ++++ ++f+ + tr|H2Q1G2|H2Q1G2_PANTR 351 IGKQEDEKVVNMHPQDFKTLFETI 374 799999*************99987 PP >> tr|H2QQ96|H2QQ96_PANTR Protein arginine methyltransferase 9 OS=Pan troglodytes OX=9598 GN=PRMT9 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.6 0.0 5.4e-07 0.0021 16 94 .. 152 234 .. 146 258 .. 0.87 Alignments for each domain: == domain 1 score: 17.6 bits; conditional E-value: 5.4e-07 HHHHHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-..SEEEEE-S-HHHHHHHHCHTT...TCCCEEEE-S-CCGHHHHHHH CS RrnaAD 16 kkvineIldklnleekltvleigpGkGilTeeLakra..kqvvvlElderLakllqkkls...kdeklkvvhqDflkfetpeev 94 ++++n +++k+ ++ vl+ig+G+Gil+ +ak a +v + El + +++l + + + +k++h+ l+ e+p+ + tr|H2QQ96|H2QQ96_PANTR 152 NTIYNAAIQKAVCLGSKSVLDIGAGTGILSM-FAKKAgaHSVYACELSKTMYELACDVVAankMEAGIKLLHTKSLDIEIPKHI 234 6799*************************97.56666589999**********9999877776788899999999999999965 PP >> tr|H2QHP4|H2QHP4_PANTR tRNA methyltransferase 61B OS=Pan troglodytes OX=9598 GN=TRMT61B PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.7 0.0 8.2e-06 0.032 17 69 .. 232 287 .. 228 321 .. 0.82 Alignments for each domain: == domain 1 score: 13.7 bits; conditional E-value: 8.2e-06 HHHHHHHHHH-GGTTSEEEEES-TTTCCHHHHHCC-...SEEEEE-S-HHHHHHHH CS RrnaAD 17 kvineIldklnleekltvleigpGkGilTeeLakra...kqvvvlElderLakllq 69 k in Il+ +++++ +tvle g+G G ++ L k + v ++E +++ l + tr|H2QHP4|H2QHP4_PANTR 232 KDINMILSMMDINPGDTVLEAGSGSGGMSLFLSKAVgsqGRVISFEVRKDHHDLAK 287 6699****************************999987777899998876444444 PP >> tr|A0A2I3RRF6|A0A2I3RRF6_PANTR Ribosomal RNA adenine dimethylase domain containing 1 OS=Pan troglodytes OX=9598 GN=R # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.6 0.0 3.5e-05 0.13 32 70 .. 155 194 .. 143 218 .. 0.85 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 3.5e-05 SEEEEES-TTTCCHHHHHCC-.SEEEEE-S-HHHHHHHHC CS RrnaAD 32 ltvleigpGkGilTeeLakra.kqvvvlElderLakllqk 70 ++v+++g+G+G+l++ +a + v ++E d+rL+++ q tr|A0A2I3RRF6|A0A2I3RRF6_PANTR 155 TQVVDVGSGQGHLSRFMALGLgLMVKSIEGDQRLVERAQR 194 79***************888879999**********9997 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (265 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 646 (0.0280772); expected 460.2 (0.02) Passed bias filter: 434 (0.018863); expected 460.2 (0.02) Passed Vit filter: 56 (0.00243394); expected 23.0 (0.001) Passed Fwd filter: 6 (0.000260779); expected 0.2 (1e-05) Initial search space (Z): 23008 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.25u 0.01s 00:00:00.26 Elapsed: 00:00:00.08 # Mc/sec: 39558.72 // [ok]