find_cas9ranges.Rd
Find cas9 ranges in targetranges
find_cas9ranges(targetranges, verbose = TRUE)
targetranges | GenomicRanges::GRanges |
---|---|
verbose | logical(1) |
GenomicRanges::GRanges
require(magrittr) bedfile <- system.file('extdata/SRF_sites.bed', package='multicrispr') bsgenome <- BSgenome.Mmusculus.UCSC.mm10::Mmusculus targetranges <- read_bed(bedfile, bsgenome) %>% flank_fourways()#>#>#>#>#>#>#>#>#>#>find_cas9ranges(targetranges)#>#>#>#> GRanges object with 33518 ranges and 1 metadata column: #> seqnames ranges strand | seqs #> <Rle> <IRanges> <Rle> | <character> #> [1] chr1 4712436-4712458 + | CTCTAGAGTCACAGAACTTATGG #> [2] chr1 4712480-4712502 + | AATTTATTATGACTTACAGTCGG #> [3] chr1 4712524-4712546 + | TTCAGTTGCAGCTGTGAAAATGG #> [4] chr1 4712566-4712588 + | TTACTCAGTCTCACGCAAGCAGG #> [5] chr1 4712650-4712672 + | GAAGGTGTAGCCCAGATTAAAGG #> ... ... ... ... . ... #> [33514] chrY 89126327-89126349 - | AGCGTGGACACGCCCACCCTCGG #> [33515] chrY 89126684-89126706 - | CTTCCAGGCCCTAATAGCCGTGG #> [33516] chrY 89126655-89126677 - | TCTGGGCATAAGCTCCCCAATGG #> [33517] chrY 89126584-89126606 - | AAAGTCCAATCTCTCTGCTCTGG #> [33518] chrY 89126559-89126581 - | TCAAAGCAGCTGCATTTGCTCGG #> ------- #> seqinfo: 66 sequences (1 circular) from mm10 genome