bamdst is a lightweight tool to stat the depth coverage of target regions of bam file(s).
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

41 lines
1.3KB

  1. task bamdst {
  2. File bam
  3. File bai
  4. File bed
  5. String sample
  6. String docker
  7. String cluster_config
  8. String disk_size
  9. command <<<
  10. set -o pipefail
  11. set -e
  12. mkdir -p /cromwell_root/tmp/bamdst
  13. bamdst -p ${bed} -o /cromwell_root/tmp/bamdst ${bam}
  14. cd /cromwell_root/tmp/bamdst
  15. cp chromosomes.report ${sample}_chromosomes.report
  16. cp coverage.report ${sample}_coverage.report
  17. cp depth_distribution.plot ${sample}_depth_distribution.plot
  18. cp insertsize.plot ${sample}_insertsize.plot
  19. cp region.tsv.gz ${sample}_region.tsv.gz
  20. cp uncover.bed ${sample}_uncover.bed
  21. >>>
  22. runtime {
  23. docker: docker
  24. cluster: cluster_config
  25. systemDisk: "cloud_ssd 40"
  26. dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  27. }
  28. output {
  29. File chromosomes_report = "${sample}_chromosomes.report"
  30. File coverage_report = "${sample}_coverage.report"
  31. File depth_distribution_plot = "${sample}_depth_distribution.plot"
  32. File depth_tsv_gz = "${sample}_depth.tsv.gz"
  33. File insertsize_plot = "${sample}_insertsize.plot"
  34. File region_tsv_gz = "${sample}_region.tsv.gz"
  35. File uncover_bed = "${sample}_uncover.bed"
  36. }
  37. }