For RNA-seq data, from bams to AS Information files.
No puede seleccionar más de 25 temas Los temas deben comenzar con una letra o número, pueden incluir guiones ('-') y pueden tener hasta 35 caracteres de largo.

samtools.wdl 647B

1234567891011121314151617181920212223242526272829303132333435
  1. task samtools {
  2. File bam_file
  3. String bam_id=basename(bam_file)
  4. String bam_index = bam_id + ".bai"
  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}
  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_index
  24. }
  25. }