Dear all,
From bisulfite sequencing data (downloaded from NGSmethDB ) I've calculated some CpG sites showing that the number of their methylated cytosines (methylC) is more than the number of their corresponding coverages (numReads). Thus it will give the methylation value greater than 1.0. Is this a feasible solution? why?
my example is:
##########################################################
chrom chromStart chromEnd context numReads methylC meth
16144246 chr21 46872388 46872389 CG 27 28 _______ 1.0370370
16144370 chr21 46876524 46876525 CG 9 10 _______ 1.1111111
16144659 chr21 46887294 46887295 CG 13 9 _______ 0.6923077
16144722 chr21 46888779 46888780 CG 5 6 _______ 1.2000000
16144783 chr21 46890849 46890850 CG 5 6 _______ 1.2000000
16144875 chr21 46892957 46892958 CG 9 9 _______ 1.0000000
16144890 chr21 46893253 46893254 CG 18 18 _______ 1.0000000
16144898 chr21 46893401 46893402 CG 15 14 _______ 0.9333333
##########################################################
Thank you for the help!
From bisulfite sequencing data (downloaded from NGSmethDB ) I've calculated some CpG sites showing that the number of their methylated cytosines (methylC) is more than the number of their corresponding coverages (numReads). Thus it will give the methylation value greater than 1.0. Is this a feasible solution? why?
my example is:
##########################################################
chrom chromStart chromEnd context numReads methylC meth
16144246 chr21 46872388 46872389 CG 27 28 _______ 1.0370370
16144370 chr21 46876524 46876525 CG 9 10 _______ 1.1111111
16144659 chr21 46887294 46887295 CG 13 9 _______ 0.6923077
16144722 chr21 46888779 46888780 CG 5 6 _______ 1.2000000
16144783 chr21 46890849 46890850 CG 5 6 _______ 1.2000000
16144875 chr21 46892957 46892958 CG 9 9 _______ 1.0000000
16144890 chr21 46893253 46893254 CG 18 18 _______ 1.0000000
16144898 chr21 46893401 46893402 CG 15 14 _______ 0.9333333
##########################################################
Thank you for the help!
Comment