diff --git a/software/ivar/consensus/main.nf b/software/ivar/consensus/main.nf index a0288365..b30e8c8b 100644 --- a/software/ivar/consensus/main.nf +++ b/software/ivar/consensus/main.nf @@ -20,6 +20,7 @@ process IVAR_CONSENSUS { input: tuple val(meta), path(bam) + path(fasta) output: tuple val(meta), path("*.fa") , emit: fasta @@ -33,6 +34,7 @@ process IVAR_CONSENSUS { def save_mpileup = params.save_mpileup ? "tee ${prefix}.mpileup |" : "" """ samtools mpileup \\ + --fasta-ref $fasta -aa -A -d 0 -Q 0 \\ $options.args2 \\ $bam | \\