-
Notifications
You must be signed in to change notification settings - Fork 17
/
03-quantification.cwl
145 lines (145 loc) · 3.77 KB
/
03-quantification.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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
class: Workflow
cwlVersion: v1.0
doc: DNase-seq 03 quantification
requirements:
- class: ScatterFeatureRequirement
- class: StepInputExpressionRequirement
- class: InlineJavascriptRequirement
inputs:
nthreads:
default: 1
type: int
input_pileup_bedgraphs:
type: File[]
input_bam_files:
type: File[]
input_peak_xls_files:
type: File[]
input_read_count_dedup_files:
type: File[]
input_genome_sizes:
type: File
steps:
bedsort_scaled_bdg:
run: ../quant/bedSort.cwl
in:
bed_file: scale-bedgraph/bedgraph_scaled
scatter: bed_file
out:
- bed_file_sorted
bedsort_genomecov:
run: ../quant/bedSort.cwl
in:
bed_file: bedtools_genomecov/output_bedfile
scatter: bed_file
out:
- bed_file_sorted
bdg2bw-extend-norm:
run: ../quant/bedGraphToBigWig.cwl
in:
output_suffix:
valueFrom: .fragment_extended.bw
genome_sizes: input_genome_sizes
bed_graph: bedsort_scaled_bdg/bed_file_sorted
scatter: bed_graph
out:
- output_bigwig
scale-bedgraph:
run: ../peak_calling/scale-bedgraph.cwl
in:
bedgraph_file: input_pileup_bedgraphs
read_count_file: input_read_count_dedup_files
scatterMethod: dotproduct
scatter:
- bedgraph_file
- read_count_file
out:
- bedgraph_scaled
bedsort_clipped_bedfile:
run: ../quant/bedSort.cwl
in:
bed_file: clip-off-chrom/bed_file_clipped
scatter: bed_file
out:
- bed_file_sorted
bdg2bw-raw:
run: ../quant/bedGraphToBigWig.cwl
in:
output_suffix:
valueFrom: .raw.bw
genome_sizes: input_genome_sizes
bed_graph: bedsort_genomecov/bed_file_sorted
scatter: bed_graph
out:
- output_bigwig
bedtools_genomecov:
run: ../map/bedtools-genomecov.cwl
in:
bg:
valueFrom: ${return true}
g: input_genome_sizes
ibam: input_bam_files
scatter: ibam
out:
- output_bedfile
bamcoverage:
run: ../quant/deeptools-bamcoverage.cwl
in:
numberOfProcessors: nthreads
extendReads:
valueFrom: ${return 200}
normalizeUsing:
valueFrom: RPKM
binSize:
valueFrom: ${return 50}
bam: input_bam_files
output_suffix:
valueFrom: .norm.bw
scatter: bam
out:
- output_bam_coverage
bdg2bw-extend:
run: ../quant/bedGraphToBigWig.cwl
in:
output_suffix:
valueFrom: .fragment_extended.bw
genome_sizes: input_genome_sizes
bed_graph: bedsort_clipped_bedfile/bed_file_sorted
scatter: bed_graph
out:
- output_bigwig
clip-off-chrom:
run: ../quant/bedClip.cwl
in:
bed_file: extend-reads/stdoutfile
genome_sizes: input_genome_sizes
scatter: bed_file
out:
- bed_file_clipped
extend-reads:
run: ../quant/bedtools-slop.cwl
in:
i: input_pileup_bedgraphs
b:
valueFrom: ${return 0}
g: input_genome_sizes
scatter: i
out:
- stdoutfile
outputs:
bigwig_raw_files:
doc: Raw reads bigWig (signal) files
type: File[]
outputSource: bdg2bw-raw/output_bigwig
bigwig_norm_files:
doc: Normalized reads bigWig (signal) files
type: File[]
outputSource: bamcoverage/output_bam_coverage
bigwig_extended_norm_files:
doc: Normalized fragment extended reads bigWig (signal) files
type: File[]
outputSource: bdg2bw-extend-norm/output_bigwig
bigwig_extended_files:
doc: Fragment extended reads bigWig (signal) files
type: File[]
outputSource: bdg2bw-extend/output_bigwig