Citrus Sinensis ID: 026734

---

Local Sequence Feature Prediction

Residue Number Marker
--------10--------20--------30--------40--------50--------60--------70--------80--------90-------100-------110-------120-------130-------140-------150-------160-------170-------180-------190-------200-------210-------220-------230----
Protein Sequence
MGHMHGFFPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQAG
Secondary Structure (PSIPRED)
ccccccccccccccccccccccccccccccccccccccccccccccccEEEcccccccccHHHHcccccccHHHHccccccccccccccccccccccccccHHHHHHHHcccccccccccccHHHHHHHccccccccccccccEEEEEccEEEEEcccccEEEEccccccccccccccccccccHHHHHHHccccccccccccccHHHHHHHHcccEEHHHHHHHHHHHHHHcc
Secondary Structure Prediction (SSPRO)
cccccccccccccHHHHHHHccccccccccccccccccccccccccccEEEcccccccccHHHHHHccccHHHHHHHcccccccccccccccccccccccccccHHHHHccccccHHHccccHHHHHHHHHHHcEEccccccEEEEEccccEEEcccccEEEccccccHHHHHHHHHcccccccHHHHHHHccccccccccccccHHHHHHHHcccccHHHHHHHHHHHHHHcc
Coil and Loop (DISEMBL)
mghmhgffpphvnptlgkakrkplSFRCQLLNnqllrspspllgnddmtvitkhdqpddswdkllescdpverflefsnsgdqvntetgflstgtaavddfgaaeetvevdrcvsidhflvPLSLIERWGsamdivypdskrcccftksyYRYVKGTGsllatvqpknkgkasslkeqhlryftprevanlhsfpgdfqfphhlSLRQRYALLGNSLSIAVVAPLLQYLFAQAG
Flexible Loop (DISEMBL)
mghmhgffpphvnptlgkAKRKPLSFRCQLLNNqllrspspllGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLatvqpknkgkasslkEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQAG
Low Complexity Region (SEG)
MGHMHGFFPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQAG
Disordered region (IsUnstruct)
***********************LSFRCQLLNNQLL************************WDKLLESCDPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATV**************HLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA***
Disordered Region (DISOPRED)
*****GFFPPHVNPTLGKAKRKP******LLNNQLLRSPSPLLG*********************ESCDPVERFLEFSNSGDQVN********************ETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKG***SLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQA*
Disordered Region (DISEMBL)
MGHMHGFFPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQ**********KEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQAG
Disordered Region (DISPRO)
****HGFFPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSN***********************AAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQA*
Transmembrane Helix (TMHMM)
oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Transmembrane Helix (HMMTOP)
ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooHHHHHHHHHHHHHHHHHHHHiiii
Transmembrane Helix (MEMSAT)
ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooohhhhhhhhhhhhhhhhhhhiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
TM Helix, Signal Peptide (MEMSAT_SVM)
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiihhhhhhhhhhhhhhhhhooooo
TM Helix, Signal Peptide (Phobius)
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiHHHHHHHHHHHHHHHHHHHHHoooo
Signal Peptide (SignalP HMM Mode)
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
Signal Peptide (SignalP NN Mode)
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
Coiled Coils (COILS)
MGHMHGFFPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQAG
Positional Conservation
no confident homologs detected

---

Close Homologs for Annotation Transfer

Close Homologs in SWISS-PROT Database Detected by BLAST

Original result of BLAST against SWISS-PROT Database

ID	Length	Definition	RBH(Q2H)	RBH(H2Q)	Q cover	H cover	Identity	E-value
Query	234	2.2.26 [Sep-21-2011]
O14717	391	tRNA (cytosine(38)-C(5))-	yes	no	0.491	0.294	0.437	4e-22
Q7YS61	391	tRNA (cytosine-5-)-methyl	yes	no	0.470	0.281	0.475	5e-22
Q4G073	391	tRNA (cytosine(38)-C(5))-	yes	no	0.478	0.286	0.44	2e-21
O55055	415	tRNA (cytosine(38)-C(5))-	yes	no	0.470	0.265	0.449	4e-20
Q54JH6	379	DNA (cytosine-5)-methyltr	yes	no	0.478	0.295	0.4	2e-16
P40999	330	tRNA (cytosine(38)-C(5))-	yes	no	0.517	0.366	0.373	4e-16

>sp|O14717|TRDMT_HUMAN tRNA (cytosine(38)-C(5))-methyltransferase OS=Homo sapiens GN=TRDMT1 PE=1 SV=1

 Score =  104 bits (260), Expect = 4e-22,   Method: Compositional matrix adjust.
 Identities = 56/128 (43%), Positives = 82/128 (64%), Gaps = 13/128 (10%)

Query: 116 IDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKNKGKA 172
           ++ +L+P   + R+   +DIV P  +R  CFTK Y  Y++GTGS+L T   VQ +N  K+
Sbjct: 263 VNQYLLPPKSLLRYALLLDIVQPTCRRSVCFTKGYGSYIEGTGSVLQTAEDVQVENIYKS 322

Query: 173 S---SLKEQ-------HLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVV 222
               S +EQ        LRYFTP+E+ANL  FP +F FP  ++++QRY LLGNSL++ VV
Sbjct: 323 LTNLSQEEQITKLLILKLRYFTPKEIANLLGFPPEFGFPEKITVKQRYRLLGNSLNVHVV 382

Query: 223 APLLQYLF 230
           A L++ L+
Sbjct: 383 AKLIKILY 390

Specifically methylates cytosine 38 in the anticodon loop of tRNA(Asp).
Homo sapiens (taxid: 9606)
EC: 2EC: .EC: 1EC: .EC: 1EC: .EC: 2EC: 0EC: 4

>sp|Q7YS61|TRDMT_BOVIN tRNA (cytosine-5-)-methyltransferase OS=Bos taurus GN=TRDMT1 PE=2 SV=1

 Score =  104 bits (260), Expect = 5e-22,   Method: Compositional matrix adjust.
 Identities = 59/124 (47%), Positives = 82/124 (66%), Gaps = 14/124 (11%)

Query: 119 FLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKNKGKA-SS 174
           FL P SL+ R+   +DIV P S+R  CFTK Y RY++GTGS+L T   VQ +N  K+ +S
Sbjct: 267 FLPPKSLL-RYALLLDIVKPTSRRSMCFTKGYGRYIEGTGSVLQTTEDVQIENIYKSLTS 325

Query: 175 LKEQ---------HLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPL 225
           L ++          LR+FTP+E+ANL  FP +F FP   +++QRY LLGNSL++ VVA L
Sbjct: 326 LSQEEKIMRLSMLQLRFFTPKEIANLLGFPPEFGFPEMTTVKQRYRLLGNSLNVHVVAKL 385

Query: 226 LQYL 229
           ++ L
Sbjct: 386 IKIL 389

Specifically methylates cytosine 38 in the anticodon loop of tRNA(Asp).

Bos taurus (taxid: 9913)

EC: 2

EC: .

EC: 1

EC: .

EC: 1

EC: .

EC: 2

EC: 0

EC: 4

>sp|Q4G073|TRDMT_RAT tRNA (cytosine(38)-C(5))-methyltransferase OS=Rattus norvegicus GN=Trdmt1 PE=2 SV=1

 Score =  102 bits (254), Expect = 2e-21,   Method: Compositional matrix adjust.
 Identities = 55/125 (44%), Positives = 76/125 (60%), Gaps = 13/125 (10%)

Query: 118 HFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKN------ 168
            +L+P   + R+   +DIV P S+R  CFTK Y  Y++GTGS+L T   VQ +N      
Sbjct: 265 QYLLPAKSLLRYALLLDIVKPTSRRSMCFTKGYGSYIEGTGSVLQTAEDVQIENIYKSLP 324

Query: 169 ----KGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAP 224
               + K + L    LRYFTP+E+ANL  FP +F FP   +++QRY LLGNSL++ VV+ 
Sbjct: 325 DLPPEEKIAKLSMLKLRYFTPKEIANLLGFPPEFGFPEKTTVKQRYRLLGNSLNVHVVSK 384

Query: 225 LLQYL 229
           LL  L
Sbjct: 385 LLTVL 389

Specifically methylates cytosine 38 in the anticodon loop of tRNA(Asp).

Rattus norvegicus (taxid: 10116)

EC: 2

EC: .

EC: 1

EC: .

EC: 1

EC: .

EC: 2

EC: 0

EC: 4

>sp|O55055|TRDMT_MOUSE tRNA (cytosine(38)-C(5))-methyltransferase OS=Mus musculus GN=Trdmt1 PE=2 SV=2

 Score = 98.2 bits (243), Expect = 4e-20,   Method: Compositional matrix adjust.
 Identities = 58/129 (44%), Positives = 77/129 (59%), Gaps = 19/129 (14%)

Query: 117 DHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQ----------- 165
           D +L+P  L+ R+   +DIV P S+R  CFTK Y  Y++GTGS+L   +           
Sbjct: 264 DEYLLPPKLLLRYALLLDIVKPTSRRSMCFTKGYGSYIEGTGSVLQAAEDAQIENIYKSL 323

Query: 166 ----PKNK-GKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIA 220
               P+ K  K S LK   LRYFTP+E+ANL  FP +F FP   +++QRY LLGNSL++ 
Sbjct: 324 PDLPPEEKIAKLSMLK---LRYFTPKEIANLQGFPPEFGFPEKTTVKQRYRLLGNSLNVH 380

Query: 221 VVAPLLQYL 229
           VVA LL  L
Sbjct: 381 VVAKLLTVL 389

Specifically methylates cytosine 38 in the anticodon loop of tRNA(Asp).

Mus musculus (taxid: 10090)

EC: 2

EC: .

EC: 1

EC: .

EC: 1

EC: .

EC: 2

EC: 0

EC: 4

>sp|Q54JH6|CMT1_DICDI DNA (cytosine-5)-methyltransferase OS=Dictyostelium discoideum GN=dnmA PE=1 SV=1

 Score = 85.9 bits (211), Expect = 2e-16,   Method: Compositional matrix adjust.
 Identities = 48/120 (40%), Positives = 71/120 (59%), Gaps = 8/120 (6%)

Query: 117 DHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKA---- 172
           + + VP +L+   G   DI   DSK   C TKSY ++++GTGS+   +Q  N  KA    
Sbjct: 245 EKYKVPHNLLLSKGMLFDIKQKDSKTSNCVTKSYGKFIEGTGSI---IQMDNNFKADIND 301

Query: 173 -SSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA 231
             SL    LRYF+P+E+  LH FP +F+F   L+  Q Y L+GNSL++ +V+ LL+ L +
Sbjct: 302 NKSLIPLKLRYFSPKEITRLHGFPEEFKFSPKLTTIQCYRLIGNSLNVKIVSELLKVLVS 361

Involved in epigenetic gene silencing. Methylates specific cytosine residues in the retrotransposons DIRS-1 and Skipper.

Dictyostelium discoideum (taxid: 44689)

EC: 2

EC: .

EC: 1

EC: .

EC: 1

EC: .

EC: 3

EC: 7

>sp|P40999|PMT1M_SCHPO tRNA (cytosine(38)-C(5))-methyltransferase OS=Schizosaccharomyces pombe (strain 972 / ATCC 24843) GN=pmt1 PE=1 SV=1

 Score = 85.1 bits (209), Expect = 4e-16,   Method: Compositional matrix adjust.
 Identities = 47/126 (37%), Positives = 71/126 (56%), Gaps = 5/126 (3%)

Query: 108 VEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPK 167
           +E++R  S   ++V  S++ +WG   DIV PDS  CCCFT+ Y   V+G GS+L     +
Sbjct: 202 LEIERDWS--SYMVLESVLNKWGHQFDIVKPDSSSCCCFTRGYTHLVQGAGSILQMSDHE 259

Query: 168 NKGKASSLKEQ--HLRYFTPREVANLHSFPGDFQFPH-HLSLRQRYALLGNSLSIAVVAP 224
           N  +          LRYFT REVA L  FP   ++   +++ +  Y LLGNS+++ VV+ 
Sbjct: 260 NTHEQFERNRMALQLRYFTAREVARLMGFPESLEWSKSNVTEKCMYRLLGNSINVKVVSY 319

Query: 225 LLQYLF 230
           L+  L 
Sbjct: 320 LISLLL 325

Does not have a cytosine-5 methyltransferase activity due to the insertion of a Ser residue between the Pro-Cys motif found at the active site of C5 MTases. When this serine is deleted it becomes catalytically active and recognizes and methylates the sequence CC[AT]GG.

Schizosaccharomyces pombe (strain 972 / ATCC 24843) (taxid: 284812)

Close Homologs in the Non-Redundant Database Detected by BLAST

Original result of BLAST against Nonredundant Database

