Hi,
According to the last example of this figure reads that fall in overlapping regions of 2 features are discarded in every mode of htseq-count. I understand how this is reasonable in DE analysis. However, I am now interested in raw counts of features in my experimental setup as they are just predictions and I want to count a read that falls in an ambiguous region for both features. Is there any other tool or R package out there that allows this?
Thanks
+1 for featureCounts, it works well and is well documented.