Skip to content

CpG Methylation Calculation from bam file #43

@Irfanwustl

Description

@Irfanwustl

Hi,
I am trying to understand which read to consider file for methylation value calculation from the bam file generated by BISCUIT? Also how to figure out which is "C strand" and "G strand" from the generated bam? For example, in the following biscuit tview image, there are 14 non duplicated reads (all has mapq=60) for the highlighted CpG.

biscuit

However, the generated bedgraph file is saying that the total depth is 10.
Screen Shot 2020-09-17 at 3 18 26 PM

Will you please explain what I am missing here? I have also attached the corresponding sample sam file (for some reason can not attach the bam file) if you want to reproduce it. Thanks

corres_pos_cgr1-1222701_nodup.txt

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions