/
04-quantification-se-revstranded.cwl
154 lines (154 loc) · 3.97 KB
/
04-quantification-se-revstranded.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
#!/usr/bin/env cwl-runner
class: Workflow
cwlVersion: v1.0
doc: RNA-seq 04 quantification
requirements:
- class: ScatterFeatureRequirement
- class: StepInputExpressionRequirement
- class: InlineJavascriptRequirement
- class: SubworkflowFeatureRequirement
inputs:
input_bam_files:
type: File[]
nthreads:
default: 1
type: int
rsem_reference_files:
doc: RSEM genome reference files - generated with the rsem-prepare-reference command
type: Directory
annotation_file:
doc: GTF annotation file
type: File
input_transcripts_bam_files:
type: File[]
input_genome_sizes:
type: File
outputs:
featurecounts_counts:
doc: Normalized fragment extended reads bigWig (signal) files
type: File[]
outputSource: featurecounts/output_files
rsem_isoforms_files:
doc: RSEM isoforms files
type: File[]
outputSource: rsem-calc-expr/isoforms
rsem_genes_files:
doc: RSEM genes files
type: File[]
outputSource: rsem-calc-expr/genes
bw_raw_files:
doc: Raw bigWig files.
type: File[]
outputSource: bdg2bw-raw/output_bigwig
bw_norm_files:
doc: Normalized by RPKM bigWig files.
type: File[]
outputSource: bamcoverage/output_bam_coverage
steps:
basename:
run: ../utils/basename.cwl
in:
file_path:
source: input_bam_files
valueFrom: $(self.basename)
sep:
valueFrom: .Aligned.out.sorted
scatter: file_path
out:
- basename
featurecounts:
run: ../quant/subread-featurecounts.cwl
in:
B:
valueFrom: ${return true}
g:
valueFrom: gene_id
output_filename:
source: basename/basename
valueFrom: $(self + ".featurecounts.counts.txt")
s:
valueFrom: ${return 2}
t:
valueFrom: exon
annotation_file: annotation_file
T: nthreads
input_files:
source: input_bam_files
valueFrom: ${if (Array.isArray(self)) return self; return [self]; }
scatterMethod: dotproduct
scatter:
- input_files
- output_filename
out:
- output_files
rsem-calc-expr:
run: ../quant/rsem-calculate-expression.cwl
in:
reference_name:
source: rsem_reference_files
valueFrom: |
${
var trans_file_str = self.listing.map(function(e){return e.location}).filter(function(e){return e.match(/\.transcripts\.fa$/)})[0];
return trans_file_str.match(/.*[\\\\\\/](.*)\.transcripts\.fa$/)[1];
}
reference_files: rsem_reference_files
no-bam-output:
valueFrom: ${return true}
quiet:
valueFrom: ${return true}
seed:
valueFrom: ${return 1234}
sample_name:
source: basename/basename
valueFrom: $(self + ".rsem")
bam: input_transcripts_bam_files
num-threads: nthreads
scatterMethod: dotproduct
scatter:
- bam
- sample_name
out:
- isoforms
- genes
- rsem_stat
bedsort_genomecov:
run: ../quant/bedSort.cwl
scatter: bed_file
in:
bed_file: bedtools_genomecov/output_bedfile
out:
- bed_file_sorted
bdg2bw-raw:
run: ../quant/bedGraphToBigWig.cwl
scatter: bed_graph
in:
output_suffix:
valueFrom: .raw.bw
genome_sizes: input_genome_sizes
bed_graph: bedsort_genomecov/bed_file_sorted
out:
- output_bigwig
bedtools_genomecov:
run: ../map/bedtools-genomecov.cwl
scatter: ibam
in:
bg:
valueFrom: ${return true}
g: input_genome_sizes
ibam: input_bam_files
out:
- output_bedfile
bamcoverage:
run: ../quant/deeptools-bamcoverage.cwl
scatter: bam
in:
binSize:
valueFrom: ${return 1}
numberOfProcessors: nthreads
bam: input_bam_files
output_suffix:
valueFrom: .norm.bw
normalizeUsing:
valueFrom: RPKM
out:
- output_bam_coverage