MultiQC changes

assets/multiqc_config.yml

@@ -45,16 +45,193 @@ module_order:
         - "*.final.out"
   - custom_content
 
+# Other MultiQC config stuff here
+custom_data:
+  mapping:
+    parent_id: mapping
+    parent_name: "Mapping"
+    file_format: "tsv"
+    section_name: "Mapping"
+    description: "The mapping metrics for each experiment"
+    plot_type: "bargraph"
+  dedup_reads:
+    parent_id: dedup
+    parent_name: "Deduplication"
+    file_format: "tsv"
+    section_name: "Reads"
+    description: "The number of reads before and after PCR deduplication for each experiment"
+    plot_type: "bargraph"
+    pconfig:
+      ylab: "Count"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+  dedup_ratio:
+    parent_id: dedup
+    parent_name: "Deduplication"
+    file_format: "tsv"
+    section_name: "Ratio"
+    description: "The PCR deduplication ratio for each experiment"
+    plot_type: "bargraph"
+    pconfig:
+      ylab: "Ratio"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+  dedup_mean_umis:
+    parent_id: dedup
+    parent_name: "Deduplication"
+    file_format: "tsv"
+    section_name: "Mean UMIs"
+    description: "Mean number of unique UMIs per position for each experiment"
+    plot_type: "bargraph"
+    pconfig:
+      ylab: "Mean number"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+  crosslinks_counts:
+    parent_id: crosslinks
+    parent_name: "Crosslinks"
+    file_format: "tsv"
+    section_name: "Counts"
+    description: "The number of crosslinks or crosslink sites for each experiment"
+    plot_type: "bargraph"
+    pconfig:
+      ylab: "Count"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+  crosslinks_ratio:
+    parent_id: crosslinks
+    parent_name: "Crosslinks"
+    file_format: "tsv"
+    section_name: "Ratios"
+    description: "The ratio of number of cDNA mapping to crosslink positions for each experiment"
+    #plot_type: 'bargraph'
+    pconfig:
+      ylab: "Count"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+      tt_decimals: 2
+  peaks_counts:
+    parent_id: peaks
+    parent_name: "Peaks"
+    file_format: "tsv"
+    section_name: "Counts"
+    description: "The total number of peaks called by each peak caller"
+    plot_type: "bargraph"
+    pconfig:
+      ylab: "Number of peaks"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+  xlinks_in_peaks:
+    parent_id: peaks
+    parent_name: "Peaks"
+    file_format: "tsv"
+    section_name: "Crosslinks positions in peaks"
+    description: "The total percentage of crosslinks within peaks for each peak caller"
+    #plot_type: 'bargraph'
+    pconfig:
+      ylab: "Percentage of crosslinks"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+      tt_decimals: 2
+      tt_suffix: "%"
+  xlinksites_in_peaks:
+    parent_id: peaks
+    parent_name: "Peaks"
+    file_format: "tsv"
+    section_name: "Crosslinks positions in peaks"
+    description: "The total percentage of crosslink sites within peaks for each peak caller"
+    #plot_type: 'bargraph'
+    pconfig:
+      ylab: "Percentage of crosslink sites"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+      tt_decimals: 2
+      tt_suffix: "%"
+  peaks_xlinksite_coverage:
+    parent_id: peaks
+    parent_name: "Peaks"
+    file_format: "tsv"
+    section_name: "Peak-crosslink coverage"
+    description: "The total percentage of nucleotides within peaks covered by a crosslink site"
+    plot_type: "bargraph"
+    pconfig:
+      ylab: "Percentage of nucleotides within peaks"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+      tt_decimals: 2
+      tt_suffix: "%"
+  summary_type:
+    parent_id: Summary
+    parent_name: "Summary"
+    file_format: "tsv"
+    section_name: "Percentage of cDNA premap"
+    description: "The total percentage of cDNA summary mapped"
+    #plot_type: 'bargraph'
+    pconfig:
+      ylab: "Type"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+  summary_subtype:
+    parent_id: Summary
+    parent_name: "Summary"
+    file_format: "tsv"
+    section_name: "Percentage of cDNA premap subtypes"
+    description: "The total percentage of cDNA subtypes mapped"
+    #plot_type: 'bargraph'
+    pconfig:
+      ylab: "Type"
+      #stacking: False
+      cpswitch: False
+      tt_percentages: False
+
+sp:
+  mapping:
+    fn: "mapping.tsv"
+  dedup_reads:
+    fn: "dedup_reads.tsv"
+  dedup_ratio:
+    fn: "dedup_ratio.tsv"
+  dedup_mean_umis:
+    fn: "dedup_mean_umis.tsv"
+  crosslinks_counts:
+    fn: "xlinks_counts.tsv"
+  crosslinks_ratio:
+    fn: "xlinks_ratio.tsv"
+  peaks:
+    fn: "total_peaks.tsv"
+  xlinks_in_peaks:
+    fn: "xlinks_in_peaks.tsv"
+  xlinksites_in_peaks:
+    fn: "xlinksites_in_peaks.tsv"
+  peaks_xlinksite_coverage:
+    fn: "peaks_xlinksite_coverage.tsv"
+  summary_type:
+    fn: "summary_type_metrics.tsv"
+  summary_subtype:
+    fn: '"summary_subtype_metrics.tsv'
+
 custom_content:
   order:
