Skip to content

Commit d488dde

Browse files
committed
Clearer intersect step handling
1 parent 4aae304 commit d488dde

File tree

1 file changed

+8
-7
lines changed

1 file changed

+8
-7
lines changed

main.nf

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -203,17 +203,14 @@ process mask_beds {
203203

204204

205205
// To make intersect files we need to combine them into one channel with
206-
// toSortedList() (fermi is before manta in alphabet). And also figure out if we
207-
// have one or two files, therefore the tap and count_vcfs.
208-
masked_vcfs.tap { count_vcfs_tmp }
209-
.tap { masked_vcfs }
206+
// toSortedList() (fermi is before manta in alphabet).
207+
masked_vcfs.tap { masked_vcfs }
208+
.filter( ~/manta|fermikit/ )
210209
.toSortedList().set { intersect_input }
211-
count_vcfs_tmp.count().set { count_vcfs }
212210

213211
process intersect_files {
214212
input:
215213
set file(fermi_vcf), file(manta_vcf) from intersect_input
216-
val nvcfs from count_vcfs
217214
output:
218215
file "combined_masked.vcf" into intersections
219216

@@ -224,7 +221,7 @@ process intersect_files {
224221
module 'bioinfo-tools'
225222
module "$params.modules.bedtools"
226223

227-
when: nvcfs == 2
224+
when: 'make_intersect' in workflowSteps
228225

229226
script:
230227
"""
@@ -458,5 +455,9 @@ def processWorkflowSteps(steps) {
458455
workflowSteps.push( 'fastq' )
459456
}
460457

458+
if ('manta' in workflowSteps && 'fermikit' in workflowSteps) {
459+
workflowSteps.push( 'make_intersect' )
460+
}
461+
461462
return workflowSteps
462463
}

0 commit comments

Comments
 (0)