Get the sliding windows that overlap a GRanges object.
getWinOverlapGRanges(x, seqInfo, winWidth = 1000L, winStep = 100L, nbOverlapBases = 1)
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
a data frame that contains some key information of the sequences
the length of the sliding window, 1000 by default.
the step length to sliding the window, 100 by default.
a window is considered to overlap with a range of
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
This method makes use of the method
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
nbOverlapBases to specify the minimum number of bases that a
window should overlap with a range of
x, instead of using proprotion
A list of two logical vectors (for positive and negative strand) defining which windows that overlap with the given GRanges object.
1 2 3 4 5 6 7
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.