@@ -224,15 +224,29 @@ if (!params.skip_realignment) {
224
224
}
225
225
}
226
226
else {
227
- realigned_bams = deduplicated_bams
227
+ process skipRealignmentAroundindels {
228
+ cpus 1
229
+ memory ' 16m'
230
+ tag " ${ name} "
231
+
232
+ input:
233
+ set name, bam_name, type, file(bam), file(bai) from deduplicated_bams
234
+
235
+ output:
236
+ set val(name), val(bam_name), val(type), file(" ${ bam} " ), file(" ${ bai} " ) into realigned_bams
237
+
238
+ """
239
+ echo "ZZZZZ..."
240
+ """
241
+ }
228
242
}
229
243
230
244
if (! params. skip_bqsr) {
231
245
process baseQualityScoreRecalibration {
232
246
cpus 3
233
247
memory ' 4g'
234
248
module ' java/1.8.0' // GATK requires Java 8
235
- publishDir " ${ publish_dir} " , mode: " move "
249
+ publishDir " ${ publish_dir} " , mode: " copy "
236
250
tag " ${ name} "
237
251
238
252
input:
@@ -266,21 +280,21 @@ else {
266
280
process createOutput {
267
281
cpus 1
268
282
memory ' 1g'
269
- publishDir " ${ publish_dir} " , mode: " move "
283
+ publishDir " ${ publish_dir} " , mode: " copy "
270
284
tag " ${ name} "
271
285
272
286
input:
273
287
set name, bam_name, type, file(bam), file(bai) from realigned_bams
274
288
275
289
output:
276
290
set val(" ${ name} " ), val(" ${ type} " ), val(" ${ publish_dir} /${ bam_name} .preprocessed.bam" ) into recalibrated_bams
277
- file " ${ bam_name} .preprocessed.bam" into recalibrated_bam
278
- file " ${ bam_name} .preprocessed.bai" into recalibrated_bai
291
+ file " ${ bam_name} .preprocessed.bam" into recalibrated_bam
292
+ file " ${ bam_name} .preprocessed.bai" into recalibrated_bai
279
293
280
- """
281
- mv ${ bam} ${ bam_name} .preprocessed.bam
282
- mv ${ bai} ${ bam_name} .preprocessed.bai
283
- """
294
+ """
295
+ cp ${ bam} ${ bam_name} .preprocessed.bam
296
+ cp ${ bai} ${ bam_name} .preprocessed.bai
297
+ """
284
298
}
285
299
}
286
300
0 commit comments