contextseqs.Rd
Get [-3, +3] contextseqs for given cas9ranges
contextseqs(cas9ranges)
cas9ranges | GenomicRanges::GRanges |
---|
require(magrittr) bedfile <- system.file('extdata/SRF_sites.bed', package = 'multicrispr') bsgenome <- BSgenome.Mmusculus.UCSC.mm10::Mmusculus cas9ranges <- read_bed(bedfile, bsgenome) %>% flank_fourways() %>% find_cas9ranges()#>#>#>#>#>#>#>#>#>#>#>#>#>#> [1] "CTCTAGAGTCACAGAACTTATGG" "AATTTATTATGACTTACAGTCGG" #> [3] "TTCAGTTGCAGCTGTGAAAATGG"cas9ranges[1:3] %>% contextseqs()#> [1] "GGTTCTCTAGAGTCACAGAACTTATGGACA" "AAGGAATTTATTATGACTTACAGTCGGCAG" #> [3] "ATTGTTCAGTTGCAGCTGTGAAAATGGAAG"