-
Notifications
You must be signed in to change notification settings - Fork 17
/
Copy path04-peakcall-se.cwl
167 lines (167 loc) · 4.85 KB
/
04-peakcall-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
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
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
#!/usr/bin/env cwl-runner
class: Workflow
cwlVersion: v1.0
doc: 'ATAC-seq 04 quantification - SE'
requirements:
- class: ScatterFeatureRequirement
- class: StepInputExpressionRequirement
- class: InlineJavascriptRequirement
inputs:
input_bam_files:
type: File[]
genome_effective_size:
default: hs
doc: Effective genome size used by MACS2. It can be numeric or a shortcuts:'hs' for human (2.7e9), 'mm' for mouse (1.87e9), 'ce' for C. elegans (9e7) and 'dm' for fruitfly (1.2e8), Default:hs
type: string
input_genome_sizes:
doc: Two column tab-delimited file with chromosome size information
type: File
as_narrowPeak_file:
doc: Definition narrowPeak file in AutoSql format (used in bedToBigBed)
type: File
nthreads:
default: 1
type: int
steps:
spp:
run: ../spp/spp.cwl
scatterMethod: dotproduct
scatter:
- input_bam
in:
input_bam: input_bam_files
nthreads: nthreads
savp:
valueFrom: ${return true}
out:
- output_spp_cross_corr
- output_spp_cross_corr_plot
extract-peak-frag-length:
run: ../spp/extract-best-frag-length.cwl
scatter: input_spp_txt_file
in:
input_spp_txt_file: spp/output_spp_cross_corr
out:
- output_best_frag_length
peak-calling:
run: ../peak_calling/macs2-callpeak.cwl
scatterMethod: dotproduct
scatter:
- treatment
in:
q:
valueFrom: ${return 0.1}
bdg:
valueFrom: ${return true}
treatment:
source: input_bam_files
valueFrom: $([self])
g: genome_effective_size
format:
valueFrom: BAM
extsize:
valueFrom: ${return 200}
nomodel:
valueFrom: ${return true}
shift:
valueFrom: ${return -100}
out:
- output_peak_file
- output_peak_summits_file
- output_ext_frag_bdg_file
- output_peak_xls_file
count-reads-filtered:
run: ../peak_calling/count-reads-after-filtering.cwl
in:
peak_xls_file: peak-calling/output_peak_xls_file
scatter: peak_xls_file
out:
- read_count_file
count-peaks:
run: ../utils/count-with-output-suffix.cwl
in:
output_suffix:
valueFrom: .peak_count.within_replicate.txt
input_file: peak-calling/output_peak_file
scatter: input_file
out:
- output_counts
filter-reads-in-peaks:
run: ../peak_calling/samtools-filter-in-bedfile.cwl
scatterMethod: dotproduct
scatter:
- input_bam_file
- input_bedfile
in:
input_bam_file: input_bam_files
input_bedfile: peak-calling/output_peak_file
out:
- filtered_file
extract-count-reads-in-peaks:
run: ../peak_calling/samtools-extract-number-mapped-reads.cwl
scatter: input_bam_file
in:
output_suffix:
valueFrom: .read_count.within_replicate.txt
input_bam_file: filter-reads-in-peaks/filtered_file
out:
- output_read_count
trunk-peak-score:
run: ../utils/trunk-peak-score.cwl
scatter: peaks
in:
peaks: peak-calling/output_peak_file
out:
- trunked_scores_peaks
peaks-bed-to-bigbed:
run: ../quant/bedToBigBed.cwl
in:
type:
valueFrom: bed6+4
as: as_narrowPeak_file
genome_sizes: input_genome_sizes
bed: trunk-peak-score/trunked_scores_peaks
scatter: bed
out:
- bigbed
outputs:
output_filtered_read_count_file:
doc: Filtered read count reported by MACS2
type: File[]
outputSource: count-reads-filtered/read_count_file
output_peak_summits_file:
doc: File containing peak summits
type: File[]
outputSource: peak-calling/output_peak_summits_file
output_peak_file:
doc: peakshift/phantomPeak results file
type: File[]
outputSource: peak-calling/output_peak_file
output_peak_count_within_replicate:
doc: Peak counts within replicate
type: File[]
outputSource: count-peaks/output_counts
output_spp_cross_corr_plot:
doc: peakshift/phantomPeak results file
type: File[]
outputSource: spp/output_spp_cross_corr_plot
output_spp_x_cross_corr:
doc: peakshift/phantomPeak results file
type: File[]
outputSource: spp/output_spp_cross_corr
output_extended_peak_file:
doc: peakshift/phantomPeak extended fragment results file
type: File[]
outputSource: peak-calling/output_ext_frag_bdg_file
output_peak_xls_file:
doc: Peak calling report file (*_peaks.xls file produced by MACS2)
type: File[]
outputSource: peak-calling/output_peak_xls_file
output_peak_bigbed_file:
doc: Peaks in bigBed format
type: File[]
outputSource: peaks-bed-to-bigbed/bigbed
output_read_in_peak_count_within_replicate:
doc: Reads peak counts within replicate
type: File[]
outputSource: extract-count-reads-in-peaks/output_read_count