I am using a pipeline (Paleomix) which mark duplicates in BAM files using Picard MarkDuplicates, but without removing them from the outputs.
I wish to compare the coverage in these files both with and without PCR duplicates filtering.
Would you know any tool or command adapted to this task?
I have tried to use different samtools depth scripts, but they all give me different results.
Thanks for your help!