|
kkonganti@1
|
1 params {
|
|
kkonganti@1
|
2 workflow_conceived_by = 'Kranti Konganti'
|
|
kkonganti@1
|
3 workflow_built_by = 'Kranti Konganti'
|
|
kkonganti@1
|
4 workflow_version = '0.5.0'
|
|
kkonganti@1
|
5 bcs_root_dbdir = '/tool/tool-data/cfsan_bettercallsal_db/latest'
|
|
kkonganti@1
|
6 bcs_db_mode = 'snp'
|
|
kkonganti@1
|
7 bcs_db_mode_index = (params.bcs_db_mode ==~ /snp/ ? '_cluster' : '_serotype')
|
|
kkonganti@1
|
8 bcs_thresholds = 'strict'
|
|
kkonganti@1
|
9 bcs_concat_pe = true
|
|
kkonganti@1
|
10 bbmerge_run = (params.bcs_concat_pe ? false : true)
|
|
kkonganti@1
|
11 bbmerge_adapters = false
|
|
kkonganti@1
|
12 bbmerge_reads = -1
|
|
kkonganti@1
|
13 bbmerge_ziplevel = 1
|
|
kkonganti@1
|
14 bbmerge_ordered = false
|
|
kkonganti@1
|
15 bbmerge_qtrim = true
|
|
kkonganti@1
|
16 bbmerge_qtrim2 = false
|
|
kkonganti@1
|
17 bbmerge_trimq = 10
|
|
kkonganti@1
|
18 bbmerge_minlength = 1
|
|
kkonganti@1
|
19 bbmerge_tbo = false
|
|
kkonganti@1
|
20 bbmerge_minavgquality = (params.bcs_thresholds ==~ /strict/ ? 30 : 20)
|
|
kkonganti@1
|
21 bbmerge_trimpolya = true
|
|
kkonganti@1
|
22 bbmerge_pfilter = (params.bcs_thresholds ==~ /strict/ ? 1 : false)
|
|
kkonganti@1
|
23 bbmerge_ouq = false
|
|
kkonganti@1
|
24 bbmerge_owq = true
|
|
kkonganti@1
|
25 bbmerge_strict = false
|
|
kkonganti@1
|
26 bbmerge_verystrict = false
|
|
kkonganti@1
|
27 bbmerge_ultrastrict = (params.bcs_thresholds ==~ /strict/ ? true : false)
|
|
kkonganti@1
|
28 bbmerge_maxstrict = false
|
|
kkonganti@1
|
29 bbmerge_loose = false
|
|
kkonganti@1
|
30 bbmerge_veryloose = false
|
|
kkonganti@1
|
31 bbmerge_ultraloose = false
|
|
kkonganti@1
|
32 bbmerge_maxloose = false
|
|
kkonganti@1
|
33 bbmerge_fast = false
|
|
kkonganti@1
|
34 bbmerge_k = 60
|
|
kkonganti@1
|
35 bbmerge_prealloc = true
|
|
kkonganti@1
|
36 fastp_run = true
|
|
kkonganti@1
|
37 fastp_failed_out = false
|
|
kkonganti@1
|
38 fastp_merged_out = false
|
|
kkonganti@1
|
39 fastp_overlapped_out = false
|
|
kkonganti@1
|
40 fastp_6 = false
|
|
kkonganti@1
|
41 fastp_reads_to_process = 0
|
|
kkonganti@1
|
42 fastp_fix_mgi_id = false
|
|
kkonganti@1
|
43 fastp_A = false
|
|
kkonganti@1
|
44 fastp_f = 0
|
|
kkonganti@1
|
45 fastp_t = 0
|
|
kkonganti@1
|
46 fastp_b = 0
|
|
kkonganti@1
|
47 fastp_F = 0
|
|
kkonganti@1
|
48 fastp_T = 0
|
|
kkonganti@1
|
49 fastp_B = 0
|
|
kkonganti@1
|
50 fastp_dedup = true
|
|
kkonganti@1
|
51 fastp_dup_calc_accuracy = 6
|
|
kkonganti@1
|
52 fastp_poly_g_min_len = 10
|
|
kkonganti@1
|
53 fastp_G = true
|
|
kkonganti@1
|
54 fastp_x = false
|
|
kkonganti@1
|
55 fastp_poly_x_min_len = 10
|
|
kkonganti@1
|
56 fastp_cut_front = true
|
|
kkonganti@1
|
57 fastp_cut_tail = false
|
|
kkonganti@1
|
58 fastp_cut_right = true
|
|
kkonganti@1
|
59 fastp_W = 20
|
|
kkonganti@1
|
60 fastp_M = 30
|
|
kkonganti@1
|
61 fastp_q = 30
|
|
kkonganti@1
|
62 fastp_u = 40
|
|
kkonganti@1
|
63 fastp_n = 5
|
|
kkonganti@1
|
64 fastp_e = 0
|
|
kkonganti@1
|
65 fastp_l = 35
|
|
kkonganti@1
|
66 fastp_max_len = 0
|
|
kkonganti@1
|
67 fastp_y = true
|
|
kkonganti@1
|
68 fastp_Y = 30
|
|
kkonganti@1
|
69 fastp_U = false
|
|
kkonganti@1
|
70 fastp_umi_loc = false
|
|
kkonganti@1
|
71 fastp_umi_len = false
|
|
kkonganti@1
|
72 fastp_umi_prefix = false
|
|
kkonganti@1
|
73 fastp_umi_skip = false
|
|
kkonganti@1
|
74 fastp_p = true
|
|
kkonganti@1
|
75 fastp_P = 20
|
|
kkonganti@1
|
76 fastp_use_custom_adapters = false
|
|
kkonganti@1
|
77 fastp_adapter_fasta = (params.fastp_use_custom_adapters ? "${projectDir}"
|
|
kkonganti@1
|
78 + File.separator
|
|
kkonganti@1
|
79 + 'assets'
|
|
kkonganti@1
|
80 + File.separator
|
|
kkonganti@1
|
81 + 'adaptors.fa' : false)
|
|
kkonganti@1
|
82 mashscreen_run = true
|
|
kkonganti@1
|
83 mash_sketch = (params.bcs_root_dbdir
|
|
kkonganti@1
|
84 + File.separator
|
|
kkonganti@1
|
85 + 'mash_sketch'
|
|
kkonganti@1
|
86 + File.separator
|
|
kkonganti@1
|
87 + 'msh.k21.1000h.'
|
|
kkonganti@1
|
88 + params.bcs_db_mode
|
|
kkonganti@1
|
89 + '.msh')
|
|
kkonganti@1
|
90 mashscreen_w = false
|
|
kkonganti@1
|
91 mashscreen_i = false
|
|
kkonganti@1
|
92 mashscreen_v = false
|
|
kkonganti@1
|
93 tuspy_run = true
|
|
kkonganti@1
|
94 tuspy_s = false
|
|
kkonganti@1
|
95 tuspy_m = false
|
|
kkonganti@1
|
96 tuspy_ps = (params.bcs_root_dbdir
|
|
kkonganti@1
|
97 + File.separator
|
|
kkonganti@1
|
98 + 'index_metadata'
|
|
kkonganti@1
|
99 + File.separator
|
|
kkonganti@1
|
100 + 'per_'
|
|
kkonganti@1
|
101 + params.bcs_db_mode
|
|
kkonganti@1
|
102 + params.bcs_db_mode_index
|
|
kkonganti@1
|
103 + '.ACC2SERO.pickle')
|
|
kkonganti@1
|
104 tuspy_gd = (params.bcs_root_dbdir
|
|
kkonganti@1
|
105 + File.separator
|
|
kkonganti@1
|
106 + 'scaffold_genomes')
|
|
kkonganti@1
|
107 tuspy_gds = '_scaffolded_genomic.fna.gz'
|
|
kkonganti@1
|
108 tuspy_n = 10
|
|
kkonganti@1
|
109 sourmashsketch_run = true
|
|
kkonganti@1
|
110 sourmashsketch_mode = 'dna'
|
|
kkonganti@1
|
111 sourmashsketch_file = false
|
|
kkonganti@1
|
112 sourmashsketch_f = false
|
|
kkonganti@1
|
113 sourmashsketch_merge = false
|
|
kkonganti@1
|
114 sourmashsketch_name = false
|
|
kkonganti@1
|
115 sourmashsketch_p = 'abund,scaled=1000,k=51,k=61,k=71'
|
|
kkonganti@1
|
116 sourmashsketch_singleton = true
|
|
kkonganti@1
|
117 sourmashsketch_randomize = false
|
|
kkonganti@1
|
118 sourmashgather_run = (params.sourmashsketch_run ?: false)
|
|
kkonganti@1
|
119 sourmashgather_n = false
|
|
kkonganti@1
|
120 sourmashgather_thr_bp = false
|
|
kkonganti@1
|
121 sourmashgather_ignoreabn = false
|
|
kkonganti@1
|
122 sourmashgather_prefetch = false
|
|
kkonganti@1
|
123 sourmashgather_noprefetch = false
|
|
kkonganti@1
|
124 sourmashgather_ani_ci = true
|
|
kkonganti@1
|
125 sourmashgather_k = 71
|
|
kkonganti@1
|
126 sourmashgather_protein = false
|
|
kkonganti@1
|
127 sourmashgather_noprotein = false
|
|
kkonganti@1
|
128 sourmashgather_dayhoff = false
|
|
kkonganti@1
|
129 sourmashgather_nodayhoff = false
|
|
kkonganti@1
|
130 sourmashgather_hp = false
|
|
kkonganti@1
|
131 sourmashgather_nohp = false
|
|
kkonganti@1
|
132 sourmashgather_dna = true
|
|
kkonganti@1
|
133 sourmashgather_nodna = false
|
|
kkonganti@1
|
134 sourmashgather_scaled = false
|
|
kkonganti@1
|
135 sourmashgather_inc_pat = false
|
|
kkonganti@1
|
136 sourmashgather_exc_pat = false
|
|
kkonganti@1
|
137 sourmashsearch_run = (params.sourmashsketch_run && !params.sourmashgather_run ?: false)
|
|
kkonganti@1
|
138 sourmashsearch_n = false
|
|
kkonganti@1
|
139 sourmashsearch_thr = '0'
|
|
kkonganti@1
|
140 sourmashsearch_contain = false
|
|
kkonganti@1
|
141 sourmashsearch_maxcontain = true
|
|
kkonganti@1
|
142 sourmashsearch_ignoreabn = true
|
|
kkonganti@1
|
143 sourmashsearch_ani_ci = false
|
|
kkonganti@1
|
144 sourmashsearch_k = 71
|
|
kkonganti@1
|
145 sourmashsearch_protein = false
|
|
kkonganti@1
|
146 sourmashsearch_noprotein = false
|
|
kkonganti@1
|
147 sourmashsearch_dayhoff = false
|
|
kkonganti@1
|
148 sourmashsearch_nodayhoff = false
|
|
kkonganti@1
|
149 sourmashsearch_hp = false
|
|
kkonganti@1
|
150 sourmashsearch_nohp = false
|
|
kkonganti@1
|
151 sourmashsearch_dna = true
|
|
kkonganti@1
|
152 sourmashsearch_nodna = false
|
|
kkonganti@1
|
153 sourmashsearch_scaled = false
|
|
kkonganti@1
|
154 sourmashsearch_inc_pat = false
|
|
kkonganti@1
|
155 sourmashsearch_exc_pat = false
|
|
kkonganti@1
|
156 sfhpy_run = true
|
|
kkonganti@1
|
157 sfhpy_fcn = (params.sourmashsearch_run ? 'similarity' : 'f_match')
|
|
kkonganti@1
|
158 sfhpy_fcv = (params.sourmashsearch_run ? 0.9 : 0.1)
|
|
kkonganti@1
|
159 sfhpy_gt = true
|
|
kkonganti@1
|
160 sfhpy_lt = false
|
|
kkonganti@1
|
161 kmaindex_run = true
|
|
kkonganti@1
|
162 kmaindex_t_db = false
|
|
kkonganti@1
|
163 kmaindex_k = 31
|
|
kkonganti@1
|
164 kmaindex_m = false
|
|
kkonganti@1
|
165 kmaindex_hc = false
|
|
kkonganti@1
|
166 kmaindex_ML = false
|
|
kkonganti@1
|
167 kmaindex_ME = false
|
|
kkonganti@1
|
168 kmaindex_Sparse = false
|
|
kkonganti@1
|
169 kmaindex_ht = false
|
|
kkonganti@1
|
170 kmaindex_hq = false
|
|
kkonganti@1
|
171 kmaindex_and = false
|
|
kkonganti@1
|
172 kmaindex_nbp = false
|
|
kkonganti@1
|
173 kmaalign_run = true
|
|
kkonganti@1
|
174 kmaalign_ignorequals = false
|
|
kkonganti@1
|
175 kmaalign_int = false
|
|
kkonganti@1
|
176 kmaalign_ef = false
|
|
kkonganti@1
|
177 kmaalign_vcf = false
|
|
kkonganti@1
|
178 kmaalign_sam = false
|
|
kkonganti@1
|
179 kmaalign_nc = true
|
|
kkonganti@1
|
180 kmaalign_na = true
|
|
kkonganti@1
|
181 kmaalign_nf = true
|
|
kkonganti@1
|
182 kmaalign_a = false
|
|
kkonganti@1
|
183 kmaalign_and = true
|
|
kkonganti@1
|
184 kmaalign_oa = false
|
|
kkonganti@1
|
185 kmaalign_bc = false
|
|
kkonganti@1
|
186 kmaalign_bcNano = false
|
|
kkonganti@1
|
187 kmaalign_bcd = false
|
|
kkonganti@1
|
188 kmaalign_bcg = false
|
|
kkonganti@1
|
189 kmaalign_ID = (params.bcs_thresholds ==~ /strict/ ? 10.0 : 5.0)
|
|
kkonganti@1
|
190 kmaalign_md = false
|
|
kkonganti@1
|
191 kmaalign_dense = false
|
|
kkonganti@1
|
192 kmaalign_ref_fsa = false
|
|
kkonganti@1
|
193 kmaalign_Mt1 = false
|
|
kkonganti@1
|
194 kmaalign_1t1 = false
|
|
kkonganti@1
|
195 kmaalign_mrs = (params.bcs_thresholds ==~ /strict/ ? 0.99 : 0.90)
|
|
kkonganti@1
|
196 kmaalign_mrc = (params.bcs_thresholds ==~ /strict/ ? 0.99 : 0.90)
|
|
kkonganti@1
|
197 kmaalign_mp = (params.bcs_thresholds ==~ /strict/ ? 30 : 20)
|
|
kkonganti@1
|
198 kmaalign_eq = (params.bcs_thresholds ==~ /strict/ ? 30 : 20)
|
|
kkonganti@1
|
199 kmaalign_mp = (params.kmaalign_ignorequals ? 0 : params.kmaalign_mp)
|
|
kkonganti@1
|
200 kmaalign_eq = (params.kmaalign_ignorequals ? 0 : params.kmaalign_eq)
|
|
kkonganti@1
|
201 kmaalign_mq = false
|
|
kkonganti@1
|
202 kmaalign_5p = false
|
|
kkonganti@1
|
203 kmaalign_3p = false
|
|
kkonganti@1
|
204 kmaalign_apm = false
|
|
kkonganti@1
|
205 kmaalign_cge = false
|
|
kkonganti@1
|
206 kmaalign_mrs = (params.sourmashgather_run ? false : params.kmaalign_mrs)
|
|
kkonganti@1
|
207 kmaalign_and = (params.sourmashgather_run ? false : true)
|
|
kkonganti@1
|
208 kmaalign_ID = (params.sourmashgather_run ? false : params.kmaalign_ID)
|
|
kkonganti@1
|
209 salmonidx_run = true
|
|
kkonganti@1
|
210 salmonidx_k = false
|
|
kkonganti@1
|
211 salmonidx_gencode = false
|
|
kkonganti@1
|
212 salmonidx_features = false
|
|
kkonganti@1
|
213 salmonidx_keepDuplicates = false
|
|
kkonganti@1
|
214 salmonidx_keepFixedFasta = false
|
|
kkonganti@1
|
215 salmonidx_filterSize = false
|
|
kkonganti@1
|
216 salmonidx_sparse = false
|
|
kkonganti@1
|
217 salmonidx_n = false
|
|
kkonganti@1
|
218 salmonidx_decoys = false
|
|
kkonganti@1
|
219 salmonalign_libtype = 'SF'
|
|
kkonganti@1
|
220 gsrpy_run = true
|
|
kkonganti@1
|
221 gsrpy_url = true
|
|
kkonganti@1
|
222 gsrpy_snp_clus_metadata = (params.bcs_root_dbdir
|
|
kkonganti@1
|
223 + File.separator
|
|
kkonganti@1
|
224 + 'download_pdg_metadata'
|
|
kkonganti@1
|
225 + File.separator
|
|
kkonganti@1
|
226 + getPDGRelease(params.bcs_root_dbdir)
|
|
kkonganti@1
|
227 + File.separator
|
|
kkonganti@1
|
228 + getPDGRelease(params.bcs_root_dbdir)
|
|
kkonganti@1
|
229 + '.reference_target.cluster_list.tsv')
|
|
kkonganti@1
|
230 sourmashcompare_mode = 'ani,containment'
|
|
kkonganti@1
|
231 fq_filter_by_len = 0
|
|
kkonganti@1
|
232 fq_suffix = '.fastq.gz'
|
|
kkonganti@1
|
233 fq2_suffix = '_R2_001.fastq.gz'
|
|
kkonganti@1
|
234 }
|
|
kkonganti@1
|
235
|
|
kkonganti@1
|
236 /*
|
|
kkonganti@1
|
237 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
|
kkonganti@1
|
238 METHOD TO GET PDG RELEASE AND SET REST OF THE METADATA FILES
|
|
kkonganti@1
|
239 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
|
kkonganti@1
|
240 */
|
|
kkonganti@1
|
241
|
|
kkonganti@1
|
242 def getPDGRelease(db_path) {
|
|
kkonganti@1
|
243 db_path_obj = new File( db_path )
|
|
kkonganti@1
|
244 def pdg_release = []
|
|
kkonganti@1
|
245
|
|
kkonganti@1
|
246 db_path_obj.eachDirRecurse {
|
|
kkonganti@1
|
247 it.toString() =~ /PDG/ \
|
|
kkonganti@1
|
248 ? pdg_release << it.toString().find(/PDG\d+\.\d+/) \
|
|
kkonganti@1
|
249 : null
|
|
kkonganti@1
|
250 }
|
|
kkonganti@1
|
251
|
|
kkonganti@1
|
252 return pdg_release.unique().join('').toString()
|
|
kkonganti@1
|
253 }
|
