Which read should be considered to calculate methylation value?
0
1
Entering edit mode
3.6 years ago
irfanwustl ▴ 80

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 • 721 views
ADD COMMENT

Login before adding your answer.

Traffic: 1778 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6