# Snakemake Fastqc: "Multiple run or shell keywords in rule run_fastqc."

I am trying to check the quality of RNA-Seq data from Illumina using fastqc in snakemake in a conda environment. I get the error "Multiple run or shell keywords in rule run_fastqc".

Snakefile:

rule run_fastqc:
input:
"/[RNA fastq file path]/well01_S001_L001_R1_001.fastq.gz"
output:
"./[out path]/well01_S001_L001_R1_001.fastqc.html",
"./[out path]/well01_S001_L001_R1_001.fastqc.zip"
wrapper:
"0.79.0/bio/fastqc"
shell:
"fastqc {input} --outdir=./[out path]/"


I've also tried:

...
shell:
"fastqc {input} -o ./[out path]/"
...
shell:
"fastqc -o ./[out path]/ {input}"
...
shell:
"fastqc {input}"
"fastqc -o ./[out path]/"
...
shell:
"fastqc {input} > -o ./[out path]/"


Code:

snakemake --cores 1 well01_S001_L001_R1_001.fastqc.html well01_S001_L001_R1_001.fastqc.zip


Output:

SyntaxError in line 9 of /[path to Snakefile]/Snakefile:
Multiple run or shell keywords in rule run_fastqc. (Snakefile, line 9)


When I run fastqc alone (not in snakemake), it works and I get the two output files. How can I get snakemake to work?

Thank you for your time and help!

## 1 Answer

From a quick look at the manual, it looks like you are supposed to use just the wrapper, not combine the wrapper with a shell. Can you show me an example of such combined usage?

The error seems to suggest that the wrapper already has a shell directive that you shell conflicts with.

• Thank you! In Output, I also had to remove the "./[out path]/" before the output file name, I guess the "--outdir" in Shell is sufficient. I got the files but still received the error: "Job Missing files after 5 seconds: well01_S001_L001_R1_001.fastqc.html well01_S001_L001_R1_001.fastqc.zip This might be due to filesystem latency. If that is the case, consider to increase the wait time with --latency-wait. Job id: 0 completed successfully, but some output files are missing. 0" and I also waited for 60s and got the same error. Nov 11 '21 at 21:06
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 0, "mathjax_asciimath": 1, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0, "math_score": 0.49565795063972473, "perplexity": 6659.434638879001}, "config": {"markdown_headings": true, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": false}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "plain"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "warc_path": "s3://commoncrawl/crawl-data/CC-MAIN-2022-05/segments/1642320304961.89/warc/CC-MAIN-20220126192506-20220126222506-00228.warc.gz"}