# 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/DNA_methylase.hmm.txt # target sequence database: proteomes/Pan_troglodytes.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: DNA_methylase [M=335] Accession: PF00145.17 Description: C-5 cytosine-specific DNA methylase 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.8e-47 161.0 0.0 5.8e-47 161.0 0.0 3.2 3 tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR DNA (cytosine-5)-methyltransfera 3.9e-42 145.1 0.0 4.9e-42 144.8 0.0 1.1 1 tr|H2R2B5|H2R2B5_PANTR TRDMT1 isoform 3 OS=Pan troglody 2.4e-16 60.3 0.3 1.9e-11 44.3 0.0 2.8 3 tr|A0A2I3SYU8|A0A2I3SYU8_PANTR DNA methyltransferase 3 beta OS= 5.2e-16 59.3 0.0 1.1e-11 45.0 0.0 2.3 2 tr|A0A2I3S8P5|A0A2I3S8P5_PANTR DNA methyltransferase 3 alpha OS 0.0013 18.5 0.0 0.0025 17.6 0.0 1.4 1 tr|A0A2I3SPD6|A0A2I3SPD6_PANTR Methyltransferase like 5 OS=Pan ------ inclusion threshold ------ 0.039 13.7 3.6 0.073 12.8 3.4 1.7 1 tr|A0A2I3TA51|A0A2I3TA51_PANTR EF-hand calcium binding domain 1 Domain annotation for each sequence (and alignments): >> tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR DNA (cytosine-5)-methyltransferase OS=Pan troglodytes OX=9598 GN=DNMT1 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -5.1 4.8 5.3 2e+04 258 294 .. 170 228 .. 26 233 .. 0.52 2 ? -19.7 22.4 6 2.3e+04 177 293 .. 194 347 .. 171 367 .. 0.36 3 ! 161.0 0.0 1.5e-50 5.8e-47 2 334 .. 1156 1612 .. 1155 1613 .. 0.75 Alignments for each domain: == domain 1 score: -5.1 bits; conditional E-value: 5.3 -..........................-B-TTS---EEEETTEEE--T.TTGGG--B-BHHHH CS DNA_methylase 258 e..........................klekelaeagklqknkerfkpviekkriRrLtprEa 294 + e e a++ +++k++ + ++++ Rr t rE+ tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 170 SpritrkstrqttitshfakgpakrkpQEESERAKSDESIKEEDK----DQDEKRRRVTSRER 228 033322222223333334444433333444444555555555443....45688888888886 PP == domain 2 score: -19.7 bits; conditional E-value: 6 S---GGG.....TSGGGTTS-EEBTTTB--GGGSSSTT-EE--.......S---H..HHTS--................ CS DNA_methylase 177 rfpkkkdleasktlrdlleekseldeektalsdklvekeekikkttkaseegypk..elsrnr................ 237 r+p+++ +e k ++++ee ++ de++ +++ +++ ++ ++g ++ e+ r+ tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 194 RKPQEE-SERAKSDESIKEEDKDQDEKRRRVTSRERVARPLPAEEPERAKSGTRTekEEERDEkeekrlrsqtkeptpk 271 444444.233444444444433333333333333333333333333333333332211111111111211222112221 PP ................B--TTS...-B----S-STT-.........................-B-TTS---EEEETTEE CS DNA_methylase 238 ................idkveelkglsatrrksgrge.........................klekelaeagklqknke 275 + + k+ e tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 272 qklkeepdrearagvqA-------------------DededgdekdekkhrsqpkdlaakrrP----------EEKEPE 321 11112222222333331...................022222222222222211112222221..........222222 PP E--T.........TTGGG--B-BHHH CS DNA_methylase 276 rfkp........viekkriRrLtprE 293 + +p ++++ R+ tp+E tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 322 KVNPqisdekdeDEKEEKRRKTTPKE 347 22223333333334555555555555 PP == domain 3 score: 161.0 bits; conditional E-value: 1.5e-50 EEEEES-TT-HHHHHHHHTT.-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-............GS--... CS DNA_methylase 2 kfidLfaGiGGfrlgleaag.fecvaaseidkkaaktYeanfk..kvsigDitkidkk............klpe..d 61 +++d f+G+GG++ g+++ag ++ +a+e+++ aa+ ++ n++ v ++D + + k lp+ d tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 1156 RTLDVFSGCGGLSEGFHQAGiSDTLWAIEMWDPAAQAFRLNNPgsTVFTEDCNILLKLvmagettnsrgqRLPQkgD 1232 789*****************999*********************745555555555446666777777779998667 PP -SEEEE---GGGTECTSE..E--TTSGGGHHHHHHHHHHHHH--SEEEEEEEGGGGGGCCSHHHHHHHHHHHHHTEE CS DNA_methylase 62 iDlLvgGfPCqpfSvaGk..qkgfedtrgtLffeiiriikekkPkvfllENVkgLlkhdkgktlkvilekleelGYe 136 +++L+gG+PCq+fS ++ + +++L + + + + ++P +fllENV++ ++ +++ +lk l l +GY+ tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 1233 VEMLCGGPPCQGFSGMNRfnSRTYSKFKNSLVVSFLSYCDYYRPRFFLLENVRNFVSFKRSMVLKLTLRCLVRMGYQ 1309 **************887654666677777************************************************ PP EEEEEEEGGGTT-SB--EEEEEEEEEGGG---...............-------SS---GGG..............T CS DNA_methylase 137 vsykllnasdyGvpQnReRvfivgirkdleln...............vlvkvkefrfpkkkd.........leaskt 189 +++ +l+a +yGv+Q R+R +i++ +++l ++v v +++f ++ ++++t tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 1310 CTFGVLQAGQYGVAQTRRRAIILAAAPGEKLPlfpeplhvfapracqLSVVVDDKKFVSNITrlssgpfrtITVRDT 1386 ***************************99888888888777777766333333444433333566676555555556 PP SGGGTTS-EEBTTTB--GGGSSSTT-EE--.......S---HHHTS--............................. CS DNA_methylase 190 lrdlleekseldeektalsdklvekeekikkttkaseegypkelsrnr............................. 237 + dl+e + + + + ++ e +++ +++ + +y+ + r+ tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 1387 MSDLPEVR-----NGASALEISYNGEPQSWFQRQLRGAQYQ-PILRDHickdmsalvaarmrhiplapgsdwrdlpn 1457 66665543.....2222222233333333333333333444.22222234444666677788999999999999999 PP ....B--TTS...-B----S-STT-.................................................... CS DNA_methylase 238 ....idkveelkglsatrrksgrge.................................................... 258 + ++ ++++l+ t+++ ++g+ tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 1458 ievrLSDGTMARKLRYTHHDRKNGRsssgalrgvcscveggasgkacdpaarqfntlipwclphtgnrhnhwaglyg 1534 98884444444555555554444445566678888888888888888888887777777777777777777777777 PP ......-B-TTS---EEEETTEEE--T.TTGGG--B-BHHHHHHHHT--TT-----SBHHHHHHHHTS--HHHHHHH CS DNA_methylase 259 ......klekelaeagklqknkerfkpviekkriRrLtprEaarlqGfPddfefgskksaaykqiGnsVpvpvieai 329 ++ +++ ++++k+++ v++++++R ++rE+ar qGfPd++ + ++ ++ +q+Gn+Vp+p+++ai tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 1535 rlewdgFFSTTVTNPEPMGKQGR----VLHPEQHRVVSVRECARSQGFPDTYRLFGNILDKHRQVGNAVPPPLAKAI 1607 66666555666666666666666....9************************************************* PP HHHHH CS DNA_methylase 330 akeik 334 + eik tr|A0A2I3T9Z6|A0A2I3T9Z6_PANTR 1608 GLEIK 1612 99986 PP >> tr|H2R2B5|H2R2B5_PANTR TRDMT1 isoform 3 OS=Pan troglodytes OX=9598 GN=TRDMT1 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 144.8 0.0 1.3e-45 4.9e-42 2 332 .. 5 388 .. 4 391 .] 0.74 Alignments for each domain: == domain 1 score: 144.8 bits; conditional E-value: 1.3e-45 EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-GS--..-SEEEE---GGGTECTSEE--T CS DNA_methylase 2 kfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfk..kvsigDitkidkkklpe.diDlLvgGfPCqpfSvaGkqkgf 83 ++++L +G+GG++++l++ ++vaa++ + a ++Y+ nf+ + + i+ i+ +++ ++D+++ ++PCqpf G q ++ tr|H2R2B5|H2R2B5_PANTR 5 RVLELYSGVGGMHHALRESCipAQVVAAIDVNTVANEVYKYNFPhtQLLAKTIEGITLEEFDRlSFDMILMSPPCQPFTRIGRQGDM 91 899***************99999*********************777777778888888999999********************** PP TSGGGHHHHHHHHHHHHH-..-SEEEEEEEGGGGGGCCSHHHHHHHHHHHHHTEEEEEEEEEGGGTT-SB--EEEEEEEEEGGG--- CS DNA_methylase 84 edtrgtLffeiiriikekk..PkvfllENVkgLlkhdkgktlkvilekleelGYevsykllnasdyGvpQnReRvfivgirkdleln 168 d r + f +i i+ + Pk++llENVkg + t + ++++ e++G++++ ll G+p +R R f+++ +++++l tr|H2R2B5|H2R2B5_PANTR 92 TDSRTNSFLYILDILPRLQklPKYILLENVKGFEVS---STRDLLIQTIENCGFQYQEFLLSPTSLGIPNSRLRYFLIAKLQSEPLP 175 ***************998888***********9888...567889***************************************999 PP ........-------SS---GGG.....TSGGGTTS-EE.....BTTTB--GGGSSSTT-EE--.......S---H..HHTS--... CS DNA_methylase 169 ........vlvkvkefrfpkkkdleasktlrdlleekse.....ldeektalsdklvekeekikkttkaseegypk..elsrnr... 237 ++++ e+ +p+k+ ++ ++++e+++e + + + +d+++ k e+ ++ + +++++ + + + tr|H2R2B5|H2R2B5_PANTR 176 fqapgqvlMEFPKIESVHPQKYA---MDVENKIQEKNVEpnisfDGSIQCSGKDAILFKLETAEEIHRKNQQDSDLsvK----Mlkd 255 99999977666666666666654...33444444443333344433334444444444444444333333333331222....1344 PP .......................................B--TTS...-B----S-STT-.....-B-TTS---EEEETTEEE--T. CS DNA_methylase 238 .......................................idkveelkglsatrrksgrge.....klekelaeagklqknkerfkpv 280 + g+ ++++ +g+ + ++++++ ++ +++e+ tr|H2R2B5|H2R2B5_PANTR 256 fleddtdvnqyllppksllryallldivqptcrrsvcftK-------GYGSYIEGTGSVLqtaedVQVENIYKSLTNLSQEEQIT-K 334 4677778888888888888888888888888855444333.......44444444443223322222333333333333333333.3 PP TTGGG--B-BHHHHHHHHT--TT-----SBH..HHHHHHHTS--HHHHHHHHHH CS DNA_methylase 281 iekkriRrLtprEaarlqGfPddfefgskks..aaykqiGnsVpvpvieaiake 332 + + R +tp+E a l GfP +f f++k + ++y+++Gns++v v+ +++k tr|H2R2B5|H2R2B5_PANTR 335 LLILKLRYFTPKEIANLLGFPPEFGFPEKITvkQRYRLLGNSLNVHVVAKLIKI 388 555678*******************************************99985 PP >> tr|A0A2I3SYU8|A0A2I3SYU8_PANTR DNA methyltransferase 3 beta OS=Pan troglodytes OX=9598 GN=DNMT3B PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.0 0.1 0.3 1.1e+03 219 275 .. 349 404 .. 306 428 .. 0.73 2 ! 44.3 0.0 4.9e-15 1.9e-11 1 120 [. 575 708 .. 575 725 .. 0.84 3 ! 12.9 0.0 1.7e-05 0.066 258 331 .. 769 844 .. 727 848 .. 0.77 Alignments for each domain: == domain 1 score: -1.0 bits; conditional E-value: 0.3 -.......S---HHHTS--B--TTS...-B----S-STT--B-TTS---EEEETTEE CS DNA_methylase 219 kkttkaseegypkelsrnridkveelkglsatrrksgrgeklekelaeagklqknke 275 k++++++ + + ++ r +k e+ k++++trr++ ++ ++++ ++ +l n + tr|A0A2I3SYU8|A0A2I3SYU8_PANTR 349 KPNNTQPVVNKS-KVRRAGSRKLESRKYENKTRRRTADDSATSDYCPAPKRLKTNCY 404 455555555566.67777789999999999999999988766666666555555544 PP == domain 2 score: 44.3 bits; conditional E-value: 4.9e-15 -EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-GS--...-SEEEE---GGG CS DNA_methylase 1 lkfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfk..kvsigDitkidkkklpe.d.iDlLvgGfPCqp 73 +++++Lf+Gi l l++ g + ++ase ++ +++ +++ + +D+++i+kk+++e +Dl++gG PC++ tr|A0A2I3SYU8|A0A2I3SYU8_PANTR 575 IRVLSLFDGIATGYLVLKELGikVGKYVASEVCEESIAVGTVKHEgnIKYVNDVRNITKKNIEEwGpFDLVIGGSPCND 653 5899*****************887778999999997777777776669999*********9998647************ PP TECTS.EE--TTSGGGHHHHHHHHHHHHH-.......-SEEEEEEEGGGGGGCCS CS DNA_methylase 74 fSvaG.kqkgfedtrgtLffeiiriikekk.......PkvfllENVkgLlkhdkg 120 S ++kg + g Lffe+ +++ + P + ++ENV ++ dk+ tr|A0A2I3SYU8|A0A2I3SYU8_PANTR 654 LSNVNpARKGLYEGTGRLFFEFYHLLNYSRpkegddrPFFWMFENVVAMKVGDKR 708 **87626********************9995555555668899999988776665 PP == domain 3 score: 12.9 bits; conditional E-value: 1.7e-05 --B-TTS---EEEETTEEE--T.TTGGG--B-BHHHHHHHHT--TT-----SBH..HHHHHHHTS--HHHHHHHHH CS DNA_methylase 258 eklekelaeagklqknkerfkpviekkriRrLtprEaarlqGfPddfefgskks..aaykqiGnsVpvpvieaiak 331 +k+++ ++++++ + k++ pv+ + ++ L E+ r+ GfP ++ s+ a+ k++G s +vpvi++++ tr|A0A2I3SYU8|A0A2I3SYU8_PANTR 769 KKVQTITTKSNSIKQGKNQLFPVVMNGKEDVLWCTELERIFGFPVHYTDVSNMGrgARQKLLGRSWSVPVIRHLFA 844 25556666666666666666679999*********************999888888*****************986 PP >> tr|A0A2I3S8P5|A0A2I3S8P5_PANTR DNA methyltransferase 3 alpha OS=Pan troglodytes OX=9598 GN=DNMT3A PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 45.0 0.0 2.9e-15 1.1e-11 1 120 [. 623 756 .. 623 765 .. 0.88 2 ! 12.1 0.0 3.1e-05 0.12 262 332 .. 821 893 .. 773 896 .. 0.76 Alignments for each domain: == domain 1 score: 45.0 bits; conditional E-value: 2.9e-15 -EEEEES-TT-HHHHHHHHTT..-EEEEEEE--TTTHHHHHHH-...-SEEEES-TTTS-GS--...-SEEEE---GGG CS DNA_methylase 1 lkfidLfaGiGGfrlgleaag..fecvaaseidkkaaktYeanfk..kvsigDitkidkkklpe.d.iDlLvgGfPCqp 73 +++++Lf+Gi l l+ g ++ ++ase ++ +++ +++ + gD++ +++k ++e +Dl++gG PC++ tr|A0A2I3S8P5|A0A2I3S8P5_PANTR 623 IRVLSLFDGIATGLLVLKDLGiqVDRYIASEVCEDSITVGMVRHQgkIMYVGDVRSVTQKHIQEwGpFDLVIGGSPCND 701 5899*****************99999*********999999999988999**********9998647************ PP TECTS.EE--TTSGGGHHHHHHHHHHHHH-.......-SEEEEEEEGGGGGGCCS CS DNA_methylase 74 fSvaG.kqkgfedtrgtLffeiiriikekk.......PkvfllENVkgLlkhdkg 120 S+ ++kg + g Lffe+ r++ + + P + l+ENV + +dk+ tr|A0A2I3S8P5|A0A2I3S8P5_PANTR 702 LSIVNpARKGLYEGTGRLFFEFYRLLHDARpkegddrPFFWLFENVVAMGVSDKR 756 **98736*********************995555555568899999888776664 PP == domain 2 score: 12.1 bits; conditional E-value: 3.1e-05 TTS---EEEETTEEE--T.TTGGG--B-BHHHHHHHHT--TT-----SBH..HHHHHHHTS--HHHHHHHHHH CS DNA_methylase 262 kelaeagklqknkerfkpviekkriRrLtprEaarlqGfPddfefgskks..aaykqiGnsVpvpvieaiake 332 + ++++++ + k + pv ++++ L E+ r+ GfP ++ s+ s a+ +++G s +vpvi++++ tr|A0A2I3S8P5|A0A2I3S8P5_PANTR 821 TITTRSNSIKQGKDQHFPVFMNEKEDILWCTEMERVFGFPVHYTDVSNMSrlARQRLLGRSWSVPVIRHLFAP 893 33333333333333334588899999999***************999999889999*************9875 PP >> tr|A0A2I3SPD6|A0A2I3SPD6_PANTR Methyltransferase like 5 OS=Pan troglodytes OX=9598 GN=METTL5 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.5e-07 0.0025 3 70 .. 55 128 .. 54 129 .. 0.84 Alignments for each domain: == domain 1 score: 17.6 bits; conditional E-value: 6.5e-07 EEEES-TT-HHHHHHHHTT-EEEEEEE--TTTHHHHHHH-.-SEEEES-TTTS-......GS--.-SEEEE--- CS DNA_methylase 3 fidLfaGiGGfrlgleaagfecvaaseidkkaaktYeanfkkvsigDitkidkk......klpediDlLvgGfP 70 + dL +G+G +++g ++ g ++ +id++a +++++n ++ ++ +i+ ++ + +++++D ++ +P tr|A0A2I3SPD6|A0A2I3SPD6_PANTR 55 VADLGCGCGVLSIGTAMLGAGLCVGFDIDEDALEIFNRNAEEFELTNIDMVQCDvcllsnRMSKSFDTVIMNPP 128 579**************************************888888887777666666666656777776665 PP >> tr|A0A2I3TA51|A0A2I3TA51_PANTR EF-hand calcium binding domain 12 OS=Pan troglodytes OX=9598 GN=EFCAB12 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.8 3.4 1.9e-05 0.073 144 259 .. 47 174 .. 41 333 .. 0.70 Alignments for each domain: == domain 1 score: 12.8 bits; conditional E-value: 1.9e-05 GGGTT-SB--EEEEEEEEEGGG----------SS---GGG......TSGGGTTS-EEBTTTB--GGGSSSTT-EE--.. CS DNA_methylase 144 asdyGvpQnReRvfivgirkdlelnvlvkvkefrfpkkkd.leasktlrdlleekseldeektalsdklvekeekikkt 221 +d+ pQ R+R+++v ++d++ ++++ ++ + p k + + + rd++e++ d + + +++ + +e ++ tr|A0A2I3TA51|A0A2I3TA51_PANTR 47 QKDFRLPQTRRRIIMVPRKEDQT--PINPASQPQAPPKPIpSFKVLEARDIQEQP--EDRKTWLSQRSKLRQELESFGD 121 589******************97..456666655555533344556667777776..5666666666666666666667 PP .....S---H.........HHTS--B--TTS...-B----S-STT-........- CS DNA_methylase 222 tkaseegypk.........elsrnridkveelkglsatrrksgrge........k 259 k+++e+ p+ ++ +++ + + +at r++++++ + tr|A0A2I3TA51|A0A2I3TA51_PANTR 122 VKRWLENKPSitpseakvlHMIH--EEQSAQPNASQATTRTTRKKAprlsrlsrQ 174 77777777754444444332222..244445555555555555433333333331 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (335 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 1777 (0.077234); expected 460.2 (0.02) Passed bias filter: 708 (0.0307719); expected 460.2 (0.02) Passed Vit filter: 61 (0.00265125); expected 23.0 (0.001) Passed Fwd filter: 11 (0.000478095); 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.46u 0.02s 00:00:00.48 Elapsed: 00:00:00.15 # Mc/sec: 26671.04 // [ok]