Skip to main content

Identifying focal full methylation of cell subpopulation and inferring fraction

Project description


DNA methylation patterns within a cell population from individual somatic tissues are highly heterogeneous and polymorphic. We developed CellMethy method to quantify fraction of focal full methylation of cell subpopulation (FMC fraction) and identify fully methylated regions of cell subpopulation (CellMethy) based on single base resolution DNA methylation data. This script is used for infering fraction of focal full methylation cell subpopulation.

Easy to start:

Inputfile: File sepearated by “\t” after bismark processed including read id, strand, chromosome, position of CpGc and methylation state (Z or z).

The result would be placed into input folder: outputfile.

If your fastq mapping is done with bismark and extracted methylation state, use following command:

python -f inputfile -o outputfile


-f: The file name of input file after bismark processed, include five columns: read ID, strand, chromosome, position of CpG and methylation states Z (methylated) or z (unmethylated) separated by ‘\t’.

-c: Lowest coverage cutoff in each bin, default is 10.

-b: number of CpGs in each bin, default is 5.

-o: The file name of output file showing full methylation of cell subpopulation, include five columns: chromosome, start, end, FMC fraction, and CpGs number in the region separated by “\t”.


cd ~/CellMethy-*.*.*
python ./CellMethy/bin/ -f ./data/data_file -b5 -c 10 -o myoutput

How to get Input file

If you have fastq data, you can mapping with bismark tools.


bismark ./referenceGenome --bowtie2 test.fastq -o test.sam
bismark_methylation_extractor -s --comprehensive test.sam

The output file named “CpG_context_test.txt” is inputfile of CellMethy. The format of CellMethy inputfile include read id, strand, chromosome, position of CpGc and methylation state (Z or z) separated by ‘\t’.

Read1       +       chr21   9827508 Z
Read1       -       chr21   9827503 z
Read1       -       chr21   9827484 z
Read2       +       chr21   9827434 Z
Read2       +       chr21   9827454 Z
Read2       -       chr21   9827483 z

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

CellMethy-1.2.0.tar.gz (3.2 MB view hashes)

Uploaded source

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page