## ----style, eval=TRUE, echo=FALSE, results="asis"--------------------------
BiocStyle::latex()

## ----echo=FALSE, results='hide'--------------------------------
set.seed(1)
options(width=65)
library(knitr)
opts_chunk$set(background = "#FFFFFF", comment = NA, fig.path = "")

## ----load_packages,results='hide'------------------------------
library(VariantAnnotation)
library(GenomicRanges)

## ----load_data-------------------------------------------------
data(package = "COSMIC.67")
data(cosmic_67, package = "COSMIC.67")

## ----overlap_with_tp53-----------------------------------------
tp53_range = GRanges("17", IRanges(7565097, 7590856))
vcf_path = system.file("vcf", "cosmic_67.vcf.gz", package = "COSMIC.67")
cosmic_tp53 = readVcf(vcf_path, genome = "GRCh37", ScanVcfParam(which = tp53_range))
cosmic_tp53

## ----CGC-------------------------------------------------------
data(cgc_67, package = "COSMIC.67")
head(cgc_67)

## ----echo=FALSE------------------------------------------------
sessionInfo()