Should I remove incomplete bisulfite conversion in the BSSeeker2 call step?
0
0
Entering edit mode
5.2 years ago
colin.kern ★ 1.1k

I have a question about these two arguments for the BSSeeker2 methylation calling step:

-x, --rm-SX           Removed reads with tag 'XS:i:1', which would be
                      considered as not fully converted by bisulfite
                      treatment [Default: False]
--rm-CCGG             Removed sites located in CCGG, avoiding the bias
                      introduced by artificial DNA methylation status
                      'XS:i:1', which would be considered as not fully
                      converted by bisulfite treatment [Default: False]

Both of these seem to be removing false positives from the results, which seem like a generally good thing to do, but their default value is "False". Usually the defaults for a program are calibrated to the most general use case, but I'm not sure if that's the case here. Is it standard practice to use these arguments? Why would I want to not use these?

I have RRBS data, not WGBS, if that makes any difference.

methylation rrbs • 882 views
ADD COMMENT

Login before adding your answer.

Traffic: 2063 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6