Take a gene name or ID, search the counts matrix for all clusters that fall within that gene's boundaries (found from the exon table).

make_gene_plot(
  gene_name,
  clusterID = NULL,
  cluster_list = NULL,
  introns = NULL,
  introns_to_plot = NULL,
  counts,
  exons_table = NULL,
  len = 500,
  min_exon_length = 0.5,
  main_title = NA,
  snp_pos = NA,
  snp = NA,
  summary_func = colSums,
  legend_title = "Mean counts",
  debug = F
)