For RNA-seq data, from bams to AS Information files.
Du kannst nicht mehr als 25 Themen auswählen Themen müssen entweder mit einem Buchstaben oder einer Ziffer beginnen. Sie können Bindestriche („-“) enthalten und bis zu 35 Zeichen lang sein.

47 Zeilen
1.2KB

  1. task spladder {
  2. String sample_id
  3. File bam
  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} \
  13. --annotation ${reference_gtf_file} \
  14. --ignore-mismatches \
  15. --confidence 2 \
  16. --readlen 150 \
  17. --sparse bam y \
  18. --parallel 8 \
  19. --outdir ${sample_id}/spladder_out
  20. find . -depth > fileList.txt
  21. >>>
  22. runtime {
  23. docker: docker
  24. cluster: cluster
  25. systemDisk: "cloud_ssd 40"
  26. dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  27. }
  28. output {
  29. File fileList = "fileList.txt"
  30. Array[File] AS_gff = glob("${sample_id}/spladder_out/*.gff3")
  31. Array[File] AS_pickle = glob("${sample_id}/spladder_out/*.pickle")
  32. Array[File] AS_gz = glob("${sample_id}/spladder_out/*.gz")
  33. Array[File] AS_hdf5 = glob("${sample_id}/spladder_out/*.hdf5")
  34. Array[File] spladder = glob("${sample_id}/spladder_out/spladder/*")
  35. Array[File] tmp = glob("${sample_id}/spladder_out/tmp/*")
  36. }
  37. }