GI	Length	Definition	Q cover	H cover	Identity	E-value
Query	234
359476106	393	PREDICTED: tRNA (cytosine-5-)-methyltran	0.905	0.539	0.699	1e-77
296081999	366	unnamed protein product [Vitis vinifera]	0.799	0.510	0.663	1e-71
224092438	350	DNA methyltransferase [Populus trichocar	0.799	0.534	0.630	2e-68
297812763	383	DNA methyltransferase-2 [Arabidopsis lyr	0.871	0.532	0.643	1e-67
21553420	383	DNA methyltransferase PMT1-like protein	0.876	0.535	0.641	4e-67
18420929	383	DNA methyltransferase-2 [Arabidopsis tha	0.876	0.535	0.641	6e-67
449441930	385	PREDICTED: tRNA (cytosine(38)-C(5))-meth	0.888	0.540	0.613	2e-64
255548992	404	DNA (cytosine-5)-methyltransferase, puta	0.846	0.490	0.604	3e-64
363814463	385	uncharacterized protein LOC100792567 [Gl	0.901	0.548	0.612	9e-62
357445897	378	tRNA (cytosine-5-)-methyltransferase [Me	0.901	0.558	0.586	2e-60

>gi|359476106|ref|XP_002282454.2| PREDICTED: tRNA (cytosine-5-)-methyltransferase-like [Vitis vinifera]

 Score =  295 bits (755), Expect = 1e-77,   Method: Compositional matrix adjust.
 Identities = 151/216 (69%), Positives = 171/216 (79%), Gaps = 4/216 (1%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKPLSF  QL NNQLL +PSPL G++D ++I +HDQ +  WDKLL+SC+PVERFLEF 
Sbjct: 176 AKRKPLSFENQLFNNQLLSTPSPLFGHEDDSLIDEHDQSEKIWDKLLQSCEPVERFLEFK 235

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRC--VSIDHFLVPLSLIERWGSAMDIV 136
           N+ + + TE+ +L T + + D  G  E   E DR    S D F VPLSLIERWGSAMDIV
Sbjct: 236 NTINPLETESSYLDTLSVSTDASGELENN-ERDRSDSNSWDQFSVPLSLIERWGSAMDIV 294

Query: 137 YPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPG 196
           YPDSKRCCCFTKSYYRYVKGTGSLLAT Q K K KAS L+EQ LRYFTPREVANLHSFP 
Sbjct: 295 YPDSKRCCCFTKSYYRYVKGTGSLLATGQSK-KDKASLLEEQCLRYFTPREVANLHSFPE 353

Query: 197 DFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
           DF FP H++LRQRYALLGNSLS+AVVAPLL YLF Q
Sbjct: 354 DFHFPQHVNLRQRYALLGNSLSVAVVAPLLTYLFTQ 389

Source: Vitis vinifera

Species: Vitis vinifera

Genus: Vitis

Family: Vitaceae

Order: Vitales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|296081999|emb|CBI21004.3| unnamed protein product [Vitis vinifera]

 Score =  275 bits (703), Expect = 1e-71,   Method: Compositional matrix adjust.
 Identities = 142/214 (66%), Positives = 159/214 (74%), Gaps = 27/214 (12%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKPLSF  QL NNQLL +PSPL G++D ++I +HDQ +  WDKLL+SC+PVERFLEF 
Sbjct: 176 AKRKPLSFENQLFNNQLLSTPSPLFGHEDDSLIDEHDQSEKIWDKLLQSCEPVERFLEFK 235

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
           N+ + + T+                           S D F VPLSLIERWGSAMDIVYP
Sbjct: 236 NTINPLETDN--------------------------SWDQFSVPLSLIERWGSAMDIVYP 269

Query: 139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
           DSKRCCCFTKSYYRYVKGTGSLLAT Q K K KAS L+EQ LRYFTPREVANLHSFP DF
Sbjct: 270 DSKRCCCFTKSYYRYVKGTGSLLATGQSK-KDKASLLEEQCLRYFTPREVANLHSFPEDF 328

Query: 199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
            FP H++LRQRYALLGNSLS+AVVAPLL YLF Q
Sbjct: 329 HFPQHVNLRQRYALLGNSLSVAVVAPLLTYLFTQ 362

Source: Vitis vinifera

Species: Vitis vinifera

Genus: Vitis

Family: Vitaceae

Order: Vitales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|224092438|ref|XP_002309610.1| DNA methyltransferase [Populus trichocarpa] gi|222855586|gb|EEE93133.1| DNA methyltransferase [Populus trichocarpa]

 Score =  264 bits (674), Expect = 2e-68,   Method: Compositional matrix adjust.
 Identities = 135/214 (63%), Positives = 155/214 (72%), Gaps = 27/214 (12%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKPLSF  ++ NN+LL SPSPL  +++  V    DQP  +WD L+ SC+PVERFLEF 
Sbjct: 164 AKRKPLSFHHEVCNNKLLWSPSPLFEHNENNVTDGWDQPPVNWDTLIHSCEPVERFLEFK 223

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
           +S ++V   TG                             ++VPLSLIERW SAMDIVYP
Sbjct: 224 SSSNRVTGGTG---------------------------TDYIVPLSLIERWESAMDIVYP 256

Query: 139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
           DSKRCCCFTKSYYRYVKGTGSLLAT+QP  KGK SSLK+Q LRYFTPREVANLHSFP DF
Sbjct: 257 DSKRCCCFTKSYYRYVKGTGSLLATIQPNIKGKESSLKKQGLRYFTPREVANLHSFPDDF 316

Query: 199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
           QFP  +SLRQ YALLGNSLS+AVVAPLL+YLF+Q
Sbjct: 317 QFPEAISLRQCYALLGNSLSVAVVAPLLRYLFSQ 350

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|297812763|ref|XP_002874265.1| DNA methyltransferase-2 [Arabidopsis lyrata subsp. lyrata] gi|297320102|gb|EFH50524.1| DNA methyltransferase-2 [Arabidopsis lyrata subsp. lyrata]

 Score =  261 bits (668), Expect = 1e-67,   Method: Compositional matrix adjust.
 Identities = 137/213 (64%), Positives = 161/213 (75%), Gaps = 9/213 (4%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKPL F+ Q  NN+LL SPSPL G DD   + K  Q ++  DKLL+ C+PVE+FLE  
Sbjct: 177 AKRKPLPFKSQHSNNKLLWSPSPLYGCDDQVEVGKC-QAEEGLDKLLQFCEPVEKFLE-- 233

Query: 79  NSGDQVNTETGFLS-TGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVY 137
            S   V+ E   +  +   + D +G  +++V      S+  +LVP+SLIERWG+AMDIVY
Sbjct: 234 -SAAHVDGEPSSVDDSENGSKDCYGQEDDSVPD----SVHQYLVPVSLIERWGNAMDIVY 288

Query: 138 PDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGD 197
           PDSKRCCCFTKSYYRYVKGTGSLLATVQPK KGK S LKEQ LRYFTPREVAN HSFP D
Sbjct: 289 PDSKRCCCFTKSYYRYVKGTGSLLATVQPKIKGKESCLKEQRLRYFTPREVANFHSFPED 348

Query: 198 FQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
           F FP H+SLRQRYA+LGNSLS+AVVAPLL+YLF
Sbjct: 349 FVFPKHISLRQRYAMLGNSLSVAVVAPLLRYLF 381

Source: Arabidopsis lyrata subsp. lyrata

Species: Arabidopsis lyrata

Genus: Arabidopsis

Family: Brassicaceae

Order: Brassicales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|21553420|gb|AAM62513.1| DNA methyltransferase PMT1-like protein [Arabidopsis thaliana]

 Score =  260 bits (664), Expect = 4e-67,   Method: Compositional matrix adjust.
 Identities = 136/212 (64%), Positives = 154/212 (72%), Gaps = 7/212 (3%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKPL F+ Q  NN+LL SP PL G DD     K  Q ++  DKLLE C+PVE+FLE +
Sbjct: 177 AKRKPLPFKSQHSNNKLLWSPDPLYGRDDQVEFGKC-QAEEGLDKLLEFCEPVEKFLELA 235

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
                V+ E   +        D    E     D   S+  +LVP+SLIERWG+AMDIVYP
Sbjct: 236 A---HVDGEPSSVDDSENGSKDCCGQEGDSVPD---SVHQYLVPVSLIERWGNAMDIVYP 289

Query: 139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
           DSKRCCCFTKSYYRYVKGTGSLLATVQPK KGK S LKEQ LRYFTPREVAN HSFP DF
Sbjct: 290 DSKRCCCFTKSYYRYVKGTGSLLATVQPKIKGKESCLKEQRLRYFTPREVANFHSFPEDF 349

Query: 199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
           +FP H+SLRQRYA+LGNSLS+AVVAPLL+YLF
Sbjct: 350 EFPKHISLRQRYAMLGNSLSVAVVAPLLRYLF 381

Source: Arabidopsis thaliana

Species: Arabidopsis thaliana

Genus: Arabidopsis

Family: Brassicaceae

Order: Brassicales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|18420929|ref|NP_568474.1| DNA methyltransferase-2 [Arabidopsis thaliana] gi|332006070|gb|AED93453.1| DNA methyltransferase-2 [Arabidopsis thaliana]

 Score =  259 bits (662), Expect = 6e-67,   Method: Compositional matrix adjust.
 Identities = 136/212 (64%), Positives = 153/212 (72%), Gaps = 7/212 (3%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKPL F+ Q  NN+LL SP PL G DD     K  Q ++  DKLLE C PVE+FLE +
Sbjct: 177 AKRKPLPFKSQHSNNKLLWSPDPLYGRDDQVEFGKC-QAEEGLDKLLEFCKPVEKFLELA 235

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
                V+ E   +        D    E     D   S+  +LVP+SLIERWG+AMDIVYP
Sbjct: 236 A---HVDGEPSSVDDSENGSKDCCGQEGDSVPD---SVHQYLVPVSLIERWGNAMDIVYP 289

Query: 139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
           DSKRCCCFTKSYYRYVKGTGSLLATVQPK KGK S LKEQ LRYFTPREVAN HSFP DF
Sbjct: 290 DSKRCCCFTKSYYRYVKGTGSLLATVQPKIKGKESCLKEQRLRYFTPREVANFHSFPEDF 349

Query: 199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
           +FP H+SLRQRYA+LGNSLS+AVVAPLL+YLF
Sbjct: 350 EFPKHISLRQRYAMLGNSLSVAVVAPLLRYLF 381

Source: Arabidopsis thaliana

Species: Arabidopsis thaliana

Genus: Arabidopsis

Family: Brassicaceae

Order: Brassicales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|449441930|ref|XP_004138735.1| PREDICTED: tRNA (cytosine(38)-C(5))-methyltransferase-like [Cucumis sativus] gi|449514978|ref|XP_004164528.1| PREDICTED: tRNA (cytosine(38)-C(5))-methyltransferase-like [Cucumis sativus]

 Score =  251 bits (640), Expect = 2e-64,   Method: Compositional matrix adjust.
 Identities = 132/215 (61%), Positives = 156/215 (72%), Gaps = 7/215 (3%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGND-DMTVITKHDQPDDSWDKLLESCDPVERFLEF 77
           AKRKPL F  +L NNQLL SP+PLL +D ++  I  H     S DKLL SC+P+ RFL+ 
Sbjct: 175 AKRKPLFFEKELYNNQLLWSPNPLLDSDANLESIETHGSQATS-DKLLLSCEPICRFLDH 233

Query: 78  SNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVY 137
           SN  +    +   + +  A   +  A E+ +E       + +LVP SLI+RWGSAMDIVY
Sbjct: 234 SNHQELSFYDASSIQSEVALEKNKDADEQVIE-----PFNQYLVPSSLIDRWGSAMDIVY 288

Query: 138 PDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGD 197
           PDS+RCCCFTKSYYRYVKGTGSLLA  + +  GKA SLKEQ LRYFTPREVANLHSFP D
Sbjct: 289 PDSRRCCCFTKSYYRYVKGTGSLLAPFETEGMGKAHSLKEQKLRYFTPREVANLHSFPED 348

Query: 198 FQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
           FQFP H+ LRQRYALLGNSLSIAVVAPLL+YLF +
Sbjct: 349 FQFPQHIGLRQRYALLGNSLSIAVVAPLLRYLFTE 383

Source: Cucumis sativus

Species: Cucumis sativus

Genus: Cucumis

Family: Cucurbitaceae

Order: Cucurbitales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|255548992|ref|XP_002515552.1| DNA (cytosine-5)-methyltransferase, putative [Ricinus communis] gi|223545496|gb|EEF47001.1| DNA (cytosine-5)-methyltransferase, putative [Ricinus communis]

 Score =  250 bits (639), Expect = 3e-64,   Method: Compositional matrix adjust.
 Identities = 127/210 (60%), Positives = 154/210 (73%), Gaps = 12/210 (5%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKPL F+ + +N+QL+ S  PL G+D+ TV  ++DQ  +SWDKLL+SC+PV  FLEF 
Sbjct: 175 AKRKPLLFQHKYINHQLVWSLGPLFGHDEGTVADRYDQSPESWDKLLQSCEPVVSFLEFK 234

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
           NS + V+ ET F+    A  +D G  E T E +   ++D ++VP +LIERWGSAMDIVYP
Sbjct: 235 NSSNGVDAETDFV----AITNDSGDLENTGEGN-LTTMDQYIVPSNLIERWGSAMDIVYP 289

Query: 139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
           DS RCCCFTKSYYRYVKGTGSLLATVQPK KGK SSL E+ LRYFTPREVANLHSFP DF
Sbjct: 290 DSNRCCCFTKSYYRYVKGTGSLLATVQPKEKGKTSSLHEKGLRYFTPREVANLHSFPKDF 349

Query: 199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQY 228
            FP H+SLRQR       LSI+ V  +  +
Sbjct: 350 HFPEHISLRQR-------LSISFVGSMCMF 372

Source: Ricinus communis

Species: Ricinus communis

Genus: Ricinus

Family: Euphorbiaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|363814463|ref|NP_001242866.1| uncharacterized protein LOC100792567 [Glycine max] gi|255640237|gb|ACU20409.1| unknown [Glycine max]

 Score =  242 bits (618), Expect = 9e-62,   Method: Compositional matrix adjust.
 Identities = 131/214 (61%), Positives = 150/214 (70%), Gaps = 3/214 (1%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKP SF  + LN+QL++SP PL  + D T   + D   +    LL+SC P+E FLE  
Sbjct: 174 AKRKPSSFVNECLNSQLIQSPPPLFEHFD-TAADEDDSSKEDRQNLLQSCQPIENFLELK 232

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
           N    +  E+  L T  +  D      +    D   S+D + V  SLIERWGSAMD+VYP
Sbjct: 233 NPSSDIAVESVALRTDLSN-DAPRTLGKNNGHDEYESLDQYFVHPSLIERWGSAMDVVYP 291

Query: 139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
           DSKRCCCFTKSYYRYVKGTGSLLATVQP  + K +SLKEQ LRYFTPREVANLHSFP  F
Sbjct: 292 DSKRCCCFTKSYYRYVKGTGSLLATVQPVKRDK-TSLKEQCLRYFTPREVANLHSFPEGF 350

Query: 199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
            FP H+SLRQRYALLGNSLSIAVVAPLLQYLF Q
Sbjct: 351 VFPEHISLRQRYALLGNSLSIAVVAPLLQYLFTQ 384

Source: Glycine max

Species: Glycine max

Genus: Glycine

Family: Fabaceae

Order: Fabales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|357445897|ref|XP_003593226.1| tRNA (cytosine-5-)-methyltransferase [Medicago truncatula] gi|355482274|gb|AES63477.1| tRNA (cytosine-5-)-methyltransferase [Medicago truncatula]

 Score =  237 bits (605), Expect = 2e-60,   Method: Compositional matrix adjust.
 Identities = 126/215 (58%), Positives = 151/215 (70%), Gaps = 4/215 (1%)

Query: 19  AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
           AKRKP SF    LN QL++SP PL  + + T   + D   +    LL+SC P+E+FL   
Sbjct: 168 AKRKPSSFLNGCLNRQLIQSPRPLFEHFN-TAPKEDDLSLEDSQNLLQSCQPIEKFLVLK 226

Query: 79  NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
           N     N E+   ++ T+  +D      T       ++D + V  SL+ERWGSAMD+VYP
Sbjct: 227 NPNSDTNVESA--ASTTSLSNDTSRTSGTDNDHEYDTLDKYYVHPSLLERWGSAMDVVYP 284

Query: 139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
           DSKRCCCFTKSYYRYVKGTGSLLATVQP  + K +SLKEQ LRYFTPREVANLHSFP DF
Sbjct: 285 DSKRCCCFTKSYYRYVKGTGSLLATVQPVKRDK-TSLKEQCLRYFTPREVANLHSFPEDF 343

Query: 199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQA 233
           +FP H+SL+QRYALLGNSLSIAVVA LLQYLF +A
Sbjct: 344 EFPEHISLKQRYALLGNSLSIAVVASLLQYLFTEA 378

Source: Medicago truncatula

Species: Medicago truncatula

Genus: Medicago

Family: Fabaceae

Order: Fabales

Phylum: Streptophyta

Superkingdom: Eukaryota

---

Prediction of Gene Ontology (GO) Terms

Close Homologs with Gene Ontology terms Detected by BLAST

Original result of BLAST against Gene Ontology (AMIGO)

ID	Length	Definition	Q cover	H cover	Identity	E-value
Query	234
TAIR|locus:2179469	383	DNMT2 "DNA methyltransferase-2	0.876	0.535	0.646	1.1e-64
UNIPROTKB|E7EMI8	310	TRDMT1 "tRNA (cytosine(38)-C(5	0.876	0.661	0.340	1.9e-21
UNIPROTKB|F1NJ02	398	TRDMT1 "Uncharacterized protei	0.529	0.311	0.442	1.9e-21
UNIPROTKB|O14717	391	TRDMT1 "tRNA (cytosine(38)-C(5	0.876	0.524	0.340	2.8e-21
UNIPROTKB|Q7YS61	391	TRDMT1 "tRNA (cytosine-5-)-met	0.470	0.281	0.475	3.7e-21
UNIPROTKB|E2RM86	391	TRDMT1 "Uncharacterized protei	0.871	0.521	0.338	4.8e-21
MGI|MGI:1274787	415	Trdmt1 "tRNA aspartic acid met	0.470	0.265	0.449	1.2e-20
ZFIN|ZDB-GENE-041008-138	381	trdmt1 "tRNA aspartic acid met	0.5	0.307	0.407	4.5e-20
RGD|1306292	391	Trdmt1 "tRNA aspartic acid met	0.474	0.283	0.443	5.3e-20
FB|FBgn0028707	345	Mt2 "Methyltransferase 2" [Dro	0.235	0.159	0.545	8.2e-17

TAIR|locus:2179469 DNMT2 "DNA methyltransferase-2" [Arabidopsis thaliana (taxid:3702)]

 Score = 659 (237.0 bits), Expect = 1.1e-64, P = 1.1e-64
 Identities = 137/212 (64%), Positives = 156/212 (73%)

Query:    19 AKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFS 78
             AKRKPL F+ Q  NN+LL SP PL G DD     K  Q ++  DKLLE C PVE+FLE +
Sbjct:   177 AKRKPLPFKSQHSNNKLLWSPDPLYGRDDQVEFGKC-QAEEGLDKLLEFCKPVEKFLELA 235

Query:    79 NSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYP 138
                D   +       G+   D  G   ++V  D   S+  +LVP+SLIERWG+AMDIVYP
Sbjct:   236 AHVDGEPSSVDDSENGSK--DCCGQEGDSVP-D---SVHQYLVPVSLIERWGNAMDIVYP 289

Query:   139 DSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDF 198
             DSKRCCCFTKSYYRYVKGTGSLLATVQPK KGK S LKEQ LRYFTPREVAN HSFP DF
Sbjct:   290 DSKRCCCFTKSYYRYVKGTGSLLATVQPKIKGKESCLKEQRLRYFTPREVANFHSFPEDF 349

Query:   199 QFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
             +FP H+SLRQRYA+LGNSLS+AVVAPLL+YLF
Sbjct:   350 EFPKHISLRQRYAMLGNSLSVAVVAPLLRYLF 381

GO:0003677 "DNA binding" evidence=IEA;ISS

GO:0005634 "nucleus" evidence=ISM

GO:0006306 "DNA methylation" evidence=IEA;ISS

UNIPROTKB|E7EMI8 TRDMT1 "tRNA (cytosine(38)-C(5))-methyltransferase" [Homo sapiens (taxid:9606)]

 Score = 251 (93.4 bits), Expect = 1.9e-21, P = 1.9e-21
 Identities = 80/235 (34%), Positives = 126/235 (53%)

Query:    16 LGKAKRKPLSFRCQ---LLNNQLLRSPSPLLGNDDMT--VITKHDQPDDSWDKLLESC-- 68
             + K + +PL F+     L+    + S  P     D+   +  K+ +P+ S+D  ++    
Sbjct:    85 IAKLQSEPLPFQAPGQVLMEFPKIESVHPQKYAMDVENKIQEKNVEPNISFDGSIQCSGK 144

Query:    69 DPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIER 128
             D +   LE +    + N +   LS     + DF   E+  +V++     + L P SL+ R
Sbjct:   145 DAILFKLETAEEIHRKNQQDSDLSV--KMLKDF--LEDDTDVNQ-----YLLPPKSLL-R 194

Query:   129 WGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKNKGKAS---SLKEQ---- 178
             +   +DIV P  +R  CFTK Y  Y++GTGS+L T   VQ +N  K+    S +EQ    
Sbjct:   195 YALLLDIVQPTCRRSVCFTKGYGSYIEGTGSVLQTAEDVQVENIYKSLTNLSQEEQITKL 254

Query:   179 ---HLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
                 LRYFTP+E+ANL  FP +F FP  ++++QRY LLGNSL++ VVA L++ L+
Sbjct:   255 LILKLRYFTPKEIANLLGFPPEFGFPEKITVKQRYRLLGNSLNVHVVAKLIKILY 309

GO:0003677 "DNA binding" evidence=IEA

GO:0006306 "DNA methylation" evidence=IEA

GO:0008168 "methyltransferase activity" evidence=IEA

UNIPROTKB|F1NJ02 TRDMT1 "Uncharacterized protein" [Gallus gallus (taxid:9031)]

 Score = 253 (94.1 bits), Expect = 1.9e-21, P = 1.9e-21
 Identities = 61/138 (44%), Positives = 84/138 (60%)

Query:   105 EETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT- 163
             E+ +E +      +FL P SL+ R+   +DIV P  +R  CFTK Y  YV+GTGS+L T 
Sbjct:   260 EDFLEEENEEMSQYFLAPKSLV-RYAFLLDIVKPTCRRSTCFTKGYGHYVEGTGSVLQTA 318

Query:   164 --VQPKNKGK-ASSLKEQH---------LRYFTPREVANLHSFPGDFQFPHHLSLRQRYA 211
               V+ ++  K   +L E+          LRYFTPRE+ANLH FP +F FP  +S++Q Y 
Sbjct:   319 VDVELESVFKHIENLTEEEKLMKLSTLKLRYFTPREIANLHGFPPEFGFPDKVSIKQCYR 378

Query:   212 LLGNSLSIAVVAPLLQYL 229
             LLGNSL++ VVA L+  L
Sbjct:   379 LLGNSLNVRVVAKLISVL 396

GO:0003677 "DNA binding" evidence=IEA

GO:0006306 "DNA methylation" evidence=IEA

GO:0008168 "methyltransferase activity" evidence=IEA

UNIPROTKB|O14717 TRDMT1 "tRNA (cytosine(38)-C(5))-methyltransferase" [Homo sapiens (taxid:9606)]

 Score = 251 (93.4 bits), Expect = 2.8e-21, P = 2.8e-21
 Identities = 80/235 (34%), Positives = 126/235 (53%)

Query:    16 LGKAKRKPLSFRCQ---LLNNQLLRSPSPLLGNDDMT--VITKHDQPDDSWDKLLESC-- 68
             + K + +PL F+     L+    + S  P     D+   +  K+ +P+ S+D  ++    
Sbjct:   166 IAKLQSEPLPFQAPGQVLMEFPKIESVHPQKYAMDVENKIQEKNVEPNISFDGSIQCSGK 225

Query:    69 DPVERFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIER 128
             D +   LE +    + N +   LS     + DF   E+  +V++     + L P SL+ R
Sbjct:   226 DAILFKLETAEEIHRKNQQDSDLSV--KMLKDF--LEDDTDVNQ-----YLLPPKSLL-R 275

Query:   129 WGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKNKGKAS---SLKEQ---- 178
             +   +DIV P  +R  CFTK Y  Y++GTGS+L T   VQ +N  K+    S +EQ    
Sbjct:   276 YALLLDIVQPTCRRSVCFTKGYGSYIEGTGSVLQTAEDVQVENIYKSLTNLSQEEQITKL 335

Query:   179 ---HLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
                 LRYFTP+E+ANL  FP +F FP  ++++QRY LLGNSL++ VVA L++ L+
Sbjct:   336 LILKLRYFTPKEIANLLGFPPEFGFPEKITVKQRYRLLGNSLNVHVVAKLIKILY 390

GO:0003677 "DNA binding" evidence=IEA

GO:0003723 "RNA binding" evidence=IEA

GO:0008033 "tRNA processing" evidence=IEA

GO:0001975 "response to amphetamine" evidence=IEA

GO:0005634 "nucleus" evidence=IEA

GO:0003886 "DNA (cytosine-5-)-methyltransferase activity" evidence=TAS

GO:0090116 "C-5 methylation of cytosine" evidence=TAS

UNIPROTKB|Q7YS61 TRDMT1 "tRNA (cytosine-5-)-methyltransferase" [Bos taurus (taxid:9913)]

 Score = 250 (93.1 bits), Expect = 3.7e-21, P = 3.7e-21
 Identities = 59/124 (47%), Positives = 82/124 (66%)

Query:   119 FLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKNKGKA-SS 174
             FL P SL+ R+   +DIV P S+R  CFTK Y RY++GTGS+L T   VQ +N  K+ +S
Sbjct:   267 FLPPKSLL-RYALLLDIVKPTSRRSMCFTKGYGRYIEGTGSVLQTTEDVQIENIYKSLTS 325

Query:   175 LKEQH---------LRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPL 225
             L ++          LR+FTP+E+ANL  FP +F FP   +++QRY LLGNSL++ VVA L
Sbjct:   326 LSQEEKIMRLSMLQLRFFTPKEIANLLGFPPEFGFPEMTTVKQRYRLLGNSLNVHVVAKL 385

Query:   226 LQYL 229
             ++ L
Sbjct:   386 IKIL 389

GO:0005634 "nucleus" evidence=IEA

GO:0008168 "methyltransferase activity" evidence=IEA

GO:0008033 "tRNA processing" evidence=IEA

GO:0003723 "RNA binding" evidence=IEA

GO:0006306 "DNA methylation" evidence=IEA

GO:0003677 "DNA binding" evidence=IEA

UNIPROTKB|E2RM86 TRDMT1 "Uncharacterized protein" [Canis lupus familiaris (taxid:9615)]

 Score = 249 (92.7 bits), Expect = 4.8e-21, P = 4.8e-21
 Identities = 80/236 (33%), Positives = 126/236 (53%)

Query:    16 LGKAKRKPLSFRCQ---LLNNQLLRSPSPLLGNDDMT--VITKHDQPDDSWDKLLESCDP 70
             + K + +PL F+     L+      S +P     D    +  K  +P+  +DK  + C  
Sbjct:   166 IAKLQSEPLPFQAPGQVLMEFPKTESENPEKNTIDAENKIERKKIEPNICFDKGRQ-CSG 224

Query:    71 VERFL---EFSNSGDQVNTETGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIE 127
              E  L   E     D+ + +   LS     + DF   E+ V++++     +FL P SL+ 
Sbjct:   225 KETILFKLETIQELDRKHHQDNDLSV--QMLKDF--LEDDVDMNQ-----YFLPPKSLL- 274

Query:   128 RWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKNKGKA-SSLKEQH---- 179
             R+   +DIV P  +R  CFTK Y  YV+GTGS+L T   VQ ++  K+ ++L ++     
Sbjct:   275 RYALLLDIVKPTCRRSTCFTKGYGSYVEGTGSVLQTTEDVQIEDIYKSLTNLSQEEKITK 334

Query:   180 -----LRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
                  LRYFTP+E+ANL  FP +F FP  ++++QRY LLGNSL++ VVA L++ L+
Sbjct:   335 LLMLKLRYFTPKEIANLLGFPSEFGFPEKITVKQRYRLLGNSLNVHVVAKLIKILY 390

GO:0008168 "methyltransferase activity" evidence=IEA

GO:0006306 "DNA methylation" evidence=IEA

GO:0003677 "DNA binding" evidence=IEA

MGI|MGI:1274787 Trdmt1 "tRNA aspartic acid methyltransferase 1" [Mus musculus (taxid:10090)]

 Score = 247 (92.0 bits), Expect = 1.2e-20, P = 1.2e-20
 Identities = 58/129 (44%), Positives = 77/129 (59%)

Query:   117 DHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQ----------- 165
             D +L+P  L+ R+   +DIV P S+R  CFTK Y  Y++GTGS+L   +           
Sbjct:   264 DEYLLPPKLLLRYALLLDIVKPTSRRSMCFTKGYGSYIEGTGSVLQAAEDAQIENIYKSL 323

Query:   166 ----PKNK-GKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIA 220
                 P+ K  K S LK   LRYFTP+E+ANL  FP +F FP   +++QRY LLGNSL++ 
Sbjct:   324 PDLPPEEKIAKLSMLK---LRYFTPKEIANLQGFPPEFGFPEKTTVKQRYRLLGNSLNVH 380

Query:   221 VVAPLLQYL 229
             VVA LL  L
Sbjct:   381 VVAKLLTVL 389

GO:0003677 "DNA binding" evidence=IEA

GO:0003723 "RNA binding" evidence=IEA

GO:0005634 "nucleus" evidence=IEA

GO:0006306 "DNA methylation" evidence=IEA

GO:0008033 "tRNA processing" evidence=IEA

GO:0008168 "methyltransferase activity" evidence=IEA

GO:0016740 "transferase activity" evidence=IEA

GO:0032259 "methylation" evidence=IEA

ZFIN|ZDB-GENE-041008-138 trdmt1 "tRNA aspartic acid methyltransferase 1" [Danio rerio (taxid:7955)]

 Score = 240 (89.5 bits), Expect = 4.5e-20, P = 4.5e-20
 Identities = 53/130 (40%), Positives = 77/130 (59%)

Query:   116 IDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLA--------TV--- 164
             ++ +L+P   + R+   MDIV P S+R  CFTK Y  YV+GTGS+L         TV   
Sbjct:   252 MEPYLLPPKTLLRYALVMDIVQPSSRRSVCFTKGYGHYVEGTGSVLQSCMDVDLETVFKN 311

Query:   165 --QPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVV 222
               Q   + K   L    LRYFTPRE++ L  FP  F FP H+S +Q+Y +LGNSL++ VV
Sbjct:   312 LDQLSEEDKLKQLLRLKLRYFTPREISRLMGFPDHFTFPKHISFKQQYRVLGNSLNVHVV 371

Query:   223 APLLQYLFAQ 232
             + L++ + ++
Sbjct:   372 SHLIRLMLSK 381

GO:0006306 "DNA methylation" evidence=IEA

GO:0003677 "DNA binding" evidence=IEA

GO:0008168 "methyltransferase activity" evidence=IEA

GO:0005575 "cellular_component" evidence=ND

GO:0016740 "transferase activity" evidence=IEA

GO:0032259 "methylation" evidence=IEA

RGD|1306292 Trdmt1 "tRNA aspartic acid methyltransferase 1" [Rattus norvegicus (taxid:10116)]

 Score = 240 (89.5 bits), Expect = 5.3e-20, P = 5.3e-20
 Identities = 55/124 (44%), Positives = 76/124 (61%)

Query:   119 FLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKN------- 168
             +L+P   + R+   +DIV P S+R  CFTK Y  Y++GTGS+L T   VQ +N       
Sbjct:   266 YLLPAKSLLRYALLLDIVKPTSRRSMCFTKGYGSYIEGTGSVLQTAEDVQIENIYKSLPD 325

Query:   169 ---KGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPL 225
                + K + L    LRYFTP+E+ANL  FP +F FP   +++QRY LLGNSL++ VV+ L
Sbjct:   326 LPPEEKIAKLSMLKLRYFTPKEIANLLGFPPEFGFPEKTTVKQRYRLLGNSLNVHVVSKL 385

Query:   226 LQYL 229
             L  L
Sbjct:   386 LTVL 389

GO:0001975 "response to amphetamine" evidence=IEP

GO:0003674 "molecular_function" evidence=ND

GO:0003677 "DNA binding" evidence=IEA

GO:0003723 "RNA binding" evidence=IEA

GO:0005575 "cellular_component" evidence=ND

GO:0005634 "nucleus" evidence=IEA

GO:0006306 "DNA methylation" evidence=IEA

GO:0008033 "tRNA processing" evidence=IEA

GO:0008168 "methyltransferase activity" evidence=IEA

FB|FBgn0028707 Mt2 "Methyltransferase 2" [Drosophila melanogaster (taxid:7227)]

 Score = 153 (58.9 bits), Expect = 8.2e-17, Sum P(2) = 8.2e-17
 Identities = 30/55 (54%), Positives = 40/55 (72%)

Query:   175 LKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYL 229
             L +  LRYFTPREVA L SFP +F+FP   + RQ+Y LLGNS+++ VV  L++ L
Sbjct:   288 LHQVRLRYFTPREVARLMSFPENFEFPPETTNRQKYRLLGNSINVKVVGELIKLL 342

GO:0003886 "DNA (cytosine-5-)-methyltransferase activity" evidence=ISS;IDA;NAS

GO:0008168 "methyltransferase activity" evidence=ISS

GO:0006306 "DNA methylation" evidence=ISS;IDA;NAS

GO:0000790 "nuclear chromatin" evidence=NAS

GO:0043045 "DNA methylation involved in embryo development" evidence=IMP

GO:0005634 "nucleus" evidence=IDA

GO:0003677 "DNA binding" evidence=IEA

GO:0008340 "determination of adult lifespan" evidence=IMP

GO:0030488 "tRNA methylation" evidence=IMP

GO:0016458 "gene silencing" evidence=IMP

GO:0000723 "telomere maintenance" evidence=IMP

GO:0044026 "DNA hypermethylation" evidence=IMP

GO:0006979 "response to oxidative stress" evidence=IMP

GO:0010494 "cytoplasmic stress granule" evidence=IDA

GO:0008175 "tRNA methyltransferase activity" evidence=IMP

GO:0009408 "response to heat" evidence=IMP

---

Prediction of Enzyme Commission (EC) Number

EC Number Prediction by Annotation Transfer from SWISS-PROT Entries

Original result of BLAST against SWISS-PROT

No confident hit for EC number transfering in SWISSPROT detected by BLAST

EC Number Prediction by Ezypred Server

Original result from Ezypred Server

Fail to connect to Ezypred Server

EC Number Prediction by EFICAz Software

No EC number assignment, probably not an enzyme!

---

Prediction of Functionally Associated Proteins

Functionally Associated Proteins Detected by STRING

Original result from the STRING server

Your Input:
	GSVIVG00032604001	SubName- Full=Chromosome chr4 scaffold_6, whole genome shotgun sequence; (377 aa)
		(Vitis vinifera)
Predicted Functional Partners:
	GSVIVG00035341001	RecName- Full=Cytosine-specific methyltransferase; EC=2.1.1.37; (444 aa)			•				•		0.483

---

Conserved Domains and Related Protein Families

Conserved Domains Detected by RPS-BLAST

Original result of RPS-BLAST against CDD database part I

ID	Length	Definition	E-value
Query	234
cd00315	275	cd00315, Cyt_C5_DNA_methylase, Cytosine-C5 specifi	8e-19
TIGR00675	315	TIGR00675, dcm, DNA-methyltransferase (dcm)	3e-05
pfam00145	320	pfam00145, DNA_methylase, C-5 cytosine-specific DN	7e-05
COG0270	328	COG0270, Dcm, Site-specific DNA methylase [DNA rep	1e-04

>gnl|CDD|238192 cd00315, Cyt_C5_DNA_methylase, Cytosine-C5 specific DNA methylases; Methyl transfer reactions play an important role in many aspects of biology

 Score = 82.3 bits (204), Expect = 8e-19
 Identities = 36/116 (31%), Positives = 49/116 (42%), Gaps = 11/116 (9%)

Query: 115 SIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASS 174
               F  P    +     + I  PD       T SY    KGTGS+  T           
Sbjct: 170 FFSPFPKPSEKKKTLKDILRIRDPDE-PSPTLTASYG---KGTGSVHPTAPDMI------ 219

Query: 175 LKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
            KE ++R  TPRE A L  FP DF+FP   S+ Q Y  +GNS+ + V   + + + 
Sbjct: 220 GKESNIRRLTPRECARLQGFPDDFEFPGK-SVTQAYRQIGNSVPVPVAEAIAKAIK 274

Cytosine-specific DNA methylases are found both in prokaryotes and eukaryotes. DNA methylation, or the covalent addition of a methyl group to cytosine within the context of the CpG dinucleotide, has profound effects on the mammalian genome. These effects include transcriptional repression via inhibition of transcription factor binding or the recruitment of methyl-binding proteins and their associated chromatin remodeling factors, X chromosome inactivation, imprinting and the suppression of parasitic DNA sequences. DNA methylation is also essential for proper embryonic development and is an important player in both DNA repair and genome stability. Length = 275

>gnl|CDD|233085 TIGR00675, dcm, DNA-methyltransferase (dcm)

 Score = 43.8 bits (104), Expect = 3e-05
 Identities = 19/42 (45%), Positives = 25/42 (59%), Gaps = 2/42 (4%)

Query: 181 RYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVV 222
           R  TPRE A L  FP DF+FP  +S  Q Y   GN++ + V+
Sbjct: 269 RRLTPRECARLQGFPDDFKFP--VSDSQLYKQAGNAVVVPVI 308

All proteins in this family for which functions are known are DNA-cytosine methyltransferases. This family is based on the phylogenomic analysis of JA Eisen (1999, Ph.D. Thesis, Stanford University) [DNA metabolism, DNA replication, recombination, and repair]. Length = 315

>gnl|CDD|215747 pfam00145, DNA_methylase, C-5 cytosine-specific DNA methylase

 Score = 42.7 bits (101), Expect = 7e-05
 Identities = 25/88 (28%), Positives = 31/88 (35%), Gaps = 8/88 (9%)

Query: 145 CFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHL 204
               S+     G G    T     K       +   R  TPRE A L  FP DF F    
Sbjct: 238 GKGPSFTYRKSGRGEAPKTG-ILGKNGERFRPKNI-RRLTPRECARLQGFPDDFIF--PG 293

Query: 205 SLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
           S  Q Y  +GN    AV  P+ + +   
Sbjct: 294 SKTQLYKQIGN----AVPVPVAEAIAKA 317

Length = 320

>gnl|CDD|223348 COG0270, Dcm, Site-specific DNA methylase [DNA replication, recombination, and repair]

 Score = 42.0 bits (99), Expect = 1e-04
 Identities = 17/43 (39%), Positives = 21/43 (48%), Gaps = 2/43 (4%)

Query: 181 RYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVA 223
           R  T RE A L  FP DF FP   S   +Y  +GNS+   +  
Sbjct: 274 RELTVREAARLQGFPDDFVFP--GSKTDQYRQIGNSVPPLLAE 314

Length = 328

Conserved Domains Detected by HHsearch

Original result of HHsearch against CDD database

ID	Length	Definition	Probability
Query	234
KOG0919	338	consensus C-5 cytosine-specific DNA methylase [Tra	100.0
PRK10458	467	DNA cytosine methylase; Provisional	99.87
TIGR00675	315	dcm DNA-methyltransferase (dcm). All proteins in t	99.7
PF00145	335	DNA_methylase: C-5 cytosine-specific DNA methylase	99.66
cd00315	275	Cyt_C5_DNA_methylase Cytosine-C5 specific DNA meth	99.65
COG0270	328	Dcm Site-specific DNA methylase [DNA replication,	99.43

>KOG0919 consensus C-5 cytosine-specific DNA methylase [Transcription]

Probab=100.00  E-value=1.2e-35  Score=266.66  Aligned_cols=169  Identities=34%  Similarity=0.523  Sum_probs=135.3

Q ss_pred             CCCCCCceeEeeecCCCCccccccccccccCCCCCCCCCcccccccCCCCCcchhhcccCcchHHHhhcccCCCCCcccc
Q 026734            8 FPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTE   87 (234)
Q Consensus         8 ~~~~~~ry~~lAk~~~~~f~~~~~n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~i~~~le~~~~~~~~~~~   87 (234)
                      .|-+|+||||+||.. .+|+++...   .|..+|.+     ..         ..+ -.+.++-|.+|||.+         
T Consensus       156 iPNsR~Ryy~iArl~-~~F~~~G~~---s~d~~~qF-----se---------iaq-k~g~Vk~i~d~lE~~---------  207 (338)
T KOG0919|consen  156 IPNSRYRYYCIARLG-ADFPFAGGK---SWDEMPQF-----SE---------IAQ-KQGLVKQIADILEEN---------  207 (338)
T ss_pred             CCCcchheeehhhhC-CCCCCCCCc---ccccccch-----HH---------HHH-hcchHHHHHHHHHhc---------
Confidence            599999999999976 467665533   45544443     01         111 133467899999952         


Q ss_pred             ccccccCCccccccCCccccchhcccccccccCCChhHHhhhcCcceeeccCCCCCcccccceeeeecCCcceeeecCC-
Q 026734           88 TGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQP-  166 (234)
Q Consensus        88 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~y~v~~~~l~k~g~~~DIV~p~s~~s~cFTksY~r~~~GtgSvl~~~~~-  166 (234)
                                                ..+.+|+||+++|.||+..||||.|.+.+|+||||+|++|.+|+||++....- 
T Consensus       208 --------------------------~d~s~ylvp~~vL~k~~l~~DIv~P~~srs~CFTkGYthy~eGtGSilq~~~~i  261 (338)
T KOG0919|consen  208 --------------------------VDPSDYLVPDDVLTKRVLVMDIVHPAQSRSMCFTKGYTHYTEGTGSILQLVKEI  261 (338)
T ss_pred             --------------------------CCHHHccCCHHHHHHhHhheeecccccccceEeecCccceeecchHHHHHHhhh
Confidence                                      12467999999999999999999999999999999999999999999865431 


Q ss_pred             --CC----------CCCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHH
Q 026734          167 --KN----------KGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLF  230 (234)
Q Consensus       167 --~~----------~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll  230 (234)
                        .+          ..+...++..++|||||||+|||||||.+|.||..++.+|+|+++|||++|.||..|++-|.
T Consensus       262 ~~eN~~~s~~~~~~~~~~~~l~~l~LRYFTprEvArLmgFPe~fefp~~~T~kq~YRLLGNSiNVkVV~~LIklL~  337 (338)
T KOG0919|consen  262 DTENQDASKSEKILQQRLDLLHQLRLRYFTPREVARLMGFPENFEFPPETTNKQKYRLLGNSINVKVVGELIKLLT  337 (338)
T ss_pred             cccchhHHHHHHHHHHHHHHHHHHHhhccCHHHHHHHcCCCcccCCCcchhHHHHHHHhcCcccceeHHHHHHHhc
Confidence              00          01124567789999999999999999999999988999999999999999999999998875

>PRK10458 DNA cytosine methylase; Provisional

Probab=99.87  E-value=1.6e-21  Score=187.88  Aligned_cols=114  Identities=29%  Similarity=0.443  Sum_probs=84.2

Q ss_pred             cccCCChhHHh---hh---------cCcceeeccCCCCCcccccceeeeecCCcceeeecC-CCCCCCCcccc----CCC
Q 026734          117 DHFLVPLSLIE---RW---------GSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQ-PKNKGKASSLK----EQH  179 (234)
Q Consensus       117 ~~y~v~~~~l~---k~---------g~~~DIV~p~s~~s~cFTksY~r~~~GtgSvl~~~~-~~~~~k~~~l~----~~~  179 (234)
                      ++|++++++..   +|         |+.++++.|.+.++.|+|.+| +|.+++..+++... ....+...+.+    ..+
T Consensus       322 ~ky~ls~~~~~~l~~~~~k~~~~g~g~~~~i~~~~~~~~~~~t~~~-ry~k~gs~~~i~~~~~~~~~~~~~~~~~~~~~~  400 (467)
T PRK10458        322 AKYILTPVLWKYLYRYAKKHQAKGNGFGYGLVYPNNPQSVTRTLSA-RYYKDGSEILIDRGWDMALGEKDFDDPENQQHR  400 (467)
T ss_pred             cceeeCHHHHHHHHHHHhhccccCCCcceeeeecCCCCCccccccc-ccccCCCceeeecccccccccccccccccccCC
Confidence            45777766543   33         567899999888899999994 88888544444321 01111112222    356


Q ss_pred             cccCcHHHHHHhCCC--CCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          180 LRYFTPREVANLHSF--PGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       180 lR~lTprE~aRLQgF--Pddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      +|+||||||||||||  |++|.|+..++.+|+|+|+||||+|+|+++|++.|..
T Consensus       401 ~RrLTprE~aRLqGF~~pd~~~F~~~vSdtq~Ykq~GNSV~Vpvv~aIa~~L~~  454 (467)
T PRK10458        401 PRRLTPRECARLMGFEAPGEAKFRIPVSDTQAYRQFGNSVVVPVFAAVAKLLEP  454 (467)
T ss_pred             cccCCHHHHHHhCCCCCCccccccCCCCHHHHHHHhCCcccHHHHHHHHHHHHH
Confidence            899999999999999  6667776558999999999999999999999999864

>TIGR00675 dcm DNA-methyltransferase (dcm)

Probab=99.70  E-value=8.8e-17  Score=146.89  Aligned_cols=53  Identities=36%  Similarity=0.692  Sum_probs=50.4

Q ss_pred             ccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHH
Q 026734          175 LKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYL  229 (234)
Q Consensus       175 l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~L  229 (234)
                      +|+.+.|.||+||+||||||||+|.|+  ++.+++|+||||||||+++++|++.|
T Consensus       263 ~hp~~~R~lT~RE~aRLQ~FPd~f~f~--~s~~~~~~qiGNAVPp~la~~I~~~i  315 (315)
T TIGR00675       263 VHPGRIRRLTPRECARLQGFPDDFKFP--VSDSQLYKQAGNAVVVPVIEAIAKQI  315 (315)
T ss_pred             ccCCceeeCCHHHHHHHcCCCcccEeC--CCHHHHHhhhCCcccHHHHHHHHhhC
Confidence            788999999999999999999999997  59999999999999999999999875

All proteins in this family for which functions are known are DNA-cytosine methyltransferases. This family is based on the phylogenomic analysis of JA Eisen (1999, Ph.D. Thesis, Stanford University).

>PF00145 DNA_methylase: C-5 cytosine-specific DNA methylase; InterPro: IPR001525 C-5 cytosine-specific DNA methylases (2

Probab=99.66  E-value=8e-17  Score=141.54  Aligned_cols=56  Identities=34%  Similarity=0.600  Sum_probs=47.6

Q ss_pred             cccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          174 SLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       174 ~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      .+|+.+.|+|||||+||||||||+|.|+.  +..++|+||||||||+|+++|++.|++
T Consensus       280 ~~hp~~~R~LT~rE~aRLqgFPd~~~f~g--~~~~~~~qiGNAVpp~v~~~I~~~i~~  335 (335)
T PF00145_consen  280 FIHPEQNRRLTPREAARLQGFPDDFKFPG--SKTQQYKQIGNAVPPPVAEAIAKAIKK  335 (335)
T ss_dssp             EBTTSSSCB-BHHHHHHHTTSSTTS-S-S--SHHHHHHHHHCS--HHHHHHHHHHHH-
T ss_pred             ccCCCCCCcCcHHHHHHhCCCCCceEccC--CHHHHhceECCCcCHHHHHHHHHHhhC
Confidence            68999999999999999999999999984  888999999999999999999999874

1.1.37 from EC) (C5 Mtase) are enzymes that specifically methylate the C-5 carbon of cytosines in DNA to produce C5-methylcytosine [, , ]. In mammalian cells, cytosine-specific methyltransferases methylate certain CpG sequences, which are believed to modulate gene expression and cell differentiation. In bacteria, these enzymes are a component of restriction-modification systems and serve as valuable tools for the manipulation of DNA [, ]. The structure of HhaI methyltransferase (M.HhaI) has been resolved to 2.5 A []: the molecule folds into 2 domains - a larger catalytic domain containing catalytic and cofactor binding sites, and a smaller DNA recognition domain.; GO: 0003677 DNA binding, 0006306 DNA methylation; PDB: 4DA4_A 3PT6_B 3AV6_A 3AV5_A 3AV4_A 3PT9_A 1DCT_A 3LX6_A 3ME5_A 2QRV_A ....

>cd00315 Cyt_C5_DNA_methylase Cytosine-C5 specific DNA methylases; Methyl transfer reactions play an important role in many aspects of biology

Probab=99.65  E-value=2.1e-16  Score=141.87  Aligned_cols=86  Identities=38%  Similarity=0.597  Sum_probs=66.9

Q ss_pred             eeccCCCCCcccccceeeeecCCcceeeecCCCCCCCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhc
Q 026734          135 IVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLG  214 (234)
Q Consensus       135 IV~p~s~~s~cFTksY~r~~~GtgSvl~~~~~~~~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiG  214 (234)
                      ++.+++ .+.|+|..|++   +++++.....      ....++.+.|+||+||+||||||||+|.|+.. +.+++|+|||
T Consensus       190 ~~~~~~-~~~ti~~~~~~---~~~~~~~~~~------~~~~~~~~~R~lT~rE~arlqgFPd~f~f~g~-~~~~~~~qiG  258 (275)
T cd00315         190 IRDPDE-PSPTLTASYGK---GTGSVHPTAP------DMIGKESNIRRLTPRECARLQGFPDDFEFPGK-SVTQAYRQIG  258 (275)
T ss_pred             hhcCCC-CccceecCCCC---CccccccCcc------cccccCCCCCCCCHHHHHHHcCCCCCcEEcCC-CHHHHHHhhc
Confidence            344443 56788877665   4555543211      11457788999999999999999999999732 8999999999


Q ss_pred             CcccHHHHHHHHHHHHh
Q 026734          215 NSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       215 NAV~v~vv~~La~~Ll~  231 (234)
                      |||||+++++|++.|.+
T Consensus       259 NAVp~~~~~~I~~~i~~  275 (275)
T cd00315         259 NSVPVPVAEAIAKAIKE  275 (275)
T ss_pred             CCcCHHHHHHHHHHHhC
Confidence            99999999999999863

Cytosine-specific DNA methylases are found both in prokaryotes and eukaryotes. DNA methylation, or the covalent addition of a methyl group to cytosine within the context of the CpG dinucleotide, has profound effects on the mammalian genome. These effects include transcriptional repression via inhibition of transcription factor binding or the recruitment of methyl-binding proteins and their associated chromatin remodeling factors, X chromosome inactivation, imprinting and the suppression of parasitic DNA sequences. DNA methylation is also essential for proper embryonic development and is an important player in both DNA repair and genome stability.

>COG0270 Dcm Site-specific DNA methylase [DNA replication, recombination, and repair]

Probab=99.43  E-value=2.5e-13  Score=124.73  Aligned_cols=79  Identities=27%  Similarity=0.393  Sum_probs=66.4

Q ss_pred             ccceeeeecCCcceeeecCCCCCCCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHH
Q 026734          147 TKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLL  226 (234)
Q Consensus       147 TksY~r~~~GtgSvl~~~~~~~~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La  226 (234)
                      ...|.++.+...+..+..    ......+|+...|.||++|++|||||||+|.|..  +..++|+||||||||+++++|+
T Consensus       244 ~~~~~rl~~~~~~~t~~~----~~~~~~~h~~~~r~lt~rE~arlq~fPd~~~~~g--s~~~~~~qiGnsVp~~l~~~ia  317 (328)
T COG0270         244 KGSYIRLHPDKPAPTVRG----GGNERFIHPLEDRELTVREAARLQGFPDDFVFPG--SKTDQYRQIGNSVPPLLAEAIA  317 (328)
T ss_pred             CceeEeCCCCCCCceeec----CCCcccCCCCcCCCCCHHHHHHhcCCCCceEEec--cchhhhhhccCcCCHHHHHHHH
Confidence            456667777666655541    2335678999999999999999999999999985  9999999999999999999999


Q ss_pred             HHHHh
Q 026734          227 QYLFA  231 (234)
Q Consensus       227 ~~Ll~  231 (234)
                      +.|++
T Consensus       318 ~~i~~  322 (328)
T COG0270         318 KAILK  322 (328)
T ss_pred             HHHHH
Confidence            99975

---

Homologous Structure Templates

Structure Templates Detected by BLAST

Original result of BLAST against Protein Data Bank

ID	Length	Definition	E-value
Query	234
1g55_A	343	Structure Of Human Dnmt2, An Enigmatic Dna Methyltr		3e-23
3qv2_A	327	Structure Analysis Of Entamoeba Histolytica Methylt		4e-20
4h0n_A	333	Crystal Structure Of Spodoptera Frugiperda Dnmt2 E2		6e-20

>pdb|1G55|A Chain A, Structure Of Human Dnmt2, An Enigmatic Dna Methyltransferase Homologue Length = 343

Iteration: 1
 Score =  105 bits (261), Expect = 3e-23,   Method: Compositional matrix adjust.
 Identities = 56/128 (43%), Positives = 82/128 (64%), Gaps = 13/128 (10%)
Query: 116 IDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLAT---VQPKNKGKA 172
           ++ +L+P   + R+   +DIV P  +R  CFTK Y  Y++GTGS+L T   VQ +N  K+
Sbjct: 215 VNQYLLPPKSLLRYALLLDIVQPTXRRSVCFTKGYGSYIEGTGSVLQTAEDVQVENIYKS 274
Query: 173 S---SLKEQ-------HLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVV 222
               S +EQ        LRYFTP+E+ANL  FP +F FP  ++++QRY LLGNSL++ VV
Sbjct: 275 LTNLSQEEQITKLLILKLRYFTPKEIANLLGFPPEFGFPEKITVKQRYRLLGNSLNVHVV 334
Query: 223 APLLQYLF 230
           A L++ L+
Sbjct: 335 AKLIKILY 342

>pdb|3QV2|A Chain A, Structure Analysis Of Entamoeba Histolytica Methyltransferase Ehmeth Length = 327

Iteration: 1
 Score = 94.7 bits (234), Expect = 4e-20,   Method: Compositional matrix adjust.
 Identities = 50/117 (42%), Positives = 71/117 (60%), Gaps = 3/117 (2%)
Query: 117 DHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPK--NKGKASS 174
           + + +P  LI + G   DIV  D KR CCFTKSY + V+GTGS+   ++P      KA  
Sbjct: 206 ESYSIPSDLILKKGMLFDIVGKDDKRTCCFTKSYTKIVEGTGSIYCPIEPHFIPVKKAED 265
Query: 175 LKEQHLRYFTPREVANLHSFPGDFQFP-HHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
           L  ++LRYFTP E+  +H F  +F      L+ +Q+Y  LGNS+S  V+A L++YLF
Sbjct: 266 LLNKNLRYFTPNEIKKIHGFSSNFTTQIDGLTDKQQYQCLGNSVSCFVIAQLMEYLF 322

>pdb|4H0N|A Chain A, Crystal Structure Of Spodoptera Frugiperda Dnmt2 E260a/e261a/k263a Mutant Length = 333

Iteration: 1
 Score = 94.0 bits (232), Expect = 6e-20,   Method: Compositional matrix adjust.
 Identities = 57/149 (38%), Positives = 79/149 (53%), Gaps = 22/149 (14%)
Query: 100 DFG---AAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKG 156
           DFG   + E  +E D     + FLVP  ++ R     DI Y  SKR CCFTK+Y  Y  G
Sbjct: 188 DFGVPHSLESIIEED---VDEKFLVPEKML-RCAKVFDICYKTSKRSCCFTKAYTHYADG 243
Query: 157 TGSLL-------------ATVQPKNKGK--ASSLKEQHLRYFTPREVANLHSFPGDFQFP 201
           TGS+              A  Q +  G+      KE  LRYFTP+EV  +  FP  +  P
Sbjct: 244 TGSIFTDKPREVVQKCYAAAAQNEIGGEKFVELFKELKLRYFTPKEVLMIMCFPKSYNLP 303
Query: 202 HHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
            ++S++Q Y LLGNS+++ V++ LL+ LF
Sbjct: 304 TNISMKQCYRLLGNSVNVKVISELLKILF 332

Structure Templates Detected by RPS-BLAST

Original result of RPS-BLAST against PDB70 database

ID	Length	Definition	E-value
Query	234
1g55_A	343	DNA cytosine methyltransferase DNMT2; human DNA me	1e-23
3qv2_A	327	5-cytosine DNA methyltransferase; DNMT2, ehmeth; H	3e-22
1dct_A	324	Protein (modification methylase HAEIII); enzyme, c	3e-09
3me5_A	482	Cytosine-specific methyltransferase; structural ge	2e-06
3g7u_A	376	Cytosine-specific methyltransferase; DNA-binding,	5e-06
1vt4_I	1221	APAF-1 related killer DARK; drosophila apoptosome,	7e-06
2qrv_A	295	DNA (cytosine-5)-methyltransferase 3A; DNA methylt	1e-04
2c7p_A	327	Modification methylase HHAI; DNA methyltransferase	1e-04
3swr_A	1002	DNA (cytosine-5)-methyltransferase 1; epigenetics,	5e-04
3av4_A	1330	DNA (cytosine-5)-methyltransferase 1; CXXC-type zi	9e-04

>1g55_A DNA cytosine methyltransferase DNMT2; human DNA methyltransferase homologue; HET: DNA SAH; 1.80A {Homo sapiens} SCOP: c.66.1.26 Length = 343

 Score = 96.0 bits (239), Expect = 1e-23
 Identities = 50/130 (38%), Positives = 77/130 (59%), Gaps = 13/130 (10%)

Query: 114 VSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQ-------- 165
             ++ +L+P   + R+   +DIV P  +R  CFTK Y  Y++GTGS+L T +        
Sbjct: 213 TDVNQYLLPPKSLLRYALLLDIVQPTCRRSVCFTKGYGSYIEGTGSVLQTAEDVQVENIY 272

Query: 166 --PKNKGKASSLKEQH---LRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIA 220
               N  +   + +     LRYFTP+E+ANL  FP +F FP  ++++QRY LLGNSL++ 
Sbjct: 273 KSLTNLSQEEQITKLLILKLRYFTPKEIANLLGFPPEFGFPEKITVKQRYRLLGNSLNVH 332

Query: 221 VVAPLLQYLF 230
           VVA L++ L+
Sbjct: 333 VVAKLIKILY 342

>3qv2_A 5-cytosine DNA methyltransferase; DNMT2, ehmeth; HET: SAH; 2.15A {Entamoeba histolytica} Length = 327

 Score = 91.8 bits (228), Expect = 3e-22
 Identities = 50/117 (42%), Positives = 71/117 (60%), Gaps = 3/117 (2%)

Query: 117 DHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKG--KASS 174
           + + +P  LI + G   DIV  D KR CCFTKSY + V+GTGS+   ++P      KA  
Sbjct: 206 ESYSIPSDLILKKGMLFDIVGKDDKRTCCFTKSYTKIVEGTGSIYCPIEPHFIPVKKAED 265

Query: 175 LKEQHLRYFTPREVANLHSFPGDFQFP-HHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
           L  ++LRYFTP E+  +H F  +F      L+ +Q+Y  LGNS+S  V+A L++YLF
Sbjct: 266 LLNKNLRYFTPNEIKKIHGFSSNFTTQIDGLTDKQQYQCLGNSVSCFVIAQLMEYLF 322

>1dct_A Protein (modification methylase HAEIII); enzyme, cytosine methylase, transferase/DNA complex; HET: DNA C49 5CM; 2.80A {Haemophilus influenzae biotypeaegyptius} SCOP: c.66.1.26 Length = 324

 Score = 55.1 bits (133), Expect = 3e-09
 Identities = 26/100 (26%), Positives = 38/100 (38%), Gaps = 6/100 (6%)

Query: 134 DIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHS 193
           + V   ++       S  +      + +     KN  K    KE   R  T RE A +  
Sbjct: 226 NRVRQWNEPAFTVQASGRQCQLHPQAPVMLKVSKNLNKFVEGKEHLYRRLTVRECARVQG 285

Query: 194 FPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQA 233
           FP DF F  + SL   Y ++GN+     V   L Y  A+ 
Sbjct: 286 FPDDFIFH-YESLNDGYKMIGNA-----VPVNLAYEIAKT 319

>3me5_A Cytosine-specific methyltransferase; structural genomics, protein structure initiative, NEW YORK structural genomix research consortium; 1.75A {Shigella flexneri 2A} PDB: 3lx6_A Length = 482

 Score = 47.0 bits (111), Expect = 2e-06
 Identities = 28/104 (26%), Positives = 44/104 (42%), Gaps = 6/104 (5%)

Query: 135 IVYPDSKRCCCFTKSYYRYVKGTGSLLAT----VQPKNKGKASSLKEQHLRYFTPREVAN 190
           +VYP++ +    T S   Y  G   L+         +        ++   R  TPRE A 
Sbjct: 352 MVYPNNPQSVTRTLSARYYKDGAEILIDRGWDMATGEKDFDDPLNQQHRPRRLTPRECAR 411

Query: 191 LHSF--PGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
           L  F  PG+ +F   +S  Q Y   GNS+ + V A + + L  +
Sbjct: 412 LMGFEAPGEAKFRIPVSDTQAYRQFGNSVVVPVFAAVAKLLEPK 455

>3g7u_A Cytosine-specific methyltransferase; DNA-binding, NAD-binding, structural GENO protein structure initiative, PSI; 1.75A {Escherichia coli O157} Length = 376

 Score = 45.5 bits (108), Expect = 5e-06
 Identities = 20/70 (28%), Positives = 34/70 (48%), Gaps = 10/70 (14%)

Query: 167 KNKGKASSLKEQH---LRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVA 223
           ++KG   +++  H    R  TPRE A L  FP  F+F  H++    +  +GNS     V+
Sbjct: 300 RDKGSFQAVRPIHPYHPRVITPREAARLQGFPDWFRF--HVTKWHSFRQIGNS-----VS 352

Query: 224 PLLQYLFAQA 233
           P++     + 
Sbjct: 353 PIVAEYILKG 362

>1vt4_I APAF-1 related killer DARK; drosophila apoptosome, apoptosis, programmed cell death; HET: DTP; 6.90A {Drosophila melanogaster} PDB: 3iz8_A* Length = 1221

 Score = 45.6 bits (107), Expect = 7e-06
 Identities = 43/234 (18%), Positives = 69/234 (29%), Gaps = 73/234 (31%)

Query: 11  HVNPTLGKAK-----RKPLSFRCQLLNNQLLRSPSPLLGNDDMTVI--TKHDQPD--DSW 61
           H + TL   +      K L  R Q L  ++L + +P      +++I  +  D     D+W
Sbjct: 293 HHSMTLTPDEVKSLLLKYLDCRPQDLPREVLTT-NPRR----LSIIAESIRDGLATWDNW 347

Query: 62  -----DKL---LESC----DPVE---RFLEFSNSGDQVNTETGFLSTGTAAVDDFGAAEE 106
                DKL   +ES     +P E    F   S      +  T  LS              
Sbjct: 348 KHVNCDKLTTIIESSLNVLEPAEYRKMFDRLSVFPPSAHIPTILLSL-------IWFDVI 400

Query: 107 TVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQP 166
             +V     + + L   SL+E+      I  P                         +  
Sbjct: 401 KSDVMV---VVNKLHKYSLVEKQPKESTISIPS------------------------IYL 433

Query: 167 KNKGKASSLKEQHLR----YFTPREVANLHSFP--GD---FQF-PHHLSLRQRY 210
           + K K  +    H      Y  P+   +    P   D   +    HHL   +  
Sbjct: 434 ELKVKLENEYALHRSIVDHYNIPKTFDSDDLIPPYLDQYFYSHIGHHLKNIEHP 487

>2qrv_A DNA (cytosine-5)-methyltransferase 3A; DNA methyltransferase 3A (DNMT3A) and ITS regulatory factor; HET: DNA SAH; 2.89A {Homo sapiens} Length = 295

 Score = 41.5 bits (97), Expect = 1e-04
 Identities = 13/52 (25%), Positives = 21/52 (40%)

Query: 178 QHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYL 229
           +        E+  +  FP  +    ++S   R  LLG S S+ V+  L   L
Sbjct: 237 EKEDILWCTEMERVFGFPVHYTDVSNMSRLARQRLLGRSWSVPVIRHLFAPL 288

>2c7p_A Modification methylase HHAI; DNA methyltransferase, methyltransferase, base flipping, restriction system, transferase; HET: 5CM A1P SAH EPE CIT; 1.7A {Haemophilus haemolyticus} SCOP: c.66.1.26 PDB: 10mh_A* 1m0e_A* 1mht_A* 1hmy_A* 1skm_A* 2c7o_A* 2c7q_A* 2hmy_B* 2hr1_A* 3eeo_A* 3mht_A* 4mht_A* 5mht_A* 6mht_A* 7mht_A* 8mht_A* 9mht_A* 2zcj_A* 2z6u_A* 2z6q_A* ... Length = 327

 Score = 41.6 bits (98), Expect = 1e-04
 Identities = 17/89 (19%), Positives = 26/89 (29%), Gaps = 9/89 (10%)

Query: 142 RCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFP 201
           R               G + A                  R   PRE A +  +P  ++  
Sbjct: 240 RIYSTRGIAITLSAYGGGIFAKTGGYLVNG-------KTRKLHPRECARVMGYPDSYKVH 292

Query: 202 HHLSLRQRYALLGNSLSIAVVAPLLQYLF 230
              S  Q Y   GNS+ I V+  +   + 
Sbjct: 293 P--STSQAYKQFGNSVVINVLQYIAYNIG 319

>3swr_A DNA (cytosine-5)-methyltransferase 1; epigenetics, DNA methyltransferase fold, maintenance methyla transferase; HET: DNA SFG MES; 2.49A {Homo sapiens} PDB: 3pta_A* 3pt6_A* 3pt9_A* 4da4_A* Length = 1002

 Score = 39.9 bits (92), Expect = 5e-04
 Identities = 16/74 (21%), Positives = 27/74 (36%), Gaps = 8/74 (10%)

Query: 160 LLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSI 219
            +   +P  K       EQH R  + RE A    FP  ++     ++  ++  +GN+   
Sbjct: 928 TVTNPEPMGKQGRVLHPEQH-RVVSVRECARSQGFPDTYRLFG--NILDKHRQVGNA--- 981

Query: 220 AVVAPLLQYLFAQA 233
             V P L       
Sbjct: 982 --VPPPLAKAIGLE 993

>3av4_A DNA (cytosine-5)-methyltransferase 1; CXXC-type zinc finger/C5-methyltransferase family; HET: DNA; 2.75A {Mus musculus} PDB: 3av5_A* 3av6_A* Length = 1330

 Score = 39.2 bits (90), Expect = 9e-04
 Identities = 21/102 (20%), Positives = 33/102 (32%), Gaps = 10/102 (9%)

Query: 134  DIVYPDSKRCCCFTKSYYRYVKGTGS--LLATVQPKNKGKASSLKEQHLRYFTPREVANL 191
              +     R   +   Y R          +   +P  K       EQH R  + RE A  
Sbjct: 1210 WCLPHTGNRHNHWAGLYGRLEWDGFFSTTVTNPEPMGKQGRVLHPEQH-RVVSVRECARS 1268

Query: 192  HSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQA 233
              FP  ++F    ++  R+  +GN+     V P L       
Sbjct: 1269 QGFPDSYRFFG--NILDRHRQVGNA-----VPPPLAKAIGLE 1303

Structure Templates Detected by HHsearch

Original result of HHsearch against PDB70 database

ID	Length	Definition	Probability
Query	234
3qv2_A	327	5-cytosine DNA methyltransferase; DNMT2, ehmeth; H	99.97
4h0n_A	333	DNMT2; SAH binding, transferase; HET: SAH; 2.71A {	99.96
1g55_A	343	DNA cytosine methyltransferase DNMT2; human DNA me	99.93
3me5_A	482	Cytosine-specific methyltransferase; structural ge	99.8
3g7u_A	376	Cytosine-specific methyltransferase; DNA-binding,	99.74
3ubt_Y	331	Modification methylase HAEIII; protein-DNA complex	99.72
4ft4_B	784	DNA (cytosine-5)-methyltransferase 1; chromodomain	99.72
4dkj_A	403	Cytosine-specific methyltransferase; CG-specificit	99.63
2c7p_A	327	Modification methylase HHAI; DNA methyltransferase	99.63
3swr_A	1002	DNA (cytosine-5)-methyltransferase 1; epigenetics,	99.6
3av4_A	1330	DNA (cytosine-5)-methyltransferase 1; CXXC-type zi	99.59
2qrv_A	295	DNA (cytosine-5)-methyltransferase 3A; DNA methylt	99.55

>3qv2_A 5-cytosine DNA methyltransferase; DNMT2, ehmeth; HET: SAH; 2.15A {Entamoeba histolytica}

Probab=99.97  E-value=2.4e-30  Score=236.03  Aligned_cols=156  Identities=37%  Similarity=0.604  Sum_probs=129.4

Q ss_pred             CCCCCCceeEeeecCCCCccccccccccccCCCCCCCCCcccccccCCCCCcchhhcccCcchHHHhhcccCCCCCcccc
Q 026734            8 FPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTE   87 (234)
Q Consensus         8 ~~~~~~ry~~lAk~~~~~f~~~~~n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~i~~~le~~~~~~~~~~~   87 (234)
                      .|-+|.|+|++|.|+.++||.+..                                 .....+|.+|||.          
T Consensus       166 vPQ~R~R~fivg~r~~f~fP~~~~---------------------------------~~~~~~l~d~Le~----------  202 (327)
T 3qv2_A          166 IPNSRTRYYVMARLTPFKNEIQLH---------------------------------QEKESMISNYLDN----------  202 (327)
T ss_dssp             CSBCCCEEEEEEESSCCCSCCCCC---------------------------------CCSCCCGGGGCCS----------
T ss_pred             CCccceEEEEEEEeCCCCCCCccc---------------------------------ccccccHHHHhcc----------
Confidence            799999999999998754432210                                 0123568899983          


Q ss_pred             ccccccCCccccccCCccccchhcccccccccCCChhHHhhhcCcceeeccCCCCCcccccceeeeecCCcceeeecCC-
Q 026734           88 TGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQP-  166 (234)
Q Consensus        88 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~y~v~~~~l~k~g~~~DIV~p~s~~s~cFTksY~r~~~GtgSvl~~~~~-  166 (234)
                                                +..++|.++++++.+|+..+||++|++++++|||++|++|++|+|+++..... 
T Consensus       203 --------------------------~~~~~y~l~~~~~~~~~~~~di~~~~~~~~~~~t~~y~~y~~~~gs~l~~~~~~  256 (327)
T 3qv2_A          203 --------------------------NVNESYSIPSDLILKKGMLFDIVGKDDKRTCCFTKSYTKIVEGTGSIYCPIEPH  256 (327)
T ss_dssp             --------------------------SCCGGGBCCHHHHHHHGGGSCEEETTSSCBCCCCTTTTTSSTTSCCEEESSCSS
T ss_pred             --------------------------cccccccCCHHHHHhhhcccccccccccccccccccceEEecCCCceeeccccc
Confidence                                      11366999999999999999999999999999999999999999999854321 


Q ss_pred             -CCCCCCccccCCCcccCcHHHHHHhCCCCCCeecC-CCCCHHHHHHHhcCcccHHHHHHHHHHHHhh
Q 026734          167 -KNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFP-HHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ  232 (234)
Q Consensus       167 -~~~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp-~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~~  232 (234)
                       ........+++.++|+|||||++||||||++|.|+ ..++.+|+|+|+||||+|+|+++|++.|++.
T Consensus       257 ~~~~~~~~~~~~~~~R~lt~~E~~rlqgfP~~~~~~~~~~s~~~~y~~~GNsv~v~v~~~i~~~l~~~  324 (327)
T 3qv2_A          257 FIPVKKAEDLLNKNLRYFTPNEIKKIHGFSSNFTTQIDGLTDKQQYQCLGNSVSCFVIAQLMEYLFDD  324 (327)
T ss_dssp             CCCCSSGGGGTTSCCBCCCHHHHHHHTTCCTTCCSCCTTCCHHHHHHHHHTSCCHHHHHHHHHHHTTT
T ss_pred             ccccCCceeecCCccccCcHHHHHHhCcCCHHHcCCcCCCCHHHHHHHccCccCHHHHHHHHHHHHHH
Confidence             11223556788999999999999999999999998 6689999999999999999999999999864

>4h0n_A DNMT2; SAH binding, transferase; HET: SAH; 2.71A {Spodoptera frugiperda}

Probab=99.96  E-value=2.5e-29  Score=229.66  Aligned_cols=164  Identities=34%  Similarity=0.618  Sum_probs=130.0

Q ss_pred             CCCCCCceeEeeecCCCCccccccccccccCCCCCCCCCcccccccCCCCCcchhhcccCcchHHHhhcccCCCCCcccc
Q 026734            8 FPPHVNPTLGKAKRKPLSFRCQLLNNQLLRSPSPLLGNDDMTVITKHDQPDDSWDKLLESCDPVERFLEFSNSGDQVNTE   87 (234)
Q Consensus         8 ~~~~~~ry~~lAk~~~~~f~~~~~n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~i~~~le~~~~~~~~~~~   87 (234)
                      .|-+|.|+|++|.|+.+.|.++..+..+...|.|                       ...+.+|.+|||.+         
T Consensus       155 vPQ~R~R~fiva~r~~~~~~f~~~~~~~~~~P~~-----------------------~~~~~~l~d~Le~~---------  202 (333)
T 4h0n_A          155 VPNSRLRYYCTARRNNLTWPFKRRDEIITRLPKD-----------------------FGVPHSLESIIEED---------  202 (333)
T ss_dssp             CSCCCCEEEEEEEETTSCCCSCCCSSCBSSCSSC-----------------------CCSCCCSSTTCCSS---------
T ss_pred             CCccceEEEEEEEeCCCCCCCCcccchhhhCCCC-----------------------ccccccHHHHhccC---------
Confidence            8999999999999998777665433222121221                       12256788999841         


Q ss_pred             ccccccCCccccccCCccccchhcccccccccCCChhHHhhhcCcceeeccCCCCCcccccceeeeecCCcceeeecCC-
Q 026734           88 TGFLSTGTAAVDDFGAAEETVEVDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQP-  166 (234)
Q Consensus        88 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~y~v~~~~l~k~g~~~DIV~p~s~~s~cFTksY~r~~~GtgSvl~~~~~-  166 (234)
                                                 ..++|+++++++ ||+..+|++.|.+.+++|+|++|+++.+|+|+++.+... 
T Consensus       203 ---------------------------~~~~y~~~~~~~-~~~~~~d~~~~~~~~~~~~~k~~~~~~~g~gs~~~~~~~~  254 (333)
T 4h0n_A          203 ---------------------------VDEKFLVPEKML-RCAKVFDICYKTSKRSCCFTKAYTHYADGTGSIFTDKPRE  254 (333)
T ss_dssp             ---------------------------CCGGGBCCHHHH-TTGGGCCEECTTCSCCCCCCTTBTTBSSSSCCEECSSCHH
T ss_pred             ---------------------------CcccccCCHHHH-HHHHHhccCChhhhhhhhhccccceEEeccCceecccccc
Confidence                                       136799999998 899999999999999999999999999999998754310 


Q ss_pred             ---------C-----CCCCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          167 ---------K-----NKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       167 ---------~-----~~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                               .     .......+++.++|+|||||++||||||++|+|+..++.+|+|+|+||||+|+||++|++.|++
T Consensus       255 ~~~~~~~~~~~~~~G~~~~~~~~~~~~~R~lt~~E~~rl~gfp~~~~~~~~~s~~~~y~~~GNsv~v~v~~~i~~~l~~  333 (333)
T 4h0n_A          255 VVQKCYAAAAQNEIGGEKFVELFKELKLRYFTPKEVLMIMCFPKSYNLPTNISMKQCYRLLGNSVNVKVISELLKILFE  333 (333)
T ss_dssp             HHHHHHHHGGGSCTTCHHHHHHHHTTTCBCCCHHHHHHHTTCCTTCCCCTTSCHHHHHHHHHTSCCHHHHHHHHHHHHC
T ss_pred             chhhhhcccccCCCCcccceeeccCCCcCCCCHHHHHHhCCCCccccCCCCCCHHHHHHHhCCccCHHHHHHHHHHHhC
Confidence                     0     0001223467899999999999999999999998778999999999999999999999999974

>1g55_A DNA cytosine methyltransferase DNMT2; human DNA methyltransferase homologue; HET: DNA SAH; 1.80A {Homo sapiens} SCOP: c.66.1.26

Probab=99.93  E-value=1.7e-26  Score=210.42  Aligned_cols=115  Identities=43%  Similarity=0.826  Sum_probs=97.5

Q ss_pred             cccCCChhHHhhhcCcceeeccCCCCCcccccceeeeecCCcceeeecCC-CC------------CCCCccccCCCcccC
Q 026734          117 DHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQP-KN------------KGKASSLKEQHLRYF  183 (234)
Q Consensus       117 ~~y~v~~~~l~k~g~~~DIV~p~s~~s~cFTksY~r~~~GtgSvl~~~~~-~~------------~~k~~~l~~~~lR~l  183 (234)
                      .+|.++++.+.+|...+|++.|.+...+|++++|+++..|+|+++.+... .+            ......+|+.+.|+|
T Consensus       216 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~g~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~h~~~~R~l  295 (343)
T 1g55_A          216 NQYLLPPKSLLRYALLLDIVQPTCRRSVCFTKGYGSYIEGTGSVLQTAEDVQVENIYKSLTNLSQEEQITKLLILKLRYF  295 (343)
T ss_dssp             HHHBCCHHHHHHHGGGCCEECTTCSCCCCCCTTTTTCCTTSCCEEECCSSCCHHHHHHTTTTCCHHHHHHHHHTTCEECC
T ss_pred             hhccCCHHHHHHHHHHhCCCCccccccchhhccceeeecccccccccccccchhhccccccCCCCccceeccCCCCcccc
Confidence            45889999999999999999999889999999999999999987753210 00            001234678899999


Q ss_pred             cHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          184 TPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       184 TprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      ||||+||||||||+|.|+..++.+++|+||||||||+|+++|++.|++
T Consensus       296 T~RE~aRLqgFPd~f~f~g~~s~~~~ykqiGNAVpv~v~~~I~~~l~~  343 (343)
T 1g55_A          296 TPKEIANLLGFPPEFGFPEKITVKQRYRLLGNSLNVHVVAKLIKILYE  343 (343)
T ss_dssp             CHHHHHHHTTCCTTCCCCTTSCHHHHHHHHHHSCCHHHHHHHHHHHHC
T ss_pred             CHHHHHHHcCCChhhccCCCCCHHHHHHHhcCcccHHHHHHHHHHHhC
Confidence            999999999999999998657999999999999999999999999974

>3me5_A Cytosine-specific methyltransferase; structural genomics, protein structure initiative, NEW YORK structural genomix research consortium; 1.75A {Shigella flexneri 2A} PDB: 3lx6_A

Probab=99.80  E-value=9.9e-20  Score=174.26  Aligned_cols=112  Identities=29%  Similarity=0.417  Sum_probs=76.2

Q ss_pred             cccCCChhHH---hh---------hcCcceeeccCCCCCcccccceeeeecCCcceeeecCCC-CCCCCcc----ccCCC
Q 026734          117 DHFLVPLSLI---ER---------WGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPK-NKGKASS----LKEQH  179 (234)
Q Consensus       117 ~~y~v~~~~l---~k---------~g~~~DIV~p~s~~s~cFTksY~r~~~GtgSvl~~~~~~-~~~k~~~----l~~~~  179 (234)
                      ++|+++++..   .+         .|+.++++.|++.+++|+|++ ++|..++..+++..... ..+...+    ....+
T Consensus       322 ~ky~l~~~~~~~l~~~~~~~~~~g~gf~~~i~~~~~~~~~~~Ti~-a~~~k~gs~~~i~~~~~~~~~~~~~~~~~~~~~~  400 (482)
T 3me5_A          322 AKYILTPVLWKYLYRYAKKHQARGNGFGYGMVYPNNPQSVTRTLS-ARYYKDGAEILIDRGWDMATGEKDFDDPLNQQHR  400 (482)
T ss_dssp             GGGBCCHHHHHHHHHHHHC----------CEECTTSGGGGTCCBC-CC---CCSSSEECCCCCHHHHHHCTTCTTGGGGC
T ss_pred             cccccCHHHHHHHHHHHHhhhcccCCcccceecCCcccccceeeE-EeeeccCcceeecccccccCCccccccccccCCC
Confidence            4577777652   22         255689999998889999999 56655555555432100 0000011    12357


Q ss_pred             cccCcHHHHHHhCCCC----CCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          180 LRYFTPREVANLHSFP----GDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       180 lR~lTprE~aRLQgFP----ddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      +|+|||+||+||||||    ++|.++  ++.+++|+|+||||||+|+++|++.|.+
T Consensus       401 ~R~lTprE~~rlqgFp~~~~~~~~~~--~s~~~~y~q~GNsV~v~v~~~i~~~l~~  454 (482)
T 3me5_A          401 PRRLTPRECARLMGFEAPGEAKFRIP--VSDTQAYRQFGNSVVVPVFAAVAKLLEP  454 (482)
T ss_dssp             CEECCHHHHHHHHTSSCTTCCCSCCC--SCHHHHHHHHHTSCCHHHHHHHHHHHHH
T ss_pred             cccCCHHHHHHHcCCCCccccceecc--CCHHHHHHHcCCccChHHHHHHHHHHHH
Confidence            8999999999999999    466665  6999999999999999999999999964

>3g7u_A Cytosine-specific methyltransferase; DNA-binding, NAD-binding, structural GENO protein structure initiative, PSI; 1.75A {Escherichia coli O157}

Probab=99.74  E-value=1.9e-18  Score=159.88  Aligned_cols=85  Identities=24%  Similarity=0.280  Sum_probs=69.5

Q ss_pred             ccccceeeeecCCcceeeecCCCCC----CCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHH
Q 026734          145 CFTKSYYRYVKGTGSLLATVQPKNK----GKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIA  220 (234)
Q Consensus       145 cFTksY~r~~~GtgSvl~~~~~~~~----~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~  220 (234)
                      .++..|+|+.++..|.+++..+...    ...+++|+.+.|.||+||+||||||||+|.|.  ++.+++|+||||||||+
T Consensus       277 ~~~~~~~Rl~~d~ps~Ti~~~~~~~~g~~~~~r~~HP~~~R~lTvRE~ARlQsFPD~f~f~--g~~~~~~~qIGNAVPp~  354 (376)
T 3g7u_A          277 DKISRSTRLDPNGFCPTLRAGTARDKGSFQAVRPIHPYHPRVITPREAARLQGFPDWFRFH--VTKWHSFRQIGNSVSPI  354 (376)
T ss_dssp             CTTTCCBBCCTTSCBCCC-------------CCCBCSSSSSBCCHHHHHHHHTCCTTCCCC--SSHHHHHHHHHTSCCHH
T ss_pred             ccCCceeecCcccCCCceecCcCCCCCcccCccccCCccCcCCCHHHHHHhCCCCcceEEC--CChHHhheeeecCCCHH
Confidence            4566799999999998776543211    12466999999999999999999999999997  59999999999999999


Q ss_pred             HHHHHHHHHHh
Q 026734          221 VVAPLLQYLFA  231 (234)
Q Consensus       221 vv~~La~~Ll~  231 (234)
                      ++++|++.|.+
T Consensus       355 la~aia~~I~~  365 (376)
T 3g7u_A          355 VAEYILKGLYN  365 (376)
T ss_dssp             HHHHHHHHHHH
T ss_pred             HHHHHHHHHHH
Confidence            99999999975

>3ubt_Y Modification methylase HAEIII; protein-DNA complex, DNA cytosine-5 methyltransferase, DNA B S-adenosyl methionine binding; HET: ATP 2PE; 2.50A {Haemophilus aegyptius} PDB: 1dct_A*

Probab=99.72  E-value=7e-18  Score=149.54  Aligned_cols=56  Identities=30%  Similarity=0.501  Sum_probs=50.9

Q ss_pred             ccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          175 LKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       175 l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      .++...|+|||||+||||||||+|.|+. .+.+++|+||||||||+|+++|++.|.+
T Consensus       267 ~~~~~~R~LT~rE~aRLQgFPd~f~f~~-~s~~~~ykqiGNAVpp~la~~I~~~I~~  322 (331)
T 3ubt_Y          267 GKEHLYRRLTVRECARVQGFPDDFIFHY-ESLNDGYKMIGNAVPVNLAYEIAKTIKS  322 (331)
T ss_dssp             TCGGGCCBCBHHHHHHHHTCCTTCCCCC-SBHHHHHHHHHTSCCHHHHHHHHHHHHH
T ss_pred             CCCCcCcCCCHHHHHHhCCCCCCCEeCC-CCHHHHhhhCccCccHHHHHHHHHHHHH
Confidence            3556789999999999999999999974 4899999999999999999999999975

>4ft4_B DNA (cytosine-5)-methyltransferase 1; chromodomain, BAH domain, DNA methyltransferase domain, H3K9 binding, methylation, transferase; HET: DNA MLY SAH; 2.70A {Zea mays} PDB: 4ft2_A* 4fsx_A*

Probab=99.72  E-value=5.3e-18  Score=167.59  Aligned_cols=84  Identities=19%  Similarity=0.235  Sum_probs=71.5

Q ss_pred             ccccceeeeecCCcceeeecCCCCCCCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHH
Q 026734          145 CFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAP  224 (234)
Q Consensus       145 cFTksY~r~~~GtgSvl~~~~~~~~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~  224 (234)
                      .+...|+++.++..+.+++... ......++|+.+.|+||||||||||||||+|.|.  ++.+++|+||||||||+|+++
T Consensus       656 ~~~~~y~rl~~d~ps~TItt~~-~~~~~~~iHp~~~R~LTpRE~ARLQgFPD~y~f~--Gs~~~~ykQIGNAVpp~lA~a  732 (784)
T 4ft4_B          656 KSLKPFGRLWWDETVPTVVTRA-EPHNQVIIHPTQARVLTIRENARLQGFPDYYRLF--GPIKEKYIQVGNAVAVPVARA  732 (784)
T ss_dssp             TCCCTTEECCTTCCCSCCCSCC-CSSSSEEECSSSSSBCCHHHHHHHTTCCTTCCCC--SCHHHHHHHHHHSCCHHHHHH
T ss_pred             CcCccceeccCCCcccceeccc-cCCCCeecCCCCCcCCcHHHHHHHCCCCCCCEeC--CCHHHHHhhccCCCCHHHHHH
Confidence            3556688888887777665542 3445678999999999999999999999999997  599999999999999999999


Q ss_pred             HHHHHHh
Q 026734          225 LLQYLFA  231 (234)
Q Consensus       225 La~~Ll~  231 (234)
                      |++.|.+
T Consensus       733 Ig~~I~~  739 (784)
T 4ft4_B          733 LGYCLGQ  739 (784)
T ss_dssp             HHHHHHH
T ss_pred             HHHHHHH
Confidence            9999964

>4dkj_A Cytosine-specific methyltransferase; CG-specificity, DNA intercalation, CPG sequence, cytosine C5 methylation; HET: DNA C37 5CM SAH; 2.15A {Mycoplasma penetrans}

Probab=99.63  E-value=6.4e-17  Score=151.65  Aligned_cols=56  Identities=25%  Similarity=0.387  Sum_probs=49.9

Q ss_pred             cCCCcccCcHHHHHHhCCCCC-CeecCC---CCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          176 KEQHLRYFTPREVANLHSFPG-DFQFPH---HLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       176 ~~~~lR~lTprE~aRLQgFPd-df~fp~---~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      .+.++|+||||||+||||||| +|.++.   .++.+++|+|+||||+|+|+++|++.|++
T Consensus       331 ~~~~~R~ltprE~~rlqGFpd~~~~~~~~~~~~s~~~~y~~~GNsv~v~v~~~i~~~i~~  390 (403)
T 4dkj_A          331 TQQGVRYLTPLECFKYMQFDVNDFKKVQSTNLISENKMIYIAGNSIPVKILEAIFNTLEF  390 (403)
T ss_dssp             ETTEEEECCHHHHHHHTTCCHHHHHHHHHTSCSCHHHHHHHHHTSCCHHHHHHHHHTCCC
T ss_pred             cCCCcccCCHHHHHHHcCCCHHHhhhhhccCCCCHHHHHhhcCCccCHHHHHHHHHHHHH
Confidence            356799999999999999999 698852   26999999999999999999999999874

>2c7p_A Modification methylase HHAI; DNA methyltransferase, methyltransferase, base flipping, restriction system, transferase; HET: 5CM A1P SAH EPE CIT; 1.7A {Haemophilus haemolyticus} SCOP: c.66.1.26 PDB: 10mh_A* 1m0e_A* 1mht_A* 1hmy_A* 1skm_A* 2c7o_A* 2c7q_A* 2hmy_B* 2hr1_A* 3eeo_A* 3mht_A* 4mht_A* 5mht_A* 6mht_A* 7mht_A* 8mht_A* 9mht_A* 2zcj_A* 2z6u_A* 2z6q_A* ...

Probab=99.63  E-value=2e-16  Score=143.77  Aligned_cols=56  Identities=27%  Similarity=0.399  Sum_probs=51.9

Q ss_pred             cccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          174 SLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       174 ~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      .+++.+.|+||+||+||||||||+|.|+  ++.+++|+||||||||+|+++|++.|++
T Consensus       265 ~~~~~~~R~LT~rE~aRLQgFPd~f~f~--gs~~~~ykqIGNAVp~~l~~~Ia~~i~~  320 (327)
T 2c7p_A          265 YLVNGKTRKLHPRECARVMGYPDSYKVH--PSTSQAYKQFGNSVVINVLQYIAYNIGS  320 (327)
T ss_dssp             EEETTEEEECCHHHHHHHTTCCTTSCCC--SSHHHHHHHHHHCCCHHHHHHHHHHHHH
T ss_pred             cCCCCCCcCCCHHHHHHHCCCCcCcEeC--CCHHHHHhHccCCCCHHHHHHHHHHHHH
Confidence            3457789999999999999999999998  5999999999999999999999999975

>3swr_A DNA (cytosine-5)-methyltransferase 1; epigenetics, DNA methyltransferase fold, maintenance methyla transferase; HET: DNA SFG MES; 2.49A {Homo sapiens} PDB: 3pta_A* 3pt6_A* 3pt9_A* 4da4_A*

Probab=99.60  E-value=8.8e-16  Score=157.83  Aligned_cols=86  Identities=15%  Similarity=0.127  Sum_probs=73.2

Q ss_pred             cccccceeeeecCCc-ceeeecCCCCCCCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHH
Q 026734          144 CCFTKSYYRYVKGTG-SLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVV  222 (234)
Q Consensus       144 ~cFTksY~r~~~Gtg-Svl~~~~~~~~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv  222 (234)
                      ..++..|+|+.|+.. +.+++.........+++|+.+.|.||+||+||||||||+|.|.  ++..++|+||||||||+|+
T Consensus       910 ~~~~~~ygRL~wd~p~~~~it~~~~~~~~g~~~Hp~~~R~lt~rE~arlQ~fPd~~~f~--g~~~~~~~qiGNaVp~~~~  987 (1002)
T 3swr_A          910 NHWAGLYGRLEWDGFFSTTVTNPEPMGKQGRVLHPEQHRVVSVRECARSQGFPDTYRLF--GNILDKHRQVGNAVPPPLA  987 (1002)
T ss_dssp             GGGTTTTCBCCTTSCBSSCCSSCCTTCTTCCCBCSSSSSBCCHHHHHHHTTCCTTCCCC--SSHHHHHHHHHHSCCHHHH
T ss_pred             cccCccceeeccCCCCCeeEecCCCCCCCCcccCcccccCCCHHHHHHhCCCCcceEEc--CChHHHheeeeccCCHHHH
Confidence            346778999999988 5555543222345788999999999999999999999999997  4999999999999999999


Q ss_pred             HHHHHHHHh
Q 026734          223 APLLQYLFA  231 (234)
Q Consensus       223 ~~La~~Ll~  231 (234)
                      ++|++.|.+
T Consensus       988 ~~i~~~i~~  996 (1002)
T 3swr_A          988 KAIGLEIKL  996 (1002)
T ss_dssp             HHHHHHHHH
T ss_pred             HHHHHHHHH
Confidence            999999875

>3av4_A DNA (cytosine-5)-methyltransferase 1; CXXC-type zinc finger/C5-methyltransferase family; HET: DNA; 2.75A {Mus musculus} PDB: 3av5_A* 3av6_A*

Probab=99.59  E-value=1.1e-15  Score=160.56  Aligned_cols=84  Identities=18%  Similarity=0.171  Sum_probs=67.2

Q ss_pred             cccceeeeecCCc-ceeeecCCCCCCCCccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHH
Q 026734          146 FTKSYYRYVKGTG-SLLATVQPKNKGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAP  224 (234)
Q Consensus       146 FTksY~r~~~Gtg-Svl~~~~~~~~~k~~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~  224 (234)
                      ++..|+|+.++.. +.+++..........++|+.+.|.||+||+||||||||+|.|.  ++.+++|+||||||||+|+++
T Consensus      1222 ~~~~ygRL~~d~p~sttiT~~~~~~k~g~~iHP~q~R~LTVREaARLQsFPDdF~F~--Gs~t~~yrQIGNAVPPlLAkA 1299 (1330)
T 3av4_A         1222 WAGLYGRLEWDGFFSTTVTNPEPMGKQGRVLHPEQHRVVSVRECARSQGFPDSYRFF--GNILDRHRQVGNAVPPPLAKA 1299 (1330)
T ss_dssp             GTTTTCBCCTTSCBSSCCSSCCTTSSSCCCBCSSSSSBCCHHHHHHHTTCCTTCCCC--SSHHHHHHHHHHSCCHHHHHH
T ss_pred             ccccceecccCCCCCceEeCCcccCCCCcccCccccccCCHHHHHHhcCCCCCeEEC--CCHhhhhEEeEeCcCHHHHHH
Confidence            4556888888764 3333321001123467999999999999999999999999998  489999999999999999999


Q ss_pred             HHHHHHh
Q 026734          225 LLQYLFA  231 (234)
Q Consensus       225 La~~Ll~  231 (234)
                      ||+.|.+
T Consensus      1300 IA~~I~~ 1306 (1330)
T 3av4_A         1300 IGLEIKL 1306 (1330)
T ss_dssp             HHHHHHH
T ss_pred             HHHHHHH
Confidence            9999975

>2qrv_A DNA (cytosine-5)-methyltransferase 3A; DNA methyltransferase 3A (DNMT3A) and ITS regulatory factor; HET: DNA SAH; 2.89A {Homo sapiens}

Probab=99.55  E-value=1.3e-15  Score=137.17  Aligned_cols=57  Identities=23%  Similarity=0.278  Sum_probs=52.7

Q ss_pred             ccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          175 LKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       175 l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      +++.+.|+||+||+||||||||+|.|+...+.+++|+||||||||+|+++|++.|.+
T Consensus       234 ~~~~~~R~lt~rE~arlqgFPd~~~~~~~~s~~~~~~qiGNaVpv~~~~~i~~~i~~  290 (295)
T 2qrv_A          234 FMNEKEDILWCTEMERVFGFPVHYTDVSNMSRLARQRLLGRSWSVPVIRHLFAPLKE  290 (295)
T ss_dssp             EETTEEECCCHHHHHHHHTCCTTTTCCTTCCHHHHHHHHHTSCCHHHHHHHHGGGGG
T ss_pred             ccCCCcCCCCHHHHHHHcCCCHHHeeCCCcCHHHHhccEecCcCHHHHHHHHHHHHH
Confidence            456778999999999999999999998767889999999999999999999999875

---

Homologous Structure Domains

Structure Domains Detected by RPS-BLAST

Original result of RPS-BLAST against SCOP70(version1.75) database

ID	Length	Definition	E-value
Query	234
d1g55a_	343	c.66.1.26 (A:) DNMT2 {Human (Homo sapiens) [TaxId:	5e-20
d1dcta_	324	c.66.1.26 (A:) DNA methylase HaeIII {Haemophilus a	5e-10
d2c7pa1	327	c.66.1.26 (A:1-327) DNA methylase HhaI {Haemophilu	3e-09

>d1g55a_ c.66.1.26 (A:) DNMT2 {Human (Homo sapiens) [TaxId: 9606]} Length = 343

class: Alpha and beta proteins (a/b)
fold: S-adenosyl-L-methionine-dependent methyltransferases
superfamily: S-adenosyl-L-methionine-dependent methyltransferases
family: C5 cytosine-specific DNA methylase, DCM
domain: DNMT2
species: Human (Homo sapiens) [TaxId: 9606]

 Score = 84.7 bits (208), Expect = 5e-20
 Identities = 50/134 (37%), Positives = 79/134 (58%), Gaps = 13/134 (9%)

Query: 110 VDRCVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKN- 168
           ++    ++ +L+P   + R+   +DIV P  +R  CFTK Y  Y++GTGS+L T +    
Sbjct: 209 LEDDTDVNQYLLPPKSLLRYALLLDIVQPTCRRSVCFTKGYGSYIEGTGSVLQTAEDVQV 268

Query: 169 ------------KGKASSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNS 216
                       + + + L    LRYFTP+E+ANL  FP +F FP  ++++QRY LLGNS
Sbjct: 269 ENIYKSLTNLSQEEQITKLLILKLRYFTPKEIANLLGFPPEFGFPEKITVKQRYRLLGNS 328

Query: 217 LSIAVVAPLLQYLF 230
           L++ VVA L++ L+
Sbjct: 329 LNVHVVAKLIKILY 342

>d1dcta_ c.66.1.26 (A:) DNA methylase HaeIII {Haemophilus aegyptius [TaxId: 197575]} Length = 324

class: Alpha and beta proteins (a/b)
fold: S-adenosyl-L-methionine-dependent methyltransferases
superfamily: S-adenosyl-L-methionine-dependent methyltransferases
family: C5 cytosine-specific DNA methylase, DCM
domain: DNA methylase HaeIII
species: Haemophilus aegyptius [TaxId: 197575]

 Score = 56.2 bits (134), Expect = 5e-10
 Identities = 25/99 (25%), Positives = 36/99 (36%), Gaps = 5/99 (5%)

Query: 134 DIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKASSLKEQHLRYFTPREVANLHS 193
           + V   ++       S  +      + +     KN  K    KE   R  T RE A +  
Sbjct: 226 NRVRQWNEPAFTVQASGRQCQLHPQAPVMLKVSKNLNKFVEGKEHLYRRLTVRECARVQG 285

Query: 194 FPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFAQ 232
           FP DF F  + SL   Y ++GN    AV   L   +   
Sbjct: 286 FPDDFIFH-YESLNDGYKMIGN----AVPVNLAYEIAKT 319

>d2c7pa1 c.66.1.26 (A:1-327) DNA methylase HhaI {Haemophilus haemolyticus [TaxId: 726]} Length = 327

class: Alpha and beta proteins (a/b)
fold: S-adenosyl-L-methionine-dependent methyltransferases
superfamily: S-adenosyl-L-methionine-dependent methyltransferases
family: C5 cytosine-specific DNA methylase, DCM
domain: DNA methylase HhaI
species: Haemophilus haemolyticus [TaxId: 726]

 Score = 53.8 bits (128), Expect = 3e-09
 Identities = 21/125 (16%), Positives = 42/125 (33%), Gaps = 7/125 (5%)

Query: 113 CVSIDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPKNKGKA 172
              ++H ++    +      ++   P + R     K        +   +A       G  
Sbjct: 199 DSEVEHLVIDRKDLVMTNQEIEQTTPKTVRLGIVGKGGQGERIYSTRGIAITLSAYGGGI 258

Query: 173 SSLKEQHL-----RYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQ 227
            +    +L     R   PRE A +  +P  ++     S  Q Y   GNS+ I V+  +  
Sbjct: 259 FAKTGGYLVNGKTRKLHPRECARVMGYPDSYKVHP--STSQAYKQFGNSVVINVLQYIAY 316

Query: 228 YLFAQ 232
            + + 
Sbjct: 317 NIGSS 321

Homologous Domains Detected by HHsearch

Original result of HHsearch against SCOP70(version1.75) database

ID	Length	Definition	Probability
Query	234
d1g55a_	343	DNMT2 {Human (Homo sapiens) [TaxId: 9606]}	99.9
d1dcta_	324	DNA methylase HaeIII {Haemophilus aegyptius [TaxId	99.74
d2c7pa1	327	DNA methylase HhaI {Haemophilus haemolyticus [TaxI	99.73

>d1g55a_ c.66.1.26 (A:) DNMT2 {Human (Homo sapiens) [TaxId: 9606]}

class: Alpha and beta proteins (a/b)
fold: S-adenosyl-L-methionine-dependent methyltransferases
superfamily: S-adenosyl-L-methionine-dependent methyltransferases
family: C5 cytosine-specific DNA methylase, DCM
domain: DNMT2
species: Human (Homo sapiens) [TaxId: 9606]

Probab=99.90  E-value=1.5e-24  Score=187.64  Aligned_cols=116  Identities=43%  Similarity=0.827  Sum_probs=98.0

Q ss_pred             ccccCCChhHHhhhcCcceeeccCCCCCcccccceeeeecCCcceeeecCCC-------------CCCCCccccCCCccc
Q 026734          116 IDHFLVPLSLIERWGSAMDIVYPDSKRCCCFTKSYYRYVKGTGSLLATVQPK-------------NKGKASSLKEQHLRY  182 (234)
Q Consensus       116 ~~~y~v~~~~l~k~g~~~DIV~p~s~~s~cFTksY~r~~~GtgSvl~~~~~~-------------~~~k~~~l~~~~lR~  182 (234)
                      ..+|.++.+...++....+.+.+...++.|+++.++.+..+.++.+.+....             .......+|+.+.|+
T Consensus       215 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~R~  294 (343)
T d1g55a_         215 VNQYLLPPKSLLRYALLLDIVQPTCRRSVCFTKGYGSYIEGTGSVLQTAEDVQVENIYKSLTNLSQEEQITKLLILKLRY  294 (343)
T ss_dssp             HHHHBCCHHHHHHHGGGCCEECTTCSCCCCCCTTTTTCCTTSCCEEECCSSCCHHHHHHTTTTCCHHHHHHHHHTTCEEC
T ss_pred             hhhcccChHhhhhhhhhcccccccccccchhcccccccccCCCcccccccccccCCcccccccCCCCCCcccCCCCCCcC
Confidence            3668899999999999999999999999999999998888888776543210             000123357889999


Q ss_pred             CcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          183 FTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       183 lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      |||||+||||||||+|.|+..++.+++|+||||||||+|+++|++.|++
T Consensus       295 LT~rE~arLQgFPd~f~~~~~~s~~~~ykqiGNaVp~~v~~~i~k~l~e  343 (343)
T d1g55a_         295 FTPKEIANLLGFPPEFGFPEKITVKQRYRLLGNSLNVHVVAKLIKILYE  343 (343)
T ss_dssp             CCHHHHHHHTTCCTTCCCCTTSCHHHHHHHHHHSCCHHHHHHHHHHHHC
T ss_pred             CCHHHHHHhCCCCCCCEECCCCCHHHHhhhEeCCccHHHHHHHHHHhcC
Confidence            9999999999999999999878999999999999999999999999975

>d1dcta_ c.66.1.26 (A:) DNA methylase HaeIII {Haemophilus aegyptius [TaxId: 197575]}

class: Alpha and beta proteins (a/b)
fold: S-adenosyl-L-methionine-dependent methyltransferases
superfamily: S-adenosyl-L-methionine-dependent methyltransferases
family: C5 cytosine-specific DNA methylase, DCM
domain: DNA methylase HaeIII
species: Haemophilus aegyptius [TaxId: 197575]

Probab=99.74  E-value=7.5e-19  Score=150.01  Aligned_cols=56  Identities=30%  Similarity=0.501  Sum_probs=51.2

Q ss_pred             ccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          175 LKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       175 l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      .++...|+|||||+||||||||+|.|+. .+.+++|+||||||||+|+++|++.|.+
T Consensus       267 ~~~~~~R~LT~rE~arLqgFPd~~~f~~-~s~~~~y~~iGNaV~p~v~~~I~~~I~~  322 (324)
T d1dcta_         267 GKEHLYRRLTVRECARVQGFPDDFIFHY-ESLNDGYKMIGNAVPVNLAYEIAKTIKS  322 (324)
T ss_dssp             TTGGGCCBCBHHHHHHHHTCCTTCCCCC-SBHHHHHHHHHTSCCHHHHHHHHHHHHH
T ss_pred             CCCCcCcCCCHHHHHHHCCCCCCCEeCC-CCHHHHHeeEccCccHHHHHHHHHHHHH
Confidence            3566789999999999999999999974 4899999999999999999999999975

>d2c7pa1 c.66.1.26 (A:1-327) DNA methylase HhaI {Haemophilus haemolyticus [TaxId: 726]}

class: Alpha and beta proteins (a/b)
fold: S-adenosyl-L-methionine-dependent methyltransferases
superfamily: S-adenosyl-L-methionine-dependent methyltransferases
family: C5 cytosine-specific DNA methylase, DCM
domain: DNA methylase HhaI
species: Haemophilus haemolyticus [TaxId: 726]

Probab=99.73  E-value=2.1e-18  Score=149.36  Aligned_cols=57  Identities=26%  Similarity=0.394  Sum_probs=52.9

Q ss_pred             ccccCCCcccCcHHHHHHhCCCCCCeecCCCCCHHHHHHHhcCcccHHHHHHHHHHHHh
Q 026734          173 SSLKEQHLRYFTPREVANLHSFPGDFQFPHHLSLRQRYALLGNSLSIAVVAPLLQYLFA  231 (234)
Q Consensus       173 ~~l~~~~lR~lTprE~aRLQgFPddf~fp~~~s~~q~ykqiGNAV~v~vv~~La~~Ll~  231 (234)
                      ..+++.+.|+|||||+||||||||+|.|+  ++.+++|+||||||||+|+++|++.|++
T Consensus       264 ~~~~~~~~R~lT~rE~arLqgFPd~f~~~--~s~~~~~~~iGNAVp~~v~~~I~~~I~~  320 (327)
T d2c7pa1         264 GYLVNGKTRKLHPRECARVMGYPDSYKVH--PSTSQAYKQFGNSVVINVLQYIAYNIGS  320 (327)
T ss_dssp             EEEETTEEEECCHHHHHHHTTCCTTSCCC--SSHHHHHHHHHHCCCHHHHHHHHHHHHH
T ss_pred             cccCCCcCcCCCHHHHHHhCCCCCCcEeC--CCHHHHHhHccCCCCHHHHHHHHHHHHH
Confidence            34677899999999999999999999998  4999999999999999999999999975