read ENCODE rrbs bed file

read.rrbs.bed.R

library(methylKit)
obj = read("wgEncodeHaibMethylRrbsA549Dm002p7dHaibSitesRep1.bed.gz",
                 sample.id = "test", assembly = "hg19", header = FALSE, 
                 context = "CpG", resolution = "base", 
                 pipeline = list(fraction = FALSE, chr.col = 1, start.col = 3, 
                                  end.col = 3,coverage.col = 5, strand.col = 6, 
                                  freqC.col = 11))