fraction.target {TEQC} | R Documentation |
Calculates the fraction of the reference genome that is targeted
fraction.target(targets, Offset = 0, genome = c(NA, "hg38", "hg19", "hg18"), genomesize)
targets |
|
Offset |
integer; add |
genome |
genome version targets were designed and reads aligned to. For the given options the total genome size is set automatically. For other genomes or versions, leave this option empty ('NA') and specify the genome size with option 'genomesize' |
genomesize |
integer: specify the total genome size manually. If 'genomesize' is given, option 'genome' will be ignored. |
Returns the fraction of nucleotides within the genome that were targeted.
With the output from fraction.target
and fraction.reads.target
the 'enrichment' of the target capture experiment can be calculated as
'fraction of on-target reads / fraction of target within genome'
Manuela Hummel m.hummel@dkfz.de
fraction.reads.target
, get.targets
exptPath <- system.file("extdata", package="TEQC") targetsfile <- file.path(exptPath, "ExampleSet_Targets.bed") targets <- get.targets(targetsfile, skip=0) fraction.target(targets, genome="hg19")