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.

111 line
2.6KB

  1. import "./tasks/fastp.wdl" as fastp
  2. import "./tasks/hisat2.wdl" as hisat2
  3. import "./tasks/samtools.wdl" as samtools
  4. import "./tasks/stringtie.wdl" as stringtie
  5. import "./tasks/fastqc.wdl" as fastqc
  6. import "./tasks/fastqscreen.wdl" as fastqscreen
  7. import "./tasks/qualimapBAMqc.wdl" as qualimapBAMqc
  8. import "./tasks/qualimapRNAseq.wdl" as qualimapRNAseq
  9. import "./tasks/ballgown.wdl" as ballgown
  10. import "./tasks/count.wdl" as count
  11. import "./tasks/MD5.wdl" as MD5
  12. workflow {{ project_name }} {
  13. File read1
  14. File read2
  15. File idx
  16. File screen_ref_dir
  17. File fastq_screen_conf
  18. File gtf
  19. String sample_id
  20. String fastp_docker
  21. String adapter_sequence
  22. String adapter_sequence_r2
  23. String fastp_cluster
  24. String umi_loc
  25. String idx_prefix
  26. String pen_intronlen
  27. String fastqc_cluster_config
  28. String fastqc_docker
  29. String fastqscreen_docker
  30. String fastqscreen_cluster_config
  31. String hisat2_docker
  32. String hisat2_cluster
  33. String qualimapBAMqc_docker
  34. String qualimapBAMqc_cluster_config
  35. String qualimapRNAseq_docker
  36. String qualimapRNAseq_cluster_config
  37. String samtools_docker
  38. String samtools_cluster
  39. String stringtie_docker
  40. String stringtie_cluster
  41. String multiqc_cluster_config
  42. String multiqc_docker
  43. Int multiqc_disk_size
  44. Int trim_front1
  45. Int trim_tail1
  46. Int max_len1
  47. Int trim_front2
  48. Int trim_tail2
  49. Int max_len2
  50. Int disable_adapter_trimming
  51. Int length_required
  52. Int umi_len
  53. Int UMI
  54. Int qualified_quality_phred
  55. Int length_required1
  56. Int disable_quality_filtering
  57. Int pen_cansplice
  58. Int pen_noncansplice
  59. Int min_intronlen
  60. Int max_intronlen
  61. Int maxins
  62. Int minins
  63. Int fastqc_disk_size
  64. Int fastqscreen_disk_size
  65. Int qualimapBAMqc_disk_size
  66. Int qualimapRNAseq_disk_size
  67. Int insert_size
  68. Int minimum_length_allowed_for_the_predicted_transcripts
  69. Int Junctions_no_spliced_reads
  70. Int count_length
  71. Float minimum_isoform_abundance
  72. Float maximum_fraction_of_muliplelocationmapped_reads
  73. String ballgown_docker
  74. String ballgown_cluster
  75. String disk_size
  76. String count_docker
  77. String count_cluster
  78. String docker
  79. String cluster_config
  80. call fastqc.fastqc as fastqc {
  81. input:
  82. read1 = read1,
  83. read2 = read2,
  84. docker = fastqc_docker,
  85. cluster_config = fastqc_cluster_config,
  86. disk_size = fastqc_disk_size
  87. }
  88. call fastqscreen.fastq_screen as fastqscreen {
  89. input:
  90. read1 = read1,
  91. read2 = read2,
  92. screen_ref_dir = screen_ref_dir,
  93. fastq_screen_conf = fastq_screen_conf,
  94. docker = fastqscreen_docker,
  95. cluster_config = fastqscreen_cluster_config,
  96. disk_size = fastqscreen_disk_size
  97. }
  98. call MD5.MD5 as MD5 {
  99. input:
  100. file1 = fastqc.read1_zip,
  101. file2 = fastqscreen.txt1,
  102. docker = docker,
  103. disk_size = disk_size,
  104. cluster_config=cluster_config
  105. }
  106. }