Add module bcftools/norm to sarek (in progress)

conf/modules/post_variant_calling.config

@@ -16,6 +16,18 @@
 
 process {
 
+    withName: 'GERMLINE_VCFS_NORM'{
+        ext.args = { [
+            '--multiallelics - both', //split multiallelic sites into biallelic records and both SNPs and indels should be merged separately into two records
+            '--rm-dup all'            //output only the first instance of a record which is present multiple times
+        ].join(' ') }
+        ext.when   = { params.concatenate_vcfs }
+        publishDir = [
+            mode: params.publish_dir_mode,
+                path: { "${params.outdir}/variant_calling/concat/${meta.id}/" }
+        ]
+    }
+
     withName: 'GERMLINE_VCFS_CONCAT'{
         ext.args   = { "-a" }
         ext.when   = { params.concatenate_vcfs }

subworkflows/local/post_variantcalling/main.nf

@@ -8,13 +8,14 @@ workflow POST_VARIANTCALLING {
 
     take:
     vcfs
+    fasta
     concatenate_vcfs
 
     main:
     versions = Channel.empty()
 
     if (concatenate_vcfs){
-        CONCATENATE_GERMLINE_VCFS(vcfs)
+        CONCATENATE_GERMLINE_VCFS(vcfs, fasta)
 
         vcfs = vcfs.mix(CONCATENATE_GERMLINE_VCFS.out.vcfs)
         versions = versions.mix(CONCATENATE_GERMLINE_VCFS.out.versions)

subworkflows/local/vcf_concatenate_germline/main.nf

@@ -1,42 +1,59 @@
-//
 // CONCATENATE Germline VCFs
-//
 
 // Concatenation of germline vcf-files
-include { ADD_INFO_TO_VCF                                     } from '../../../modules/local/add_info_to_vcf/main'
-include { TABIX_BGZIPTABIX as TABIX_EXT_VCF                   } from '../../../modules/nf-core/tabix/bgziptabix/main'
-include { BCFTOOLS_CONCAT  as GERMLINE_VCFS_CONCAT            } from '../../../modules/nf-core/bcftools/concat/main'
-include { BCFTOOLS_SORT    as GERMLINE_VCFS_CONCAT_SORT       } from '../../../modules/nf-core/bcftools/sort/main'
-include { TABIX_TABIX      as TABIX_GERMLINE_VCFS_CONCAT_SORT } from '../../../modules/nf-core/tabix/tabix/main'
+include { ADD_INFO_TO_VCF  } from '../../../modules/local/add_info_to_vcf/main'
+include { TABIX_BGZIPTABIX as TABIX_EXT_VCF } from '../../../modules/nf-core/tabix/bgziptabix/main'
+include { BCFTOOLS_NORM as GERMLINE_VCFS_NORM } from '../../../modules/nf-core/bcftools/norm/main'
+include { BCFTOOLS_CONCAT as GERMLINE_VCFS_CONCAT } from '../../../modules/nf-core/bcftools/concat/main'
+include { BCFTOOLS_SORT as GERMLINE_VCFS_CONCAT_SORT } from '../../../modules/nf-core/bcftools/sort/main'
+include { TABIX_TABIX as TABIX_GERMLINE_VCFS_CONCAT_SORT } from '../../../modules/nf-core/tabix/tabix/main'
 
 workflow CONCATENATE_GERMLINE_VCFS {
 
     take:
     vcfs
+    fasta
 
     main:
     versions = Channel.empty()
 
-    // Concatenate vcf-files
+    // Add additional information to VCF files
     ADD_INFO_TO_VCF(vcfs)
+
+    // Compress the VCF files with bgzip
     TABIX_EXT_VCF(ADD_INFO_TO_VCF.out.vcf)
 
+    // Normalize the VCF files with BCFTOOLS_NORM
+    GERMLINE_VCFS_NORM(vcf: ADD_INFO_TO_VCF.out.vcf, fasta: fasta)
+
+    // Compress the normalized VCF files with bgzip
+    TABIX_EXT_VCF(GERMLINE_VCFS_NORM.out.vcf)
+
+    // Index the compressed normalized VCF files
+    TABIX_GERMLINE_VCFS_CONCAT_SORT(TABIX_EXT_VCF.out.gz)
+
     // Gather vcfs and vcf-tbis for concatenating germline-vcfs
-    germline_vcfs_with_tbis = TABIX_EXT_VCF.out.gz_tbi.map{ meta, vcf, tbi -> [ meta.subMap('id'), vcf, tbi ] }.groupTuple()
+    germline_vcfs_with_tbis = TABIX_GERMLINE_VCFS_CONCAT_SORT.out.map { meta, vcf, tbi -> [meta.subMap('id'), vcf, tbi] }.groupTuple()
 
+    // Concatenate the VCF files
     GERMLINE_VCFS_CONCAT(germline_vcfs_with_tbis)
+
+    // Sort the concatenated VCF files
     GERMLINE_VCFS_CONCAT_SORT(GERMLINE_VCFS_CONCAT.out.vcf)
+
+    // Index the sorted concatenated VCF files
     TABIX_GERMLINE_VCFS_CONCAT_SORT(GERMLINE_VCFS_CONCAT_SORT.out.vcf)
 
     // Gather versions of all tools used
     versions = versions.mix(ADD_INFO_TO_VCF.out.versions)
     versions = versions.mix(TABIX_EXT_VCF.out.versions)
+    versions = versions.mix(GERMLINE_VCFS_NORM.out.versions)
     versions = versions.mix(GERMLINE_VCFS_CONCAT.out.versions)
-    versions = versions.mix(GERMLINE_VCFS_CONCAT.out.versions)
+    versions = versions.mix(GERMLINE_VCFS_CONCAT_SORT.out.versions)
     versions = versions.mix(TABIX_GERMLINE_VCFS_CONCAT_SORT.out.versions)
 
     emit:
-    vcfs = germline_vcfs_with_tbis // post processed vcfs
-
+    vcfs = TABIX_GERMLINE_VCFS_CONCAT_SORT.out.gz_tbi // post-processed VCFs
     versions // channel: [ versions.yml ]
 }
+

workflows/sarek/main.nf

@@ -794,7 +794,8 @@ workflow SAREK {
 
         // POST VARIANTCALLING
         POST_VARIANTCALLING(BAM_VARIANT_CALLING_GERMLINE_ALL.out.vcf_all,
-                            params.concatenate_vcfs)
+                fasta,
+                params.concatenate_vcfs)
 
         // Gather vcf files for annotation and QC
         vcf_to_annotate = Channel.empty()