For RNA-seq data, from bams to AS Information files.
Nelze vybrat více než 25 témat Téma musí začínat písmenem nebo číslem, může obsahovat pomlčky („-“) a může být dlouhé až 35 znaků.

před 4 roky
před 4 roky
před 4 roky
před 4 roky
123456789101112131415161718192021222324252627282930313233343536373839404142434445464748
  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. }