rework read_utils.split_reads to use numeric file indexes and work on gzipped inputs and outputs