Mercurial > repos > kkonganti > cfsan_bettercallsal

comparison 0.5.0/workflows/conf/bettercallsal.config @ 1:365849f031fd

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