uniting methylRawDB objects with regional resolution and setting destrand=TRUE deletes the original objects

[alexg9010]

methylKit/R/methylDBFunctions.R

Line 747 in 4735b6f

on.exit(unlink(c(getDBPath(object),paste0(getDBPath(object),".tbi"))),add = TRUE)

[alexg9010]

library(methylKit)
data(methylKit)

mdb <- makeMethylDB(methylRawList.obj)     

# works for "base" -resolution
unite(mdb,destrand=TRUE)
mdb # is fine

# breaks for "region" -resolution
mdbt <- tileMethylCounts(mdb) 
# works
unite(mdbt)    

# deletes tabix files of mdbt
unite(mdbt,destrand=TRUE)    

mdbt # breaks !!