Questions tagged [snakemake]
Snakemake is a workflow management system with a Python-style specification language.
snakemake
1,798
questions
0
votes
1
answer
32
views
How Do I Combine Parameters Conditionally?
Suppose I have the following input FASTQ files:
$ ls -1 inputs/
CM0009619-ABB_L01_Read1_Sample_Library_XY102.fastq.gz
CM0009619-ABB_L01_Read1_Sample_Library_XY110.fastq.gz
CM0009619-...
1
vote
1
answer
36
views
Snakemake running as an AWS Batch or an AWS Fargate task raises MissingInputException on the inputs stored on a S3 bucket
We have a Dockerized Snakemake pipeline with the input data stored on a S3 bucket snakemake-bucket:
Snakefile:
rule bwa_map:
input:
"data/genome.fa"
output:
"...
0
votes
0
answers
14
views
Missing Output Exception Error in Snakemake with output directory modified
I'm trying to run a workflow on snakemake. I have to automate a couple of steps which are depending all on python scripts or pipelines already made. My rule Gene_flow_between _species has to run once ...
0
votes
0
answers
12
views
Missing Output Exception Error in Snakemake with output directory
I'm trying to run a workflow on snakemake. I have to automate a couple of steps which are depending all on python scripts or pipelines already made. My rule Gene_flow_between _species has to run once ...
0
votes
1
answer
18
views
Multiple variables in input and output using wildcards snakemake
rule all:
input:
expand("tissues/{id}/atac_seq/{srr}/{srr}_S1_L001_R1_001.fastq.gz", id =IDs, srr = df[df['library_type']=='Chromatin Accessibility']['Library']),
expand("...
0
votes
1
answer
36
views
How to log time duration (for each job and total) for a snakemake workflow?
If I run a snakemake workflow, the start time for each job is logged to the console, for example
[Fri Jul 12 10:39:15 2024]
How can I adapt the logging to also show the duration for each job and the ...
1
vote
1
answer
28
views
How to run snakemake on AWS now that Tibanna is not anymore available on Snakemake 8
Currently I am running snakemake pipelines on AWS using an older version of Snakemake (7.30.0). I would like to update to a more recent version of Snakemake. Unfortunately it seems that Tibanna is not ...
1
vote
1
answer
21
views
How to enable "output" to detect variables in "run" in snakemake
I am trying to make directories for the different "IDs" I have stored in the .csv file (they are stored under column "ID"). However, snakemake doesn't seem to be able to detect &...
0
votes
0
answers
19
views
Snakemake rule generating symlink in Python script reruns always with "Code has changed" despite no code changes
I have a snakemake rule that generates symbolic links in a target folder. This is done within a Python script and the script directive in the rule. This rule is always rerun with reason: Code has ...
1
vote
2
answers
30
views
Can you order the execution of all expanded jobs in single rule?
I have a snakemake pipeline that runs a genetic analysis. The genome has been split into many 'regions'. These regions can run in parallel, and therefore I've used expand(), and they all run as ...
1
vote
1
answer
34
views
config file taking it as literal string snakemake
shell:
"""
cellranger-arc count --id={wildcards.atac_srr}_{wildcards.rna_srr} \
--reference={config['path']['ref_genome']} \
--...
0
votes
1
answer
37
views
How to deal with paired tumor-normal samples in Snakemake?
I am converting a bash pipeline to process paired whole-exome sequenced tumor-normal samples into a SnakeMake workflow.
Paired samples are listed in my config file, as follows:
sample_list:
- sample:...
0
votes
0
answers
22
views
How to Set Exception Rules for Slurm Executor in Snakemake?
I am currently working with a Slurm workload manager to manage and schedule jobs on our cluster. While most of my tasks are submitted and managed by Slurm, I have a few specific tasks that I would ...
0
votes
1
answer
32
views
Snakemake, how to reference input in output of the same rule
My logic is probably flawed here, and I need your help to get the best possible approach to solving this:
I have a daisy-chained workflow where I annotate a .vcf file in multiple steps, and each step ...
0
votes
1
answer
92
views
How to debug scripts for PyPSA-EUR?
I would like to set some VisualStudioCode or PyCharm breakpoints for the script
https://github.com/PyPSA/pypsa-eur/blob/master/scripts/prepare_sector_network.py
and then run and debug it to better ...