Find target cas9 ranges with no offtargets

find_specific_cas9ranges(targetranges, mismatch = 2, verbose = TRUE)

Arguments

targetranges

GenomicRanges::GRanges

mismatch

max number of mismatches to consider

verbose

logical(1)

Value

GenomicRanges::GRanges mcols(GRanges) contains sequences and match counts: matches0 = perfect match counts matches1 = single mismatch counts matches2 = double mismatch counts

Examples

# NOT RUN {
   require(magrittr)
   bedfile <- system.file('extdata/SRF_sites.bed', package = 'multicrispr')
   bsgenome <- BSgenome.Mmusculus.UCSC.mm10::Mmusculus
   targetranges <- read_bed(bedfile, bsgenome)[1:10] %>% flank_fourways()
   targetranges %>% find_specific_cas9ranges(mismatch=0)
# }