plot_gene function

Plot gene expression