getWinOverlapGRanges: Get the sliding windows that overlap a GRanges object
In UofABioinformaticsHub/rnaCleanR: Calculate strandness information of a bam file

Description Usage Arguments Details Value Examples

Get the sliding windows that overlap a GRanges object.

1 2	getWinOverlapGRanges(x, seqInfo, winWidth = 1000L, winStep = 100L, nbOverlapBases = 1)

`x`	a GRanges object, which defines the coordinates of the ranges in the reference genome that all reads mapped to those ranges must be kept by the filtering method `filterDNA`.
`seqInfo`	a data frame that contains some key information of the sequences
`winWidth`	the length of the sliding window, 1000 by default.
`winStep`	the step length to sliding the window, 100 by default.
`nbOverlapBases`	a window is considered to overlap with a range of `x` if it overlaps with at least `nbOverlapBases` bases.

This finds the windows that overlaps the positive/negative strand of a GRanges object. The GRanges object, which is mustKeepRanges in the filterDNA method, defines the coordinates of the ranges in the reference genome that all reads mapped to those ranges must be kept by the filtering method filterDNA. This method makes use of the method getWinOverlapEachIRange by pretending each given range as the range of a read. Since the widths of x are not necessarily the same (as normal read lengths), we use nbOverlapBases to specify the minimum number of bases that a window should overlap with a range of x, instead of using proprotion as readProp in getWinOverlapEachIRange.

A list of two logical vectors (for positive and negative strand) defining which windows that overlap with the given GRanges object.

library(GenomicRanges)
x <- GRanges(seqnames = "10",ranges = IRanges(start = c(10000,15000),
end=c(20000,30000)),strand = c("+","-"))
seqInfo <- data.frame("Sequence"=10,"FirstBaseInPart"=1)
getWinOverlapGRanges(x,seqInfo)
seqInfo <- data.frame("Sequence"=10,"FirstBaseInPart"=10000000)
getWinOverlapGRanges(x,seqInfo)