forked from GoekeLab/bioinformatics-workflows
-
Notifications
You must be signed in to change notification settings - Fork 1
/
example.wdl
71 lines (59 loc) · 983 Bytes
/
example.wdl
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
version 1.0
workflow wdl_poc {
input {
File reads1
File reads2
File ref_txome
}
call FastQC {
input:
reads1 = reads1,
reads2 = reads2,
}
call SalmonIndex {
input:
ref_txome = ref_txome,
}
call SalmonAlignQuant {
input:
reads1 = reads1,
reads2 = reads2,
index = SalmonIndex.index
}
}
task FastQC {
input {
File reads1
File reads2
}
command {
mkdir fastqc_res; fastqc --quiet "${reads1}" "${reads2}" --outdir fastqc_res
}
output {
File fastqc_res = "fastqc_res"
}
}
task SalmonIndex {
input {
File ref_txome
}
command {
salmon index -t "${ref_txome}" -i index
}
output {
File index = "index"
}
}
task SalmonAlignQuant {
input {
File reads1
File reads2
File index
}
command {
salmon quant -i "${index}" -l A -1 "${reads1}" -2 "${reads2}" --validateMappings -o quant
}
output {
File quant = "quant"
}
}