# 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]