I am running featureCounts on my ATAC-seq data. I have allowed multi-overlapping reads (i.e. reads that overlap with more than one peak). However, I am unsure whether to include the --fraction option (assign fractional counts to features) or whether I should allow each overlap to count as 1.
Does anybody have thoughts on this? What is most commonly used in this situation and what are the pros/cons of each approach?