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.

4 年之前
4 年之前
4 年之前
4 年之前
4 年之前
1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950
  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. nt=$(nproc)
  14. spladder build --bams ${bam} \
  15. --annotation ${reference_gtf_file} \
  16. --ignore-mismatches \
  17. --confidence 2 \
  18. --readlen 150 \
  19. --parallel $nt \
  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. timeout: 864000
  29. }
  30. output {
  31. File fileList = "fileList.txt"
  32. Array[File] AS_gff = glob("${sample_id}/spladder_out/*.gff3")
  33. Array[File] AS_pickle = glob("${sample_id}/spladder_out/*.pickle")
  34. Array[File] AS_gz = glob("${sample_id}/spladder_out/*.gz")
  35. Array[File] AS_hdf5 = glob("${sample_id}/spladder_out/*.hdf5")
  36. Array[File] spladder = glob("${sample_id}/spladder_out/spladder/*")
  37. Array[File] tmp = glob("${sample_id}/spladder_out/tmp/*")
  38. }
  39. }