add an option to output only sites that have nonzero coverage in both samples

by default, the program outputs all sites regardless of coverage.
the default behavior is preferred as the dmr program needs all sites to correcntly count the number of CpG sites in a DMR