```r
# meth_unite_reshaped <- melt(data=meth_unite, id=c("chr", "start", "end", "strand"), value.name = "count") %>%
#  tidyr::separate(col = "variable" , into = c("variable", "sample"), sep = "(?<=[a-zA-Z])\\s*(?=[0-9])") %>%
#   dcast(chr+start+end+strand+sample ~  variable) %>%
#   drop_na(coverage) %>% 
#   mutate(percMeth = 100*(numCs/coverage)) %>% 
#   mutate(sample=as.numeric(sample))
# 
# save(meth_unite_reshaped, file = "meth_unite_reshaped") 
# save(meth_unite_reshaped, file = "../results/methylkit/meth_unite_reshaped") 
```