JunctionSeq outputs coordinates without the strand information when RNA libraries are unstranded.
We can resolve the issue by changing the IRanges::findOverlaps' parameter ignore.strand from FALSE to TRUE. As a consequence, introns from opposite strands would be grouped.
JunctionSeq outputs coordinates without the strand information when RNA libraries are unstranded.
We can resolve the issue by changing the IRanges::findOverlaps' parameter
ignore.strand
from FALSE to TRUE. As a consequence, introns from opposite strands would be grouped.