Snakemake Singularity with Local Resources / Questions about Snakemake with --use-singularity

user_14480405 0
Viewed 422

I'm starting to experiment with using containers with Snakemake, and I have a question about what needs to be pre-build into the container and what doesn't. For example:

I want to run a python script (stored in workflow_root/scripts/myScript.py, for example) in a container with a pipe in from another program. Do I need to build the python script into the container, declare it as an input file, or is that accessible from within the container (and how do I point to it)? My current rule looks something like:

rule myRule:
    params:
        sample = get_sample,
        basePath = sys.path[0]
    input:
        in1=get_in1,
        in2=get_in2
    output:
        out1 = "{runPath}/{sample}_read1_dcs.fq.gz",
        out2 = "{runPath}/{sample}_read2_dcs.fq.gz"
    priority: 50
    conda:
       "envs/myEnv.yaml"
    log:
        "{runPath}/logs/{sample}_myRule.log"
    shell:
        """
        set -e
        set -o pipefail
        set -x
        {{
        picard FastqToSam \
        F1={input.in1} \
        F2={input.in2} \
        O=/dev/stdout \
        SM={params.sample} \
        TMP_DIR=picardTempDir \
        SORT_ORDER=unsorted \
        | python3 {params.basePath}/scripts/myScript.py \
        --input /dev/stdin \
        --prefix {wildcards.sample}
        }} 2>&1 | tee -a {log}
        """

I want to run bwa, where I have a sizable user-provided reference that I need to use. Can I do this, or would I need to build that reference into the container? (I'd also like to use ensemble-VEP, which has its own sizable reference database to deal with).

I suppose what my question boils down to is: what files / locations are mounted to the container by Snakemake, and where do I find them when I'm writing rules involving shell commands? The documentation doesn't seem to be very clear on this, and it would be nice to be able to figure it out without having to do a bunch of experimentation.

singularity-container snakemake docker bioinformatics
2 Answers

I will share how I use snakemake with singularity and conda. This setup has worked very well for me for over an year now. This setup may or may not suit your purposes; so feel free to ask questions.

  • Snakemake workflow along with necessary scripts, configs and docs are in a separate git repo. Data would not be stored here.
  • Singularity container is defined at global level via singularity: directive in the workflow. I don't manually build the container; snakemake does that.
  • Conda environments are defined per rule via conda: directive. In complex projects, sometimes this may be >15-20 separate conda environments. Snakemake builds these conda environments inside the singularity container.
  • Data are kept separate from the source code. Instead the workflow config file (typiecally defined via configfile: directive) contains their path info.

An example Snakefile would look like this:

# config for the workflow
configfile: "configs/workflow_configs.yaml"

# singularity image to use
singularity: "docker://continuumio/miniconda3:4.7.12"

rule all:
   input:
       .....

rule some_job:
    input:
        .....
    output:
        .....
    conda: 
        "path/to/conda_env.yaml"
    shell:
        "....."

This setup works as long as all the tools needed for the workflow are available via conda. It is also possible to override global singularity container for a specific rule and use a different singularity container or no container at all.

Also, I don't build my own singularity containers and instead I just use the generic ones available from docker-hub. I don't have privileges to build my own singularity container in the HPC systems I work with, anyway. So this setup removes the hassle of building the image somewhere and then moving it to the HPC environment.

user_39982520

Everything under the working directory is available and accessible as expected. Things referenced through absolute paths, however, will require mounting.

So, if your script is under scripts/my_script.py, that will be available as is. But, if you use, for example, a temporary directory like /scratch/user42, then you will need to mount that with something like

snakemake --singularity-args '-B /scratch/user42:/scratch/user42'

Personally, I specify my mounts in the config.yaml for my cluster profile, so I don't ever have to bother with the extra CLI arguments.

user_5709180
Related