Calling allele-specific methylation from bisulfite sequencing reads
1
0
Entering edit mode
5.3 years ago
Jautis ▴ 500

I have mapped bisulfite sequencing reads (RRBS) and a vcf file of genetic variants (called from low coverage sequencing). Can anyone suggest a pipeline for calling allele-specific methylation?

Thanks!

bisulfite-sequencing allele-specific sequencing • 1.5k views
ADD COMMENT
0
Entering edit mode
5.3 years ago

I don't know if such a tool or pipeline exists off-the-shelf. If not, I think some scripting is required, probably using python/pysam java/htsjdk or other API to process bam files.

Essentially for each read spanning a variant and spanning one or more cytosines on the reference sequence you record whether the cytosine(s) are methylated or not (i.e. they are C or they have been converted to T) and you record whether the variant is a reference or an alternative allele. You need to record also the position of the cytosine(s) in genomic coordinates. (This strategy is not difficult but it might be a bit tricky to get right.)

ADD COMMENT

Login before adding your answer.

Traffic: 2048 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