gene_stats function

WIP