Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Keeping low number of PCR duplicates #28

Open
phoebe460 opened this issue Apr 13, 2024 · 0 comments
Open

Keeping low number of PCR duplicates #28

phoebe460 opened this issue Apr 13, 2024 · 0 comments

Comments

@phoebe460
Copy link

Hi SICER2 Developers,

First off, thank you for creating a great peak calling tool. I am planning to use this for my own analysis for ChIP-seq data actually. In that case, I am wondering if there is a way to keep PCR duplicates using a --keep-duplicates flag, similar to what has been implemented in EPIC2 but setting it to just a low number for instance like 1, which would remove the majority of PCR duplicates found but still keep a low number of duplicates? I don't see a parameter that can do this under SICER2 arguments here.

A similar thing can be done using their keep duplicates flag in MACS3 as follows:

--keep-dup

It controls the MACS3 behavior towards duplicate tags at the exact same location – the same coordination and the same strand. You can set this as auto, all, or an integer value. The auto option makes MACS3 calculate the maximum tags at the exact same location based on binomial distribution using 1e-5 as p-value cutoff; and the all option keeps every tag. If an integer is given, at most this number of tags will be kept at the same location. The default is to keep one tag at the same location. Default: 1

If you can clarify this for me or even include a --keep-duplicates flag before I start using your program, then that would be greatly appreciated.

Thank you,
Phoebe

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant