Add phenotype order for plotting lavalamps
Add xlabel, ylabel to plot_modaltlies_reduced, and only set ylim of plot_feature for violinplots, not the reduced plot