Skip to content
GitLab
Menu
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Klinische Genetica
capture-lumc
hutspot
Commits
6ba8b697
Commit
6ba8b697
authored
Mar 30, 2018
by
Sander Bollen
Browse files
add stats to tsv in pipeline
parent
d74bc3be
Changes
1
Hide whitespace changes
Inline
Side-by-side
Snakefile
View file @
6ba8b697
...
...
@@ -28,6 +28,7 @@ colpy = fsrc_dir("src", "collect_stats.py")
mpy = fsrc_dir("src", "merge_stats.py")
seq = fsrc_dir("src", "seqtk.sh")
fqpy = fsrc_dir("src", "fastqc_stats.py")
tsvpy = fsrc_dir("src", "stats_to_tsv.py")
if FASTQ_COUNT is None:
fqc = "python {0}".format(fsrc_dir("src", "fastq-count.py"))
...
...
@@ -585,13 +586,24 @@ rule merge_stats:
shell: "python {input.mpy} --vcfstats {input.vstat} {input.cols} > {output.stats}"
rule stats_tsv:
"""Convert stats.json to tsv"""
input:
stats=out_path("stats.json"),
sc=tsvpy
output:
stats=out_path("stats.tsv")
conda: "envs/collectstats.yml"
shell: "python {input.sc} -i {input.stats} > {output.stats}"
rule multiqc:
"""
Create multiQC report
Depends on stats.
json
to forcefully run at end of pipeline
Depends on stats.
tsv
to forcefully run at end of pipeline
"""
input:
stats=out_path("stats.
json
")
stats=out_path("stats.
tsv
")
params:
odir=out_path("."),
rdir=out_path("multiqc_report")
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment