This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to handle DNA methylation duplication in Bismark?

Hi

I have run Bismark for RRBS data, but I think it is not fix duplicates problem.

Anybody know how to correct duplicates problem for Bismark aligner?

I found deduplicate_bismark_alignment_output.pl which works as a post alignment step for DNA methylation analysis in Bismark pipeline.

Thanks

sequencing alignment

You mean PCR duplicates? You can find them with picard.

In the documentation it is stated:

deduplicate_bismark --bam [options] <filenames>

"This command will deduplicate the Bismark alignment BAM file and remove all reads but one which align to the the very same position and in the same orientation. This step is recommended for whole-genome bisulfite samples, but should not be used for reduced representation libraries such as RRBS, amplicon or target enrichment libraries."

0 answers

No answers yet.

Log in to answer this question.