Question: Average depth/coverage before and after PCR duplicates cleanup
0
gravatar for Begonia_pavonina
3 months ago by
Begonia_pavonina20 wrote:

Hello everyone,

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!

ADD COMMENTlink modified 8 weeks ago by Biostar ♦♦ 20 • written 3 months ago by Begonia_pavonina20

Try Mosdepth.

ADD REPLYlink written 3 months ago by Arup Ghosh2.5k

Thanks for the advice, Mosdepth looks a very interesting tool! Unfortunately, on the Github page I don't see the flag or any option for allowing/ignoring the PCR duplicates annotations in the BAM file. Would you remember the command used for it?

ADD REPLYlink written 3 months ago by Begonia_pavonina20
1
gravatar for ATpoint
3 months ago by
ATpoint36k
Germany
ATpoint36k wrote:

-F --flag <FLAG> exclude reads with any of the bits in FLAG set [default: 1796]

The duplicate flag is 1024, see http://broadinstitute.github.io/picard/explain-flags.html. I second the advice for mosdepth (unfortunately only available on Linux) because it is wickedly fast.

Edit 10.5.20: I managed to build mosdepth from source on Mac, see comment at the bottom of this issue: https://github.com/brentp/mosdepth/issues/103

ADD COMMENTlink modified 8 weeks ago • written 3 months ago by ATpoint36k

That's very precious information as well ATpoint!

I did not make the relation between this flag and the annotations in the BAM file. Definitely going to use Mosdepth. Thanks!

ADD REPLYlink written 3 months ago by Begonia_pavonina20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 674 users visited in the last hour