For RNA-seq data, from bams to AS Information files.
Nie możesz wybrać więcej, niż 25 tematów Tematy muszą się zaczynać od litery lub cyfry, mogą zawierać myślniki ('-') i mogą mieć do 35 znaków.

samtools.wdl 653B

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