Skip to content

Commit 2b5f38d

Browse files
committed
changes specific to UW DV call
1 parent 309f2ff commit 2b5f38d

File tree

2 files changed

+7
-5
lines changed

2 files changed

+7
-5
lines changed

wdl/pipelines/TechAgnostic/VariantCalling/CallSmallVariants.wdl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -81,7 +81,7 @@ workflow Work {
8181
String sex
8282
String prefix
8383

84-
Array[Pair[String, Pair[File, File]]] per_chr_bam_bai_and_id
84+
Array[Pair[String, Pair[File, File]]]? per_chr_bam_bai_and_id
8585

8686
Boolean is_ont
8787
Boolean is_r10_4_pore_or_later
@@ -213,7 +213,7 @@ workflow Work {
213213
214214
bam = bam,
215215
bai = bai,
216-
per_chr_bam_bai_and_id = per_chr_bam_bai_and_id,
216+
per_chr_bam_bai_and_id = select_first([per_chr_bam_bai_and_id]),
217217
is_ont = is_ont,
218218
219219
unphased_vcf = DV.vcf,

wdl/pipelines/TechAgnostic/VariantCalling/CallVariantsReadBased.wdl

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -96,8 +96,7 @@ workflow CallVariants {
9696
call GU.CollapseArrayOfStrings as get_zones {input: input_array = gcp_zones, joiner = " "}
9797
String wdl_parsable_zones = get_zones.collapsed
9898
99-
# needed for whatshap phasing anyway, so this can be used by SV calling
100-
call ShardWholeGenome.Split as SplitBamByChr { input: ref_dict = ref_bundle.dict, bam = bam, bai = bai, }
99+
101100
102101
######################################################################
103102
# Block for small variants handling
@@ -113,7 +112,7 @@ workflow CallVariants {
113112
sex = sex,
114113
prefix = prefix,
115114
116-
per_chr_bam_bai_and_id = SplitBamByChr.id_bam_bai_of_shards,
115+
# per_chr_bam_bai_and_id = SplitBamByChr.id_bam_bai_of_shards,
117116
118117
is_ont = is_ont,
119118
is_r10_4_pore_or_later = is_r10_4_pore_or_later,
@@ -144,6 +143,9 @@ workflow CallVariants {
144143
String sv_dir = sub(gcs_out_dir, "/$", "") + "/variants/sv"
145144
if (defined(sv_calling_options_json)) {
146145

146+
# needed for whatshap phasing anyway, so this can be used by SV calling
147+
call ShardWholeGenome.Split as SplitBamByChr { input: ref_dict = ref_bundle.dict, bam = bam, bai = bai, }
148+
147149
SVCallingConfig sv_options = read_json(select_first([sv_calling_options_json]))
148150
149151
call CallStructuralVariants.Work as SVjob {

0 commit comments

Comments
 (0)