Which read should be considered to calculate methylation value?
Entering edit mode
12 months ago
irfanwustl ▴ 20

I am using BISCUIT for WGBS data. After the full process BISCUIT provides bedgraph file with CpG and corresponding methylation value. However, I am wondering how BISCUIT figure out which read should be considered in a bam file for methylation value calculation? For example, in the biscuit tview image there are 14 non duplicated reads (all has mapq=60) for the highlighted CpG. However, the generated bedgraph file is saying that the depth is 10 (5th column).

chr1 1222700 1222702 0.901 10 C:0.670:3,G:1.000:7

I am wondering what is the technique to calculate methylation value from the raw bam file for a certain CpG?

Methylation BISCUIT bam • 323 views

Login before adding your answer.

Traffic: 2520 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6