-
Notifications
You must be signed in to change notification settings - Fork 17
/
02-trim-se.cwl
91 lines (91 loc) · 2.49 KB
/
02-trim-se.cwl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
#!/usr/bin/env cwl-runner
class: Workflow
cwlVersion: v1.0
doc: 'RNA-seq 02 trimming - reads: SE'
requirements:
- class: ScatterFeatureRequirement
- class: StepInputExpressionRequirement
- class: InlineJavascriptRequirement
inputs:
input_fastq_read1_files:
doc: Input read 1 fastq files
type: File[]
input_read1_adapters_files:
doc: Input read 1 adapters files
type: File[]
nthreads:
default: 1
doc: Number of threads
type: int
trimmomatic_java_opts:
doc: JVM arguments should be a quoted, space separated list
type: string?
quality_score:
default: -phred33
type: string
trimmomatic_jar_path:
default: /usr/share/java/trimmomatic.jar
doc: Trimmomatic Java jar file
type: string
outputs:
output_data_fastq_read1_trimmed_files:
doc: Trimmed fastq files for paired read 1
type: File[]
outputSource: trimmomatic/output_read1_trimmed_file
output_trimmed_read1_fastq_read_count:
doc: Trimmed read counts of paired read 1 fastq files
type: File[]
outputSource: count_fastq_reads_read1/output_read_count
steps:
trimmomatic:
run: ../trimmomatic/trimmomatic.cwl
scatterMethod: dotproduct
scatter:
- input_read1_fastq_file
- input_adapters_file
in:
input_read1_fastq_file: input_fastq_read1_files
phred:
valueFrom: '33'
nthreads: nthreads
minlen:
valueFrom: ${return 15}
input_adapters_file: input_read1_adapters_files
leading:
valueFrom: ${return 3}
slidingwindow:
valueFrom: 4:20
illuminaclip:
valueFrom: 2:30:15
end_mode:
valueFrom: SE
java_opts: trimmomatic_java_opts
trailing:
valueFrom: ${return 3}
trimmomatic_jar_path: trimmomatic_jar_path
out:
- output_read1_trimmed_file
extract_basename_read1:
run: ../utils/basename.cwl
scatter: file_path
in:
file_path:
source: trimmomatic/output_read1_trimmed_file
valueFrom: $(self.basename)
sep:
valueFrom: '(\.fastq.gz|\.fastq)'
do_not_escape_sep:
valueFrom: ${return true}
out:
- basename
count_fastq_reads_read1:
run: ../utils/count-fastq-reads.cwl
scatterMethod: dotproduct
scatter:
- input_fastq_file
- input_basename
in:
input_basename: extract_basename_read1/basename
input_fastq_file: trimmomatic/output_read1_trimmed_file
out:
- output_read_count