I am currently experimenting snakemake for pipelines. We will have paired samples when processing data. e.g. ChIP-seq (IP vs Input), DNA-seq (Normal vs Tumor). How to specify them in snakemake?
The trick is to use a python dictionary to pair samples and a function to define Input files.
I will show an toy example with ChIP-seq processing.
This blog by Tommy Tang is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.
