add subset ped step at start of CleanVcf (#250)

broadinstitute · Dec 21, 2021 · bf3704b · bf3704b
1 parent 629d382
commit bf3704b
Show file tree

Hide file tree

Showing 2 changed files with 22 additions and 2 deletions.
diff --git a/wdl/CleanVcf.wdl b/wdl/CleanVcf.wdl
@@ -2,6 +2,7 @@ version 1.0
 
 import "CleanVcfChromosome.wdl" as CleanVcfChromosome
 import "TasksMakeCohortVcf.wdl" as MiniTasks
+import "Utils.wdl" as util
 
 workflow CleanVcf {
   input {
@@ -43,6 +44,23 @@ workflow CleanVcf {
     RuntimeAttr? runtime_override_combine_clean_vcf_2
     RuntimeAttr? runtime_override_combine_revised_4
     RuntimeAttr? runtime_override_combine_multi_ids_4
+    RuntimeAttr? runtime_attr_ids_from_vcf
+    RuntimeAttr? runtime_attr_subset_ped
+  }
+
+  call util.GetSampleIdsFromVcf {
+    input:
+      vcf = complex_genotype_vcfs[0],
+      sv_base_mini_docker = sv_base_mini_docker,
+      runtime_attr_override = runtime_attr_ids_from_vcf
+  }
+  call util.SubsetPedFile {
+    input:
+      ped_file = merged_ped_file,
+      sample_list = GetSampleIdsFromVcf.out_file,
+      subset_name = "vcf_samples",
+      sv_base_mini_docker = sv_base_mini_docker,
+      runtime_attr_override = runtime_attr_subset_ped
   }
 
   #Scatter per chromosome
@@ -55,7 +73,7 @@ workflow CleanVcf {
         vcf=complex_genotype_vcfs[i],
         contig=contig,
         background_list=complex_resolve_background_fail_lists[i],
-        ped_file=merged_ped_file,
+        ped_file=SubsetPedFile.ped_subset_file,
         bothsides_pass_list=complex_resolve_bothside_pass_lists[i],
         allosome_fai=allosome_fai,
         prefix=cohort_name,

diff --git a/wdl/MakeCohortVcf.wdl b/wdl/MakeCohortVcf.wdl
@@ -328,7 +328,9 @@ workflow MakeCohortVcf {
       runtime_override_split_include_list=runtime_override_split_include_list,
       runtime_override_combine_clean_vcf_2=runtime_override_combine_clean_vcf_2,
       runtime_override_combine_revised_4=runtime_override_combine_revised_4,
-      runtime_override_combine_multi_ids_4=runtime_override_combine_multi_ids_4
+      runtime_override_combine_multi_ids_4=runtime_override_combine_multi_ids_4,
+      runtime_attr_ids_from_vcf=runtime_attr_ids_from_vcf,
+      runtime_attr_subset_ped=runtime_attr_subset_ped
   }
 
   Array[String] contigs = transpose(read_tsv(contig_list))[0]