Question: Rsamtools - summarizeOverlaps - count read on complementary strand
5.4 years ago by
Guest User • 12k
Guest User • 12k wrote:
Hi! I have RNAseq data which are single-end (50b), strand specific, but reads are the complement of the RNA (ie reads are always on the complementary strand to their corresponding gene). I would like to count reads on genes using summarizeOverlaps, however I do not find the way to tell him to count reads that are on the complementary strand to the exons. Is there a way? (with HTseq-count, I was using the option --stranded=reverse). -- output of sessionInfo(): none -- Sent via the guest posting facility at bioconductor.org.
ADD COMMENT • link •