Closed rcavalcante closed 9 years ago
CT SNP GRanges object has start = end and is 1-based. Input data from methylKit and from Bismark is 1-based. However, findOverlaps(...) is still not finding all the locations of C > T SNPs in the methylation data.
Realized that I neglected to pass filterSNPs parameter from readBismarkData(...) and methylSigReadData(...) to their helper functions.
CT SNP GRanges object has start = end and is 1-based. Input data from methylKit and from Bismark is 1-based. However, findOverlaps(...) is still not finding all the locations of C > T SNPs in the methylation data.