## get reads and targets
exptPath <- system.file("extdata", package="TEQC")
readsfile <- file.path(exptPath, "ExampleSet_Reads.bed")
reads <- get.reads(readsfile, idcol=4, skip=0)
targetsfile <- file.path(exptPath, "ExampleSet_Targets.bed")
targets <- get.targets(targetsfile, skip=0)
## calculate per-base coverages
Coverage <- coverage.target(reads, targets, perBase=TRUE)
## coverage uniformity plot
coverage.uniformity(Coverage)
Run the code above in your browser using DataLab