Skip to content

Instantly share code, notes, and snippets.

@mbk0asis

mbk0asis/all CpG locations in mm10

Last active August 16, 2018 04:16
Show Gist options
  • Save mbk0asis/702caf583a534f7ef531db8cc68929c3 to your computer and use it in GitHub Desktop.
Save mbk0asis/702caf583a534f7ef531db8cc68929c3 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
all CpG locations in mm10
source("https://bioconductor.org/biocLite.R")
biocLite("BSgenome.Mmusculus.UCSC.mm10")
setwd("~/BIO2/MBDseq_Mouse_Muscle_Tcell_JungHo")
library(BSgenome.Mmusculus.UCSC.mm10)
chrs <- names(Mmusculus)[1:21]
CGs <- lapply(chrs, function(x) start(matchPattern("CG", Mmusculus[[x]])))
names(CGs) <- chrs
str(CGs)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment