Configurable gene lists for PRO-seq.

PRO-seq/PRO-seq.snakefile

@@ -164,6 +164,7 @@ DE_BIOTYPES = [biotype for biotype in LFC_RANGE.keys() if biotype in COUNT_BIOTY
 #METAGENE_BIOTYPES = ["protein_coding"]
 #METAGENE_BIOTYPES = ["protein_coding", "protein_coding_5UTR", "protein_coding_CDS", "protein_coding_3UTR"]
 METAGENE_BIOTYPES = [biotype for biotype in ["protein_coding", "protein_coding_5UTR", "protein_coding_CDS", "protein_coding_3UTR"] if biotype in COUNT_BIOTYPES]
+# default id lists for MA plots
 ID_LISTS = [
     "lfc_statuses",
     "germline_specific",
@@ -171,6 +172,7 @@ ID_LISTS = [
     "spermatogenic_Ortiz_2014", "oogenic_Ortiz_2014",
     "piRNA_dependent_prot_si_22G_down4_top200", "piRNA_dependent_prot_si_22G_down4",
     "csr1_prot_si_supertargets_common"]
+ID_LISTS = config.get("maplot_gene_lists", ID_LISTS)
 aligner = config["aligner"]
 ########################
 # Genome configuration #