I'm using R packages to analysis my BS-sequencing data. I first used methylKit to find base-wise differential methylation. This result is with strand information for each base, looks like:
chr start end strand pvalue qvalue meth.diff
1 chr1 3000827 3000827 + 0.3666667 1 -9.090909
2 chr1 3000828 3000828 - 0.3264993 1 -8.029690
Then I applied eDMR to this data to find differential methylated region. But in the result, the strand information is all replaced with "*".
Is there a way to keep strand information?
Thanks in advance.