From e95c85b993555c1867651ab7d3601e3c422b3753 Mon Sep 17 00:00:00 2001
From: Blaise Li <blaise.li__git@nsup.org>
Date: Thu, 16 May 2019 17:08:06 +0200
Subject: [PATCH] Configurable gene lists for PRO-seq.

---
 PRO-seq/PRO-seq.snakefile | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/PRO-seq/PRO-seq.snakefile b/PRO-seq/PRO-seq.snakefile
index 2d31613..c992b56 100644
--- a/PRO-seq/PRO-seq.snakefile
+++ b/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 #
-- 
GitLab