nf-core / methylseq

Methylation (Bisulfite-Sequencing) analysis pipeline using Bismark or bwa-meth + MethylDackel
https://nf-co.re/methylseq
MIT License
137 stars 136 forks source link

Adding --maxVariantFrac and --minOppositeDepth to MethylDackel (Bwameth option) #404

Open LucaZanella15 opened 2 months ago

LucaZanella15 commented 2 months ago

Description of feature

Hi all,

I think that adding the possibility to provide the --maxVariantFrac option (and its --minOppositeDepth counterpart, by default set to 0) would be of great value.

The possibility to switch on and off --keepDupes (currently by default properly being on) would also be useful to account for current EM-seq library preparation technologies that produce complex libraries that have low probability of fragments having the same start and stop positions unless they are true duplicates