plot_specificity Plot Cas9 specificity of a list of multiple gRNA offtargets.

Usage

plot_specificity(x, ref = 1)

Arguments

x: A named list of GRanges object containing the results of a breakinspectoR analysis. Possibly filtered (eg. qval < .01).
ref: Which experiment from `x` should be taken as reference? either a number or a name. Defaults to the first experiment in `x`. Can be NULL, for no ref.

Value

A ggplot object.

Examples

if (FALSE) { # \dontrun{
  offtargets <- breakinspectoR(
    target   =system.file("extdata/vegfa.chr6.bed.gz", package="breakinspectoR"),
    nontarget=system.file("extdata/nontarget.chr6.bed.gz", package="breakinspectoR"),
    guide    ="GACCCCCTCCACCCCGCCTC",
    PAM      =c("NGG", "NAG"),
    bsgenome ="BSgenome.Hsapiens.UCSC.hg38",
    cutsiteFromPAM=3
  )
} # }
data(breakinspectoR_examples)

# simulate 2 different experiments by picking 25 random offtargets
spcas9 <- c(offtargets[offtargets$mismatches == 0], sample(offtargets, 25))
lz3    <- c(offtargets[offtargets$mismatches == 0], sample(offtargets, 25))
plot_specificity(list(SpCas9=spcas9, LZ3=lz3))