WIP

gene_stats()