Merge pull request #356 from MaxUlysse/dev_controlfreec

add option for controlfreec
nf-core · Mar 9, 2021 · f473afe · f473afe
2 parents f365b81 + eb0f8b4
commit f473afe
Show file tree

Hide file tree

Showing 3 changed files with 19 additions and 7 deletions.
diff --git a/main.nf b/main.nf
@@ -93,6 +93,8 @@ def helpMessage() {
                                           Default: ${params.cf_ploidy}
       --cf_window                   [int] Control-FREEC window size
                                           Default: Disabled
+      --cf_contamination_adjustment[bool] Control-FREEC contaminationAdjustment
+                                          Default: Disabled
       --generate_gvcf              [bool] Enable g.vcf output from GATK HaplotypeCaller
       --no_strelka_bp              [bool] Will not use Manta candidateSmallIndels for Strelka (not recommended by Best Practices)
       --pon                        [file] Panel-of-normals VCF (bgzipped) for GATK Mutect2 / Sentieon TNscope
@@ -427,9 +429,10 @@ if ('ascat' in tools) {
 
 if ('controlfreec' in tools) {
     summary['Control-FREEC'] = "Options"
-    if (params.cf_window)    summary['window']                 = params.cf_window
-    if (params.cf_coeff)     summary['coefficientOfVariation'] = params.cf_coeff
-    if (params.cf_ploidy)    summary['ploidy']                 = params.cf_ploidy
+    if (params.cf_window)                   summary['window']                  = params.cf_window
+    if (params.cf_coeff)                    summary['coefficientOfVariation']  = params.cf_coeff
+    if (params.cf_ploidy)                   summary['ploidy']                  = params.cf_ploidy
+    if (params.cf_contamination_adjustment) summary['contaminationAdjustment'] = params.cf_contamination_adjustment
 }
 
 if ('haplotypecaller' in tools)             summary['GVCF']       = params.generate_gvcf ? 'Yes' : 'No'
@@ -3274,7 +3277,7 @@ process ControlFREEC {
     breakPointThreshold = params.target_bed ? "1.2" : "0.8"
     breakPointType = params.target_bed ? "4" : "2"
     mappabilitystr = params.mappability ? "gemMappabilityFile = \${PWD}/${mappability}" : ""
-
+    contamination_adjustment = params.cf_contamination_adjustment ? "contaminationAdjustment = TRUE" : ""
     """
     touch ${config}
     echo "[general]" >> ${config}
@@ -3292,6 +3295,7 @@ process ControlFREEC {
     echo "${window}" >> ${config}
     echo "${coeffvar}" >> ${config}
     echo "${mappabilitystr}" >> ${config}
+    echo "${contamination_adjustment}" >> ${config}
     echo "" >> ${config}
     
     echo "[control]" >> ${config}
@@ -3358,6 +3362,7 @@ process ControlFREECSingle {
     breakPointThreshold = params.target_bed ? "1.2" : "0.8"
     breakPointType = params.target_bed ? "4" : "2"
     mappabilitystr = params.mappability ? "gemMappabilityFile = \${PWD}/${mappability}" : ""
+    contamination_adjustment = params.cf_contamination_adjustment ? "contaminationAdjustment = TRUE" : ""
 
     """
     touch ${config}
@@ -3376,6 +3381,7 @@ process ControlFREECSingle {
     echo "${window}" >> ${config}
     echo "${coeffvar}" >> ${config}
     echo "${mappabilitystr}" >> ${config}
+    echo "${contamination_adjustment}" >> ${config}
     echo "" >> ${config}
 
     echo "[sample]" >> ${config}

diff --git a/nextflow.config b/nextflow.config
@@ -49,9 +49,10 @@ params {
   // Variant Calling
   ascat_ploidy = null // Use default value
   ascat_purity = null // Use default value
-  cf_coeff = "0.05"  // default value for Control-FREEC
-  cf_ploidy = "2"     // you can use 2,3,4
-  cf_window = null      // by default we are not using this in Control-FREEC
+  cf_coeff = "0.05"                  // default value for Control-FREEC
+  cf_ploidy = "2"                    // you can use 2,3,4
+  cf_window = null                   // by default we are not using this in Control-FREEC
+  cf_contamination_adjustment = null // by default we are not using this in Control-FREEC
   generate_gvcf = null // g.vcf are not produced by HaplotypeCaller by default
   no_strelka_bp = null // Strelka will use Manta candidateSmallIndels if available
   pon = false // No default PON (Panel of Normals) file for GATK Mutect2 / Sentieon TNscope

diff --git a/nextflow_schema.json b/nextflow_schema.json
@@ -263,6 +263,11 @@
                     "fa_icon": "fas fa-wrench",
                     "description": "Overwrite Control-FREEC coefficientOfVariation"
                 },
+                "cf_contamination_adjustment": {
+                    "type": "boolean",
+                    "fa_icon": "fas fa-wrench",
+                    "description": "Overwrite Control-FREEC contaminationAdjustement"
+                },
                 "cf_ploidy": {
                     "type": "string",
                     "default": "2",