KataegisPortal is an R package that provides utilities to detect and visualize the Kataegis events from single nucleotide variants (SNVs).
library(devtools)
install_github("MeichunCai/KataegisPortal")
# Load packages
library(KataegisPortal)
library(BSgenome)
library(BSgenome.Hsapiens.UCSC.hg19)
# Read in snv data
mutData <- system.file("examples", "mutData.txt", package="KataegisPortal")
mutData <- read.table(mutData, header = TRUE,sep = "\t",as.is = TRUE)
head(mutData)
# Covert to KataegisPortal input
mutSNP = mutSNP.input(mut.data = mutData,
chr = "chr",
pos = "pos",
ref = "ref",
alt = "alt",
build = "hg19")
# Identify potential kataegis events
katPoint(mutSNP)
# Identify potential kataegis events with annonation
library(ChIPseeker)
library(TxDb.Hsapiens.UCSC.hg19.knownGene)
katPoint(mutSNP,txdb = TxDb.Hsapiens.UCSC.hg19.knownGene)
# Plot intermutation distance to display kataegis or hyper mutated genomic regions
mutDis.plot(plot.data = mutSNP, sample="Test")
mutDis.plot(plot.data = mutSNP, sample="Test", chr = "chr2")
# Plot of flanking sequence of all C>X mutations or C>X mutations within a chromosome/chromosome arm
baseSpe.plot(plot.data = mutSNP, sample="Test")
baseSpe.plot(plot.data = mutSNP, sample="Test", chr = "chr2",arm="q")