Let’s look at our snakemake “trim” rule from earlier:

rule trim:

input:

”reads/sample1.fastq”

output:

”trimmed_reads/sample1-trimmed.fastq”

shell:

cutadapt -A TCCGGGTS -o {output} {input}

Here we specified our inputs/outputs and our shell command.