calculateMappingBiasGatk4: Calculate Mapping Bias from GATK4 GenomicsDB

In PureCN: Copy number calling and SNV classification using targeted short read sequencing

Description

Function calculate mapping bias for each variant in the provided panel of normals GenomicsDB.

Usage

calculateMappingBiasGatk4(
  workspace,
  reference.genome,
  min.normals = 1,
  min.normals.betafit = 7,
  min.median.coverage.betafit = 5
)

Arguments

workspace	Path to the GenomicsDB created by GenomicsDBImport
reference.genome	Reference FASTA file.
min.normals	Minimum number of normals with heterozygous SNP for calculating position-specific mapping bias.
min.normals.betafit	Minimum number of normals with heterozygous SNP fitting a beta distribution
min.median.coverage.betafit	Minimum median coverage of normals with heterozygous SNP for fitting a beta distribution

Value

A GRanges object with mapping bias and number of normal samples with this variant.

Author(s)

Markus Riester

Examples

## Not run:  
resources_file <- system.file("extdata", "gatk4_pon_db.tgz", 
    package = "PureCN")
tmp_dir <- tempdir()
untar(resources_file, exdir = tmp_dir)
workspace <- file.path(tmp_dir, "gatk4_pon_db")
bias <- calculateMappingBiasGatk4(workspace, "hg19")
saveRDS(bias, "mapping_bias.rds")
unlink(tmp_dir, recursive=TRUE)

## End(Not run)

PureCN documentation built on Nov. 8, 2020, 5:37 p.m.