index concat files

nf-core · Apr 20, 2024 · 0ed8708 · 0ed8708
1 parent b9b4873
commit 0ed8708
Show file tree

Hide file tree

Showing 2 changed files with 21 additions and 1 deletion.
diff --git a/conf/quilt_subworkflow.config b/conf/quilt_subworkflow.config
@@ -119,4 +119,17 @@ process {
         ]
     }
 
+        withName: 'NFCORE_PHASEIMPUTE:PHASEIMPUTE:VCF_CONCATENATE_BCFTOOLS:BCFTOOLS_INDEX' {
+        ext.args     = {[
+                        "--tbi",
+                        ].join(" ").trim()}
+
+        publishDir = [
+            [
+                path: { "${params.outdir}/quilt_impute/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}/concat" },
+                mode: params.publish_dir_mode,
+            ],
+        ]
+    }
+
 }
diff --git a/subworkflows/local/vcf_concatenate_bcftools/vcf_concatenate_bcftools.nf b/subworkflows/local/vcf_concatenate_bcftools/vcf_concatenate_bcftools.nf
@@ -1,4 +1,5 @@
 include { BCFTOOLS_CONCAT } from '../../../modules/nf-core/bcftools/concat/main'
+include { BCFTOOLS_INDEX  } from '../../../modules/nf-core/bcftools/index/main'
 
 workflow VCF_CONCATENATE_BCFTOOLS {
 
@@ -17,7 +18,13 @@ workflow VCF_CONCATENATE_BCFTOOLS {
     // Ligate and concatenate chunks
     BCFTOOLS_CONCAT(ch_vcf_tbi_grouped)
 
+    // Index concatenated VCF
+    BCFTOOLS_INDEX(BCFTOOLS_CONCAT.out.vcf)
+
+    // Join VCFs and TBIs
+    ch_imputed_vcf_tbi = BCFTOOLS_CONCAT.out.vcf.join(BCFTOOLS_INDEX.out.tbi)
+
     emit:
-    ch_concat_vcf =  BCFTOOLS_CONCAT.out.vcf
+    ch_imputed_vcf_tbi                          // channel:  [ meta, vcf, tbi ]
 
     }