Skip to content
Snippets Groups Projects
Commit 34036959 authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Format multiqc

parent 24271608
Branches
Tags
4 merge requests!53Develop,!47Resolve "Add MultiQC",!46Resolve "Add MultiQC",!44Add MultiQC to output. Add version and reference material. fix README
Pipeline #4197 passed with stages
in 27 minutes and 9 seconds
......@@ -5,7 +5,7 @@ custom_logo_title: 'Bioinformatics Core Facility'
report_header_info:
- Contact E-mail: 'bicf@utsouthwestern.edu'
- Application Type: 'CellRanger_Count'
- Application Type: 'cellranger_count'
- Department: 'Bioinformatic Core Facility, Department of Bioinformatics'
......@@ -56,7 +56,7 @@ table_columns_visible:
thousandsSep_format: ''
report_section_order:
Software_Versions:
software_versions:
order: -1100
Software_References:
software_references:
order: -1200
......@@ -61,20 +61,18 @@ references = params.references
process checkDesignFile {
tag "$name"
publishDir "$outDir/misc/${task.process}/$name", mode: 'copy'
module 'python/3.6.1-2-anaconda'
input:
file designLocation
file fastqList
output:
file("design.checked.csv") into designPaths
script:
"""
hostname
ulimit -a
......@@ -118,14 +116,12 @@ process count211 {
module 'cellranger/2.1.1'
input:
set sample, file("${sample}_S1_L00?_R1_001.fastq.gz"), file("${sample}_S1_L00?_R2_001.fastq.gz") from samples211
file ref from refLocation211.first()
expectCells211
forceCells211
output:
file("**/outs/**") into outPaths211
file("*_metrics_summary.tsv") into metricsSummary211
......@@ -160,7 +156,6 @@ process count301 {
module 'cellranger/3.0.1'
input:
set sample, file("${sample}_S1_L00?_R1_001.fastq.gz"), file("${sample}_S1_L00?_R2_001.fastq.gz") from samples301
file ref from refLocation301.first()
expectCells301
......@@ -168,7 +163,6 @@ process count301 {
chemistryParam301
output:
file("**/outs/**") into outPaths301
file("*_metrics_summary.tsv") into metricsSummary301
......@@ -178,7 +172,7 @@ process count301 {
script:
if (forceCells301 == 0){
"""
hostname
hostname
ulimit -a
bash "$baseDir/scripts/filename_check.sh" -r "$ref"
cellranger count --id="$sample" --transcriptome="./$ref" --fastqs=. --sample="$sample" --expect-cells=$expectCells301 --chemistry="$chemistryParam301"
......@@ -186,7 +180,7 @@ process count301 {
"""
} else {
"""
hostname
hostname
ulimit -a
bash "$baseDir/scripts/filename_check.sh" -r "$ref"
cellranger count --id="$sample" --transcriptome="./$ref" --fastqs=. --sample="$sample" --force-cells=$forceCells301 --chemistry="$chemistryParam301"
......@@ -203,7 +197,6 @@ process count302 {
module 'cellranger/3.0.2'
input:
set sample, file("${sample}_S?_L001_R1_001.fastq.gz"), file("${sample}_S?_L001_R2_001.fastq.gz") from samples302
file ref from refLocation302.first()
expectCells302
......@@ -211,7 +204,6 @@ process count302 {
chemistryParam302
output:
file("**/outs/**") into outPaths302
file("*_metrics_summary.tsv") into metricsSummary302
......@@ -247,11 +239,9 @@ process versions {
input:
output:
file("*.yaml") into yamlPaths
script:
"""
hostname
ulimit -a
......@@ -271,16 +261,13 @@ process multiqc {
module 'multiqc/1.7'
input:
file ('*') from metricsSummary.collect()
file yamlPaths
output:
file "*" into mqcPaths
script:
"""
hostname
ulimit -a
......@@ -288,4 +275,4 @@ process multiqc {
sed -i '1s/^.*\tE/Sample\tE/' metrics_summary_mqc.tsv
multiqc -c $multiqcConf .
"""
}
}
\ No newline at end of file
......@@ -51,8 +51,9 @@ def main():
out_filename = output + '_mqc.yaml'
# Header for HTML
print('''
id: 'Software References'
print(
'''
id: 'software_references'
section_name: 'Software References'
description: 'This section describes references for the tools used.'
plot_type: 'html'
......
File mode changed from 100644 to 100755
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment