For RNA-seq data, from bams to AS Information files.
Vous ne pouvez pas sélectionner plus de 25 sujets Les noms de sujets doivent commencer par une lettre ou un nombre, peuvent contenir des tirets ('-') et peuvent comporter jusqu'à 35 caractères.

48 lignes
1.2KB

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