+    - clipqc
     - software-versions-by-process
     - software-versions-unique
-
 # Customise the module search patterns to speed up execution time
-sp:
-  samtools/stats:
-    fn: "*.stats"
-  samtools/flagstat:
-    fn: "*.flagstat"
-  samtools/idxstats:
-    fn: "*.idxstats*"
+# sp:
+#   samtools/stats:
+#     fn: "*.stats"
+#   samtools/flagstat:
+#     fn: "*.flagstat"
+#   samtools/idxstats:
+#     fn: "*.idxstats*"
+#   clipqc:
+#     fn: "*.txt"

conf/modules.config

@@ -34,6 +34,7 @@ process {
 ========================================================================================
 */
 
+
 if(params.run_genome_prep) {
     process {
         withName: '.*PREPARE_GENOME:GUNZIP_.*' {
@@ -71,7 +72,6 @@ if(params.run_genome_prep) {
                 path: { "${params.outdir}/00_genome" },
                 mode: "${params.publish_dir_mode}",
                 saveAs: { filename -> filename.equals('versions.yml') ? null : filename },
-                enabled: params.save_reference
             ]
         }
 
@@ -185,6 +185,7 @@ if(params.run_genome_prep) {
     }
 }
 
+
 /*
 ========================================================================================
     PRE-PROCESSING
@@ -1004,13 +1005,13 @@ if(params.run_reporting) {
             ]
         }
 
-        // withName: 'CLIPSEQ:CLIPSEQ_CLIPQC' {
-        //     publishDir = [
-        //         path: { "${params.outdir}/06_reports/clipqc" },
-        //         mode: "${params.publish_dir_mode}",
-        //         saveAs: { filename -> filename.equals('versions.yml') ? null : filename }
-        //     ]
-        // }
+         withName: 'NFCORE_CLIPSEQ:CLIPSEQ:CLIPQC' {
+             publishDir = [
+                 path: { "${params.outdir}/06_reports/clipqc" },
+                 mode: "${params.publish_dir_mode}",
+                 saveAs: { filename -> filename.equals('versions.yml') ? null : filename }
+             ]
+         }
 
         withName: 'NFCORE_CLIPSEQ:CLIPSEQ:MULTIQC' {
             ext.args   = params.multiqc_title ? "-v --title \"$params.multiqc_title\"" : '-v'

conf/test.config

@@ -20,9 +20,13 @@ params {
     max_time   = '6.h'
 
     // Input data
-    input  = 'https://raw.githubusercontent.com/nf-core/clipseq/refs/heads/feat-2-0/tests/test_new_samplesheet_FASTQ.csv'
+    input  = "./tests/test_new_samplesheet_FASTQ.csv"
     source = "fastq"
 
+    // Genome references
+    //fasta         = 's3://nf-core-awsmegatests/clipseq/input_data/reference/GRCh38.primary_assembly.genome.fa.gz'
+    //gtf           = 's3://nf-core-awsmegatests/clipseq/input_data/reference/gencode.v37.primary_assembly.annotation.gtf.gz'
+
     // Genome references
     fasta = "https://raw.githubusercontent.com/nf-core/test-datasets/clipseq/v_2_0/genome/yeast_MitoV.fa.gz"
     ncrna_fasta = "https://raw.githubusercontent.com/nf-core/test-datasets/clipseq/v_2_0/genome/homosapiens_smallRNA.fa.gz"
@@ -56,4 +60,8 @@ params {
 
     // Pipeline params
     umitools_bc_pattern = 'NNNNNNNNN'
+
+    // Don't call consensus
+    //consensus_peak        = false
+
 }

modules/local/clipqc/main.nf

@@ -12,6 +12,10 @@ process CLIPQC {
     path("icount/*")
     path("paraclu/*")
     path("clippy/*")
+    path("pureclip/*")
+    path("summary_type/*")
+    path("summary_subtype/*")
+    path("summary_gene/*")
 
     output:
     path "*.tsv"         , emit: tsv