PrepareLibrary: function PrepareLibrary

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/functions.R

Description

Define the Genomic Ranges and Genomic Positions covered by the AmpliSeq™ Panel to include in the study and define SNP positions to exclude from the study. Trimming amplicon ends is performed if specified. This function is mostly useful if you want to add some SNP positions which are not existing in the positions_ranges.rda file provided within the package. It is provided to be able to reconstruct positions_ranges.rda data.

Usage

1
2
3
4
5
6
PrepareLibrary(info.dir = "Info/", bed.filename = "lungcolonV2.bed.txt",
  snp.filename = "ExAC.r0.3.sites.vep.vcf.gz",
  snp.extra = c("chr2:212812097", "chr4:1807909", "chr7:140481511",
  "chr14:105246474", "chr18:48586344", "chr19:1223055"),
  output.name = "positions_ranges.rda", output.dir = info.dir,
  load.from.broad.insitute = FALSE)

Arguments

info.dir,

char, name of the folder containing the library information files (default 'Info/')

bed.filename,

char, name of a BED table (tab-delimited) describing the Panel (with first 3 columns: "chr" (ex:chr1), "start position" (ex:115252190), "end position" (ex:115252305), i.e. the Ion AmpliSeq™ Colon and Lung Cancer Research Panel v2 (default 'lungcolonV2.bed.txt' as provided in the inst/extdata/Info folder of the package).

snp.filename,

char, name of the vcf file describing known SNP positions, obtained from ftp://ftp.broadinstitute.org/pub/ExAC_release/release0.3/ExAC.r0.3.sites.vep.vcf.gz (default 'ExAC.r0.3.sites.vep.vcf.gz'). It requires a corresponding TBI file to be in the same folder (obtained from ftp://ftp.broadinstitute.org/pub/ExAC_release/release0.3/ExAC.r0.3.sites.vep.vcf.gz.tbi)

snp.extra,

a vector of char, a vector of extra known snp positions manually curated (ex:"chrN:XXXXXXXXX")

output.name,

char, filename to save pos_ind and pos_snp (default 'positions_ranges.rda')

output.dir,

char, directory where to save pos_ind and pos_snp (default info.dir)

load.from.broad.insitute,

boolean, if TRUE load snp.filename from Broad Institute ftp server otherwise use the file positions_ranges_broad.rda (default FALSE)

Value

Save the following variables in a .rda file defined by output.name in the folder defined by output.dir:

Author(s)

N. Pécuchet, P. Laurent-Puig and Y. Rozenholc

References

Analysis of base-position error rate of next-generation sequencing to detect tumor mutations in circulating DNA N. Pécuchet, Y. Rozenholc, E. Zonta, D. Pietraz, A. Didelot, P. Combe, L. Gibault, J-B. Bachet, V. Taly, E. Fabre, H. Blons, P. Laurent-Puig in Clinical Chemistry

See Also

positions_ranges,

Examples

1
2
   bad.pos = "chr7:15478"
   PrepareLibrary(info.dir='./',snp.extra=bad.pos,output.dir=paste0(tempdir(),'/'))

PlasmaMutationDetector documentation built on June 11, 2018, 9:04 a.m.