@@ -251,136 +251,3 @@ task HaplotypeCaller_GATK4_VCF {
251251 docker : select_first ([runtime_attr .docker , default_attr .docker ])
252252 }
253253}
254-
255-
256- # This task is here because merging bamout files using Picard produces an error.
257- task MergeBamouts {
258-
259- input {
260- Array [File ] bams
261- String prefix
262-
263- RuntimeAttr ? runtime_attr_override
264- }
265-
266- Int disk_size = ceil (size (bams , "GiB" ) * 2 ) + 10
267-
268- command <<<
269-
270- set -euxo pipefail
271-
272- # Make sure we use all our processors:
273- np = $(cat /proc/cpuinfo | grep ^processor | tail -n1 | awk '{print $NF+1}' )
274-
275- ithreads = ${np }
276-
277- # If the number of processors = 1, then `let` will return 1 here:
278- # So we need to turn off `set -e` for this command:
279- set +e
280- mthreads = $((np - 1 ))
281- set -e
282-
283- samtools merge -@${mthreads } ~{prefix }.bam ~{sep =" " bams }
284- samtools index -@${ithreads } ~{prefix }.bam
285- mv ~{prefix }.bam.bai ~{prefix }.bai
286- >>>
287-
288- #########################
289- RuntimeAttr default_attr = object {
290- cpu_cores : 1 ,
291- mem_gb : 4 ,
292- disk_gb : disk_size ,
293- boot_disk_gb : 25 ,
294- preemptible_tries : 1 ,
295- max_retries : 1 ,
296- docker : "us.gcr.io/broad-dsp-lrma/lr-basic:0.1.1"
297- }
298-
299- RuntimeAttr runtime_attr = select_first ([runtime_attr_override , default_attr ])
300- runtime {
301- cpu : select_first ([runtime_attr .cpu_cores , default_attr .cpu_cores ])
302- memory : select_first ([runtime_attr .mem_gb , default_attr .mem_gb ]) + " GiB"
303- disks : "local-disk " + select_first ([runtime_attr .disk_gb , default_attr .disk_gb ]) + " HDD"
304- bootDiskSizeGb : select_first ([runtime_attr .boot_disk_gb , default_attr .boot_disk_gb ])
305- preemptible : select_first ([runtime_attr .preemptible_tries , default_attr .preemptible_tries ])
306- maxRetries : select_first ([runtime_attr .max_retries , default_attr .max_retries ])
307- docker : select_first ([runtime_attr .docker , default_attr .docker ])
308- }
309-
310- output {
311- File output_bam = "~{prefix }.bam"
312- File output_bam_index = "~{prefix }.bai"
313- }
314- }
315-
316- task ReblockGVCF {
317-
318- input {
319- File gvcf
320- File gvcf_index
321-
322- File ref_fasta
323- File ref_fasta_fai
324- File ref_dict
325-
326- String prefix
327-
328- Array [Int ] gq_blocks = [20 , 30 , 40 ]
329-
330- Float ? tree_score_cutoff
331-
332- Array [String ]? annotations_to_keep
333-
334- RuntimeAttr ? runtime_attr_override
335- }
336-
337- Int disk_size = ceil ((size (gvcf , "GiB" ) * 4 ) + size (ref_fasta , "GiB" ) + size (ref_fasta_fai , "GiB" ) + size (ref_dict , "GiB" ) + 10 )
338-
339- String annotations_to_keep_arg = if defined (annotations_to_keep ) then "--annotations-to-keep" else ""
340-
341- command {
342- set -euxo pipefail
343-
344- gatk --java-options "-Xms3000m -Xmx3000m" \
345- ReblockGVCF \
346- -R ~{ref_fasta } \
347- -V ~{gvcf } \
348- -do -qual-approx \
349- -G StandardAnnotation -G StandardHCAnnotation \
350- -A AssemblyComplexity \
351- --annotate-with-num-discovered-alleles \
352- --floor-blocks \
353- -GQB ~{sep =" -GQB " gq_blocks } \
354- ~{"--tree-score-threshold-to-no-call " + tree_score_cutoff } \
355- ~{annotations_to_keep_arg } ~{sep =" --annotations-to-keep " annotations_to_keep } \
356- -O ~{prefix }.rb.g.vcf.gz
357- }
358-
359- #########################
360- RuntimeAttr default_attr = object {
361- cpu_cores : 2 ,
362- mem_gb : 4 ,
363- disk_gb : disk_size ,
364- boot_disk_gb : 25 ,
365- preemptible_tries : 1 ,
366- max_retries : 1 ,
367- docker : "broadinstitute/gatk-nightly:2024-04-16-4.5.0.0-25-g986cb1549-NIGHTLY-SNAPSHOT"
368- }
369- # TODO: Fix this docker image to a stable version after the next GATK release!
370-
371- RuntimeAttr runtime_attr = select_first ([runtime_attr_override , default_attr ])
372- runtime {
373- cpu : select_first ([runtime_attr .cpu_cores , default_attr .cpu_cores ])
374- memory : select_first ([runtime_attr .mem_gb , default_attr .mem_gb ]) + " GiB"
375- disks : "local-disk " + select_first ([runtime_attr .disk_gb , default_attr .disk_gb ]) + " SSD"
376- bootDiskSizeGb : select_first ([runtime_attr .boot_disk_gb , default_attr .boot_disk_gb ])
377- preemptible : select_first ([runtime_attr .preemptible_tries , default_attr .preemptible_tries ])
378- maxRetries : select_first ([runtime_attr .max_retries , default_attr .max_retries ])
379- docker : select_first ([runtime_attr .docker , default_attr .docker ])
380- }
381-
382- output {
383- File output_gvcf = "~{prefix }.rb.g.vcf.gz"
384- File output_gvcf_index = "~{prefix }.rb.g.vcf.gz.tbi"
385- }
386- }
0 commit comments