fix issue with the multiallelic filter + document the skip of the multiallelic filter

Author: priesgo

README.md

@@ -55,6 +55,8 @@ Optional input:
     * --skip_decompose_complex: flag indicating not to split complex variants (ie: MNVs and combinations of SNVs and indels)
     * --filter: specify the filter to apply if any (e.g.: PASS), only variants with this value will be kept
     * --input_bams: a tab-separated values file containing in each row the sample name, tumor and normal BAM files for annotation with Vafator
+    * --skip_multiallelic_filter: after VAFator annotations if any multiallelic variant is present (ie: two different 
+    mutations in the same position) only the highest VAF variant is kept unless this flag is passed
     * --snpeff_organism: the SnpEff organism name (eg: hg19, hg38, GRCh37.75, GRCh38.99)
     * --snpeff_datadir: the SnpEff data folder where the reference genomes were previously downloaded. Required if --snpeff_organism is provided
     * --snpeff_args: additional SnpEff arguments

modules/04_vafator.nf

@@ -45,7 +45,7 @@ process MULTIALLELIC_FILTER {
     tuple val(name), file(vcf)
 
     output:
-    tuple file("${vcf.baseName}.filtered_multiallelics.vcf"), emit: filtered_vcf
+    tuple val(name), file("${vcf.baseName}.filtered_multiallelics.vcf"), emit: filtered_vcf
 
     script:
     """

nextflow.config

@@ -27,7 +27,7 @@ env {
 // Capture exit codes from upstream processes when piping
 process.shell = ['/bin/bash', '-euo', 'pipefail']
 
-VERSION = '2.1.0'
+VERSION = '2.1.1'
 
 cleanup=true
 
@@ -63,6 +63,8 @@ Optional input:
     * --filter: specify a comma-separated list of filters to apply (e.g.: PASS,.), only variants with these values will be kept. If not provided all varianst are kept
     * --vcf-without-ad: indicate when the VCFs to normalize do not have the FORMAT/AD annotation
     * --input_bams: a tab-separated values file containing in each row the sample name, tumor and normal BAM files for annotation with Vafator
+    * --skip_multiallelic_filter: after VAFator annotations if any multiallelic variant is present (ie: two different
+    mutations in the same position) only the highest VAF variant is kept unless this flag is passed
     * --snpeff_organism: the SnpEff organism name (eg: hg19, hg38, GRCh37.75, GRCh38.99)
     * --snpeff_datadir: the SnpEff data folder where the reference genomes were previously downloaded. Required if --snpeff_organism is provided
     * --snpeff_args: additional SnpEff arguments