Call in `atacseq.characterize_regions_function` used "genome" instead of "genome_assembly".
Call in
atacseq.characterize_regions_functionused "genome" instead of "genome_assembly".