Add reference fasta as input

This commit is contained in:
Anders Goncalves da Silva 2021-02-12 11:26:55 -08:00
parent 6b22bb9782
commit a4bba006bd
No known key found for this signature in database
GPG key ID: 9DAECD8B113C8AFE

View file

@ -20,6 +20,7 @@ process IVAR_CONSENSUS {
input: input:
tuple val(meta), path(bam) tuple val(meta), path(bam)
path(fasta)
output: output:
tuple val(meta), path("*.fa") , emit: fasta tuple val(meta), path("*.fa") , emit: fasta
@ -33,6 +34,7 @@ process IVAR_CONSENSUS {
def save_mpileup = params.save_mpileup ? "tee ${prefix}.mpileup |" : "" def save_mpileup = params.save_mpileup ? "tee ${prefix}.mpileup |" : ""
""" """
samtools mpileup \\ samtools mpileup \\
--fasta-ref $fasta
-aa -A -d 0 -Q 0 \\ -aa -A -d 0 -Q 0 \\
$options.args2 \\ $options.args2 \\
$bam | \\ $bam | \\