If mode is set to `regions`, annotation file must be provided.

Output is JSON for genome-wide mode and TSV for regions mode.

Examples:
For genome-wide mode:
bsxplorer stats --output stats.json report.ipc

For regions mode:
bsxplorer stats --output stats.tsv --threads 12 --mode regions --format gff -a genomic.gff report.ipc
