Last chance! 50% off unlimited learning
Sale ends in
Plot read coverage for a GRanges object
An internal function used by plotStrandedCoverage to plot read counts within each fixed bin across the entire chromosome.
plotCoverage(x, plotLegend = FALSE, legend.cex = 1, ...)GRanges object with values slot saved for read counts within the corresponding ranges.
Binary indcator. If TRUE, legend will be plotted on the top left the plot. Legend is expected to be the chromsome name and length, which must be available in the GRange object argument.
Font size of the legend.
The read counts is plotted in blue bars as positive integer across the x-axis as the sorted positions across the chromosome. The plot can be used to examine the overall alignment properties for each chromosome.
P. Aboyoun, H. Pages and M. Lawrence (). GenomicRanges: Representation and manipulation of genomic intervals. R package version 1.8.9.
# Retrieve system files
extdata.dir <- system.file("extdata", package="RIPSeeker")
bamFiles <- list.files(extdata.dir, ".bam$", recursive=TRUE, full.names=TRUE)
bamFiles <- grep("PRC2", bamFiles, value=TRUE)
alignGal <- getAlignGal(bamFiles[1], reverseComplement=TRUE, genomeBuild="mm9")
alignGR <- as(alignGal, "GRanges")
alignGRList <- GRangesList(as.list(split(alignGR, seqnames(alignGR))))
binSize <- 1000
binGR <- binCount(alignGRList$chrX, binSize)
plotCoverage(binGR, plotLegend=TRUE)
Run the code above in your browser using DataLab