Performance comparison for HISAT-3N, Bismark, BS-Seeker2, and BSMAP on 78 million real whole-genome paired-end BS-seq reads