find_specific_cas9ranges.Rd
Find target cas9 ranges with no offtargets
find_specific_cas9ranges(targetranges, mismatch = 2, verbose = TRUE)
targetranges | GenomicRanges::GRanges |
---|---|
mismatch | max number of mismatches to consider |
verbose | logical(1) |
GenomicRanges::GRanges mcols(GRanges) contains sequences and match counts: matches0 = perfect match counts matches1 = single mismatch counts matches2 = double mismatch counts
# 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) # }