type_context: Retrieve context of base substitution types

Description Usage Arguments Value See Also Examples

View source: R/type_context.R

Description

A function to extract the bases 3' upstream and 5' downstream of the base substitution types.

Usage

1
type_context(vcf, ref_genome, extension = 1)

Arguments

vcf

A CollapsedVCF object

ref_genome

Reference genome

extension

The number of bases, that's extracted upstream and downstream of the base substitutions. (Default: 1).

Value

Mutation types and context character vectors in a named list

See Also

read_vcfs_as_granges, mut_context

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
## See the 'read_vcfs_as_granges()' example for how we obtained the
## following data:
vcfs <- readRDS(system.file("states/read_vcfs_as_granges_output.rds",
  package = "MutationalPatterns"
))

## Load the corresponding reference genome.
ref_genome <- "BSgenome.Hsapiens.UCSC.hg19"
library(ref_genome, character.only = TRUE)

## Get type context
type_context <- type_context(vcfs[[1]], ref_genome)

## Get larger type context
type_context_larger <- type_context(vcfs[[1]], ref_genome, extension = 2)

MutationalPatterns documentation built on Nov. 14, 2020, 2:03 a.m.