For RNA-seq data, from bams to AS Information files.
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.

31 line
528B

  1. task samtools {
  2. File bam_file
  3. String bam_id=basename(bam_file)
  4. String docker
  5. String cluster
  6. String disk_size
  7. command <<<
  8. set -o pipefail
  9. set -e
  10. cp ${bam_file} ${bam_id}
  11. samtools index ${bam_file}
  12. >>>
  13. runtime {
  14. docker: docker
  15. cluster: cluster
  16. systemDisk: "cloud_ssd 40"
  17. dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  18. }
  19. output {
  20. File bam = "$(bam_id)"
  21. File bam_index = "$(bam_id).bai"
  22. }
  23. }