-
Notifications
You must be signed in to change notification settings - Fork 37
/
Copy pathalign_MT.config
108 lines (90 loc) · 3.76 KB
/
align_MT.config
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
/*
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Config file for defining DSL2 per module options and publishing paths
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Available keys to override module options:
ext.args = Additional arguments appended to command in module.
ext.args2 = Second set of arguments appended to command in module (multi-tool modules).
ext.args3 = Third set of arguments appended to command in module (multi-tool modules).
ext.prefix = File name prefix for output files.
ext.when = Conditional clause
----------------------------------------------------------------------------------------
*/
//
// Mitochondria alignment options
//
process {
if (params.aligner.equals("bwamem2")) {
withName: '.*ALIGN_MT:BWAMEM2_MEM_MT' {
ext.args = { "-M -K 100000000 -R ${meta.read_group}" }
ext.when = { !(params.analysis_type.equals("wes")) }
ext.prefix = { "${meta.id}.sorted" }
}
}
if (params.aligner.equals("sentieon")) {
withName: '.*ALIGN_MT:SENTIEON_BWAMEM_MT' {
ext.args = { "-M -K 10000000 -R ${meta.read_group}" }
ext.when = { !(params.analysis_type.equals("wes"))}
ext.prefix = { "${meta.id}.sorted" }
}
}
withName: '.*ALIGN_MT:GATK4_MERGEBAMALIGNMENT_MT' {
ext.args = '--CREATE_INDEX true --MAX_GAPS -1 --SORT_ORDER queryname --INCLUDE_SECONDARY_ALIGNMENTS false --PAIRED_RUN false --VALIDATION_STRINGENCY LENIENT'
ext.prefix = { "${meta.id}_merged" }
}
withName: '.*ALIGN_MT:PICARD_ADDORREPLACEREADGROUPS_MT' {
ext.args = { [
"--VALIDATION_STRINGENCY LENIENT",
"--RGLB lib",
"--RGPL ILLUMINA",
"--RGPU barcode",
"--RGSM ${meta.id}"
].join(' ').trim() }
}
withName: '.*ALIGN_MT:PICARD_MARKDUPLICATES_MT' {
ext.args = '--VALIDATION_STRINGENCY LENIENT --CREATE_INDEX true'
ext.prefix = { "${meta.id}_markduplicates" }
}
withName: '.*ALIGN_MT:SAMTOOLS_SORT_MT' {
ext.prefix = { "${meta.id}_sorted" }
}
}
//
// Shifted mitochondria alignment options
//
process {
if (params.aligner.equals("bwamem2")) {
withName: '.*ALIGN_MT_SHIFT:BWAMEM2_MEM_MT' {
ext.args = { "-M -K 100000000 -R ${meta.read_group}" }
ext.when = { !(params.analysis_type.equals("wes")) }
ext.prefix = { "${meta.id}.sorted" }
}
}
if (params.aligner.equals("sentieon")) {
withName: '.*ALIGN_MT_SHIFT:SENTIEON_BWAMEM_MT' {
ext.args = { "-M -K 10000000 -R ${meta.read_group}" }
ext.when = { !(params.analysis_type.equals("wes"))}
ext.prefix = { "${meta.id}.sorted" }
}
}
withName: '.*ALIGN_MT_SHIFT:GATK4_MERGEBAMALIGNMENT_MT' {
ext.args = '--CREATE_INDEX true --MAX_GAPS -1 --SORT_ORDER queryname --INCLUDE_SECONDARY_ALIGNMENTS false --PAIRED_RUN false --VALIDATION_STRINGENCY LENIENT'
ext.prefix = { "${meta.id}_merged_shifted" }
}
withName: '.*ALIGN_MT_SHIFT:PICARD_ADDORREPLACEREADGROUPS_MT' {
ext.args = { [
"--VALIDATION_STRINGENCY LENIENT",
"--RGLB lib",
"--RGPL ${params.platform}",
"--RGPU barcode",
"--RGSM ${meta.id}"
].join(' ').trim() }
}
withName: '.*ALIGN_MT_SHIFT:PICARD_MARKDUPLICATES_MT' {
ext.args = '--VALIDATION_STRINGENCY LENIENT --CREATE_INDEX true'
ext.prefix = { "${meta.id}_markduplicates_shifted" }
}
withName: '.*ALIGN_MT_SHIFT:SAMTOOLS_SORT_MT' {
ext.prefix = { "${meta.id}_sorted_shifted" }
}
}