This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bismark: bedgraph to bigwig for DNA methylation

Hello, I have a small question. I want to convert the bismark coverage files to bigWig format. I have the sample_CHG.gz.bismark.cov files. Other than the first 3 columns which refer to the chromosome and coordinates, which of the last 3 columns I should keep as the score for the bedgraph ? Thanks.

SL4.0ch07       136     136     50      1       1
SL4.0ch07       153     153     100     2       0
SL4.0ch07       378     378     50      1       1
SL4.0ch07       380     380     100     2       0
methylation
• 3,745 views • link
updated by aconga • aconga • written by pixie@bioinfo pixie@bioinfo

2 answers

The score is the 4th column. The 5th column reports methylated reads in that cytosine and the 6th column reports unmethylated reads in that citosine.

So you should only keep the first 4 columns when transforming that file into a bedGraph file.

• 0 views • link
written by aconga • aconga

bedGraph has four tab-separated columns, that is chr-start-end-value where value is any numeric value that you want. So you could simply use either the read counts of the Cs or a percentage or fraction of the methylated CpGs as value. From the normal bedGraph to bigwig tools apply, e.g. bedgraphtobigwig from kentUtils.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this