featureCounts allowMultiOverlap - fractional or not?
0
0
Entering edit mode
5.4 years ago
Lucy ▴ 140

Hi,

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?

Many thanks,

Lucy

featureCounts ATAC-seq • 2.3k views
ADD COMMENT
0
Entering edit mode

Remember that ATAC-seq relies on the integration of a transposon into open chromatin, introducing an adapter that flanks the genomic fragments. Therefore, it is actually only the ends of the paired-end fragments that are of interest to infer chromatin accessability, and therefore, every overlap should count. Still, I would be surprised if this makes a notable difference.

ADD REPLY
0
Entering edit mode

Thank you - I am now trying to perform featureCounts with the cut sites instead.

ADD REPLY

Login before adding your answer.

Traffic: 1437 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6