Skip to content

Commit

Permalink
fixed disabling of artic normalisation
Browse files Browse the repository at this point in the history
If the disabling was activated in the pipeline (`--artic_normalize False`), no parameter was handed over to the artic pipeline, so that artic used the default value of 200.
Normalisation is now correctly turned off by passing `--normalise 0` to the artic pipeline.
  • Loading branch information
MarieLataretu committed Sep 13, 2023
1 parent aa993eb commit 7e1bf51
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions workflows/process/artic.nf
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@ process artic_medaka {
tuple val(name), path("${name}.fail.vcf"), emit: vcf_fail

script:
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : ''
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : '--normalise 0'
"""
artic minion --medaka \
--medaka-model ${params.medaka_model} \
Expand Down Expand Up @@ -89,7 +89,7 @@ process artic_medaka_custom_bed {
tuple val(name), path("${name}.coverage_mask.txt"), emit: coverage_mask
tuple val(name), path("${name}.fail.vcf"), emit: vcf_fail
script:
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : ''
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : '--normalise 0'
"""
# create a new primer dir as input for artic
mkdir -p primer_scheme/nCoV-2019
Expand Down Expand Up @@ -167,7 +167,7 @@ process artic_nanopolish {
tuple val(name), path("${name}.coverage_mask.txt"), emit: coverage_mask
tuple val(name), path("${name}.fail.vcf"), emit: vcf_fail
script:
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : ''
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : '--normalise 0'
"""
artic minion --minimap2 \
${normalise_arg} \
Expand Down Expand Up @@ -234,7 +234,7 @@ process artic_nanopolish_custom_bed {
tuple val(name), path("${name}.coverage_mask.txt"), emit: coverage_mask
tuple val(name), path("${name}.fail.vcf"), emit: vcf_fail
script:
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : ''
def normalise_arg = normalise_threshold ? "--normalise ${normalise_threshold}" : '--normalise 0'
"""
# create a new primer dir as input for artic
mkdir -p primer_scheme/nCoV-2019
Expand Down

0 comments on commit 7e1bf51

Please sign in to comment.