from fastq to bam 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.

102 lignes
2.2KB

  1. import "./tasks/fastp.wdl" as fastp
  2. import "./tasks/hisat2.wdl" as hisat2
  3. import "./tasks/samtools.wdl" as samtools
  4. workflow {{ project_name }} {
  5. String sample_id
  6. File read1
  7. File read2
  8. String adapter_sequence
  9. String adapter_sequence_r2
  10. String fastp_docker
  11. String fastp_cluster
  12. String disk_size
  13. String umi_loc
  14. Int trim_tail1
  15. Int max_len1
  16. Int trim_front2
  17. Int trim_tail2
  18. Int max_len2
  19. Int disable_adapter_trimming
  20. Int length_required
  21. Int umi_len
  22. Int UMI
  23. Int qualified_quality_phred
  24. Int length_required1
  25. Int disable_quality_filtering
  26. File idx
  27. String idx_prefix
  28. String pen_intronlen
  29. String hisat2_docker
  30. String hisat2_cluster
  31. Int trim_front1
  32. Int pen_cansplice
  33. Int pen_noncansplice
  34. Int min_intronlen
  35. Int max_intronlen
  36. Int maxins
  37. Int minins
  38. Int insert_size
  39. String samtools_docker
  40. String samtools_cluster
  41. call fastp.fastp as fastp {
  42. input:
  43. sample_id=sample_id,
  44. read1 = read1,
  45. read2 = read2,
  46. docker = fastp_docker,
  47. cluster = fastp_cluster,
  48. disk_size = disk_size,
  49. adapter_sequence = adapter_sequence,
  50. adapter_sequence_r2 = adapter_sequence_r2,
  51. umi_loc = umi_loc,
  52. trim_front1 = trim_front1,
  53. trim_tail1 = trim_tail1,
  54. max_len1 = max_len1,
  55. trim_front2 = trim_front2,
  56. trim_tail2 = trim_tail2,
  57. max_len2 = max_len2,
  58. disable_adapter_trimming = disable_adapter_trimming,
  59. length_required = length_required,
  60. umi_len = umi_len,
  61. UMI = UMI,
  62. qualified_quality_phred = qualified_quality_phred,
  63. length_required1 = length_required1,
  64. disable_quality_filtering = disable_quality_filtering
  65. }
  66. call hisat2.hisat2 as hisat2 {
  67. input:
  68. sample_id = sample_id,
  69. idx = idx,
  70. idx_prefix = idx_prefix,
  71. Trim_R1 = fastp.Trim_R1,
  72. Trim_R2 = fastp.Trim_R2,
  73. docker = hisat2_docker,
  74. cluster = hisat2_cluster,
  75. disk_size = disk_size,
  76. pen_intronlen = pen_intronlen,
  77. pen_cansplice = pen_cansplice,
  78. pen_noncansplice = pen_noncansplice,
  79. min_intronlen = min_intronlen,
  80. max_intronlen = max_intronlen,
  81. maxins = maxins,
  82. minins = minins
  83. }
  84. call samtools.samtools as samtools {
  85. input:
  86. sample_id = sample_id,
  87. sam = hisat2.sam,
  88. docker = samtools_docker,
  89. cluster = samtools_cluster,
  90. disk_size = disk_size,
  91. insert_size = insert_size
  92. }
  93. }