View source: R/sequencing_QCs.r
BaitCapture | R Documentation |
check bait capture efficiency. Expected to be ~70
BaitCapture(sampleSheet, genome, baits, cores = 1)
sampleSheet |
QuasR sample sheet |
genome |
BS genome |
baits |
Full path to bed file containing bait coordinates. If chromosome names are in e.g. "1" format, they'll be temporarily converted to "chr1" |
cores |
number of cores for parallel processing. Defaults to 1 |
bait capture efficiency
Qinput = paste0(tempdir(), "/NRF1Pair_Qinput.txt") library(BSgenome.Mmusculus.UCSC.mm10) if(file.exists(Qinput)){ # DO NOT RUN # BaitRegions = SingleMoleculeFootprintingData::EnrichmentRegions_mm10.rds() # BaitCaptureEfficiency = BaitCapture(sampleSheet = Qinput, genome = BSgenome.Mmusculus.UCSC.mm10, baits = BaitRegions) }
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.