RNA-seq workflow using STAR and DESeq2
Just a quick release with some fixes that have accumulated, so that deployment via snakedeploy
and the snakemake workflow catalog will also include those:
This is a major improvement of this workflow, bringing it up to date and extending its functionality. Many thanks to all of the people who reported problems and suggested improvements over the last 1.5 years. I'll try to name all here, but please ping me if I forgot anyone. There was input and contributions from (in no particular order): @isaacvock, @kilpert, @dlaehnemann, @jonathandmoore, @cbp44, @aryazand, @iwv2, @FelixMoelder, @johanneskoester, @svenrahmann
Major change in workflow configuration (/ configurability):
DESeq2
's flexibility in its model setup and does not make use of a particular name for the variable_of_interest
any more. Instead, there is a new, fully explained, and reasonably intuitive configuration setup in the config.yaml
for this, including linkouts to the relevant docs and an explanation in the snakemake workflow catalog docs (that have been improved for unchanged features, as well). Also, this setup and the workflow in general is now more comprehensively tested.Also lots of smaller fixes:
{sample}_{unit}
vs. {sample)-{unit}
that we had in some places, PR #61)others:
to extra:
(originally suggested by @kilpert: https://github.com/snakemake-workflows/rna-seq-star-deseq2/commit/f816550598c29ed2840e21b31a592ce87a7ebcfa)And a new feature addition (hasn't kicked in yet, this will probably come with the next PR merge):
Format workflow according to Python and Snakemake style guides (black and snakefmt).
Fixed naming of snakemake-workflow-catalog configuration file.
Added QC.