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.

42 lines
1.1KB

  1. task rmats {
  2. String sample_id
  3. File bam_file
  4. File reference_gtf_file
  5. String docker
  6. String cluster
  7. String disk_size
  8. command <<<
  9. set -o pipefail
  10. set -e
  11. mkdir -p ${sample_id}/spladder_out
  12. spladder build --bams ${bam_file} \
  13. --annotation ${reference_gtf_file} \
  14. --outdir ${sample_id}/spladder_out
  15. find . -depth > fileList.txt
  16. >>>
  17. runtime {
  18. docker: docker
  19. cluster: cluster
  20. systemDisk: "cloud_ssd 40"
  21. dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  22. }
  23. output {
  24. File fileList = "fileList.txt"
  25. Array[File] AS_gff = glob("${sample_id}/spladder_out/*.gff3")
  26. Array[File] AS_pickle = glob("${sample_id}/spladder_out/*.pickle")
  27. Array[File] AS_gz = glob("${sample_id}/spladder_out/*.gz")
  28. Array[File] AS_hdf5 = glob("${sample_id}/spladder_out/*.hdf5")
  29. Array[File] spladder = glob("${sample_id}/spladder_out/spladder/*")
  30. Array[File] tmp = glob("${sample_id}/spladder_out/tmp/*")
  31. }
  32. }