kkonganti@0: params {
kkonganti@0:     workflow_conceived_by = 'Kranti Konganti'
kkonganti@0:     workflow_built_by = 'Kranti Konganti'
kkonganti@0:     workflow_version = '0.1.0'
kkonganti@0:     ref_acc =  'GCF_003516125'
kkonganti@4:     cronology_root_dbdir = '/tool/tool-data/cfsan_cronology_db/0.1.0'
kkonganti@0:     dpubmlstpy_run = true
kkonganti@0:     dpubmlstpy_org = 'cronobacter'
kkonganti@0:     dpubmlstpy_mlsts = 'schemes/1'
kkonganti@0:     dpubmlstpy_profile = false
kkonganti@0:     dpubmlstpy_loci = false
kkonganti@0:     dpubmlstpy_suffix = false
kkonganti@0:     dpubmlstpy_akey = false
kkonganti@0:     dpubmlstpy_id = false
kkonganti@0:     fastp_run = true
kkonganti@0:     fastp_failed_out = false
kkonganti@0:     fastp_merged_out = false
kkonganti@0:     fastp_overlapped_out = false
kkonganti@0:     fastp_6 = false
kkonganti@0:     fastp_reads_to_process = 0 
kkonganti@0:     fastp_fix_mgi_id = false
kkonganti@0:     fastp_A = false
kkonganti@0:     fastp_f = 0
kkonganti@0:     fastp_t = 0
kkonganti@0:     fastp_b = 0
kkonganti@0:     fastp_F = 0
kkonganti@0:     fastp_T = 0
kkonganti@0:     fastp_B = 0
kkonganti@0:     fastp_dedup = true
kkonganti@0:     fastp_dup_calc_accuracy = 6
kkonganti@0:     fastp_poly_g_min_len = 10
kkonganti@0:     fastp_G = true
kkonganti@0:     fastp_x = false
kkonganti@0:     fastp_poly_x_min_len = 10 
kkonganti@0:     fastp_cut_front = true
kkonganti@0:     fastp_cut_tail = false
kkonganti@0:     fastp_cut_right = true
kkonganti@0:     fastp_W = 20
kkonganti@0:     fastp_M = 30
kkonganti@0:     fastp_q = 30
kkonganti@0:     fastp_u = 40
kkonganti@0:     fastp_n = 5
kkonganti@0:     fastp_e = 0
kkonganti@0:     fastp_l = 35
kkonganti@0:     fastp_max_len = 0
kkonganti@0:     fastp_y = true
kkonganti@0:     fastp_Y = 30
kkonganti@0:     fastp_U = false
kkonganti@0:     fastp_umi_loc = false
kkonganti@0:     fastp_umi_len = false
kkonganti@0:     fastp_umi_prefix = false
kkonganti@0:     fastp_umi_skip = false
kkonganti@0:     fastp_p = true
kkonganti@0:     fastp_P = 20
kkonganti@0:     fastp_use_custom_adapters = false
kkonganti@0:     fastp_adapter_fasta = (params.fastp_use_custom_adapters ? "${projectDir}"
kkonganti@0:         + File.separator
kkonganti@0:         + 'assets'
kkonganti@0:         + File.separator
kkonganti@0:         + 'adaptors.fa' : false)
kkonganti@0:     mashscreen_w = false
kkonganti@0:     mashscreen_i = false
kkonganti@0:     mashscreen_v = false
kkonganti@0:     tuspy_run = true
kkonganti@0:     tuspy_s = false
kkonganti@0:     tuspy_m = false
kkonganti@0:     tuspy_ps = (params.cronology_root_dbdir
kkonganti@0:         + File.separator
kkonganti@0:         + 'index_pdg_metadata' 
kkonganti@0:         + File.separator
kkonganti@0:         + 'IDXD_PDG_METAD.pickle')
kkonganti@0:     tuspy_gd = (params.cronology_root_dbdir 
kkonganti@0:         + File.separator 
kkonganti@0:         + 'scaffold_genomes')
kkonganti@0:     tuspy_gds = '_scaffolded_genomic.fna.gz'
kkonganti@0:     tuspy_n = 5
kkonganti@0:     spades_run = false
kkonganti@0:     spades_isolate = false
kkonganti@0:     spades_sc = false
kkonganti@0:     spades_meta = true
kkonganti@0:     spades_bio = false
kkonganti@0:     spades_corona = false
kkonganti@0:     spades_rna = false
kkonganti@0:     spades_plasmid = false
kkonganti@0:     spades_metaviral = false
kkonganti@0:     spades_metaplasmid = false
kkonganti@0:     spades_rnaviral = false
kkonganti@0:     spades_iontorrent = false
kkonganti@0:     spades_only_assembler = false
kkonganti@0:     spades_careful = false
kkonganti@0:     spades_cov_cutoff = false
kkonganti@0:     spades_k = false
kkonganti@0:     spades_hmm = false
kkonganti@0:     shovill_run = true
kkonganti@0:     shovill_depth = 150
kkonganti@0:     shovill_gsize = false
kkonganti@0:     shovill_minlen = 200
kkonganti@0:     shovill_mincov = 5
kkonganti@0:     shovill_namefmt = false
kkonganti@0:     shovill_keepfiles = false
kkonganti@0:     shovill_assembler = false
kkonganti@0:     shovill_opts = false
kkonganti@0:     shovill_kmers = false
kkonganti@0:     shovill_trim = false
kkonganti@0:     shovill_noreadcorr = false
kkonganti@0:     shovill_nostitch = false
kkonganti@0:     shovill_nocorr = false
kkonganti@0:     polypolish_run = true
kkonganti@0:     polypolish_d = 5
kkonganti@0:     polypolish_i = 0.2
kkonganti@0:     polypolish_m = 10
kkonganti@0:     polypolish_v = 0.5
kkonganti@0:     mashtree_run = true
kkonganti@0:     mashtree_fofn = true
kkonganti@0:     mashtree_trunclength = false
kkonganti@0:     mashtree_mindepth = false
kkonganti@0:     mashtree_kmerlength = false
kkonganti@0:     mashtree_sketchsize = false 
kkonganti@0:     mashtree_seed = false
kkonganti@0:     mashtree_genomesize = false
kkonganti@0:     mashtree_sigfigs = false
kkonganti@0:     mashtree_sortorder = false
kkonganti@0:     quast_run = true
kkonganti@0:     quast_min_contig = 200
kkonganti@0:     quast_split_scaffolds = false
kkonganti@0:     quast_euk = false
kkonganti@0:     quast_fungal = false
kkonganti@0:     quast_large = false
kkonganti@0:     quast_k = false
kkonganti@0:     quast_kmer_size = false 
kkonganti@0:     quast_circos = true
kkonganti@0:     quast_glimmer = true
kkonganti@0:     quast_gene_thr = '0,300,1500,3000'
kkonganti@0:     quast_rna_finding = true
kkonganti@0:     quast_ref_size = false
kkonganti@0:     quast_ctg_thr = '0,1000,5000,10000,25000,50000'
kkonganti@0:     quast_x_for_nx = 90
kkonganti@0:     quast_use_all_alns = false
kkonganti@0:     quast_min_alignment = 65
kkonganti@0:     quast_min_identity = '95.0'
kkonganti@0:     quast_ambig_usage = 'one'
kkonganti@0:     quast_ambig_score = '0.99'
kkonganti@0:     quast_strict_na = false
kkonganti@0:     quast_x = 1000
kkonganti@0:     quast_local_mis_size = 200 
kkonganti@0:     quast_sca_gap_size = 10000
kkonganti@0:     quast_unaln_part_size = 500
kkonganti@0:     quast_skip_unaln_mis_ctgs = false
kkonganti@0:     quast_fragmented = false
kkonganti@0:     quast_frag_max_ident = (params.quast_fragmented ? 200 : false)
kkonganti@0:     quast_plots_format = 'pdf'
kkonganti@0:     prodka_run = true
kkonganti@0:     prodigal_run = (params.prodka_run ?: false)
kkonganti@0:     prodigal_c = false
kkonganti@0:     prodigal_f = 'gbk'
kkonganti@0:     prodigal_g = 11
kkonganti@0:     prodigal_m = true
kkonganti@0:     prodigal_n = false
kkonganti@0:     prodigal_p = 'single'
kkonganti@0:     prodigal_t = true
kkonganti@0:     prokka_run = (params.prodka_run ?: false)
kkonganti@0:     prokka_dbdir = false
kkonganti@0:     prokka_addgenes = true
kkonganti@0:     prokka_addmrna = true
kkonganti@0:     prokka_locustag = false
kkonganti@0:     prokka_increment = false
kkonganti@0:     prokka_gffver = false
kkonganti@0:     prokka_compliant = true
kkonganti@0:     prokka_centre = 'CFSAN-FDA'
kkonganti@0:     prokka_accver = false
kkonganti@0:     prokka_genus = (params.dpubmlstpy_org.capitalize() ?: 'Cronobacter')
kkonganti@0:     prokka_species = false
kkonganti@0:     prokka_strain = false
kkonganti@0:     prokka_plasmid = false
kkonganti@0:     prokka_kingdom = false
kkonganti@0:     prokka_gcode = false
kkonganti@0:     prokka_usegenus = true
kkonganti@0:     prokka_metagenome = false
kkonganti@0:     prokka_rawproduct = false
kkonganti@0:     prokka_cdsrnaolap = false
kkonganti@0:     prokka_evalue = '1e-9'
kkonganti@0:     prokka_coverage = 80
kkonganti@0:     prokka_fast = false
kkonganti@0:     prokka_noanno = false
kkonganti@0:     prokka_mincontiglen = 200 
kkonganti@0:     prokka_rfam = false
kkonganti@0:     prokka_norrna = false
kkonganti@0:     prokka_notrna = false
kkonganti@0:     prokka_rnammer = false
kkonganti@0:     pirate_run = false
kkonganti@0:     pirate_steps = false
kkonganti@0:     pirate_features = false 
kkonganti@0:     pirate_nucl = false
kkonganti@0:     pirate_pan_opt = "'--diamond'"
kkonganti@0:     pirate_pan_off = false
kkonganti@0:     pirate_min_len = 200
kkonganti@0:     pirate_para_off = false
kkonganti@0:     pirate_para_args = false
kkonganti@0:     pirate_classify_off = false
kkonganti@0:     pirate_align = false
kkonganti@0:     pirate_rplots = false
kkonganti@0:     mashscreen_run = true
kkonganti@0:     mash_sketch = (params.cronology_root_dbdir
kkonganti@0:         + File.separator 
kkonganti@0:         + 'mash_paste' 
kkonganti@0:         + File.separator 
kkonganti@0:         + 'msh.k21.1000h.' 
kkonganti@0:         + getPDGRelease(params.cronology_root_dbdir)
kkonganti@0:         + '.msh')
kkonganti@0:     mlst_run = true
kkonganti@0:     mlst_scheme = (params.dpubmlstpy_org ?: 'cronobacter')
kkonganti@0:     mlst_legacy = true
kkonganti@0:     mlst_minid = 95
kkonganti@0:     mlst_mincov = 10
kkonganti@0:     mlst_minscore = 50
kkonganti@0:     mlst_blastdb = false
kkonganti@0:     mlst_datadir = false
kkonganti@0:     mlst_label = true
kkonganti@0:     abricate_run = true
kkonganti@0:     abricate_datadir = ("${projectDir}"
kkonganti@0:         + File.separator
kkonganti@0:         + 'assets'
kkonganti@0:         + File.separator
kkonganti@0:         + 'abricate_dbs')
kkonganti@0:     abricate_minid = 90
kkonganti@0:     abricate_mincov = 80
kkonganti@0:     abricate_summary_run = true
kkonganti@0:     fq_filter_by_len = 0
kkonganti@0:     fq_suffix = (params.fq_single_end ? '.fastq.gz' : '_R1_001.fastq.gz')
kkonganti@0:     fq2_suffix = '_R2_001.fastq.gz'
kkonganti@0: }
kkonganti@0: 
kkonganti@0: /*
kkonganti@0: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
kkonganti@0:     METHOD TO GET PDG RELEASE AND SET REST OF THE METADATA FILES
kkonganti@0: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
kkonganti@0: */
kkonganti@0: 
kkonganti@0: def getPDGRelease(db_path) {
kkonganti@0: 	db_path_obj = new File( db_path )
kkonganti@0: 	def pdg_release = []
kkonganti@0: 
kkonganti@0: 	db_path_obj.eachDirRecurse {
kkonganti@0: 		it.toString() =~ /PDG/ \
kkonganti@0:             ? pdg_release << it.toString().find(/PDG\d+\.\d+/) \
kkonganti@0:             : null
kkonganti@0: 	}
kkonganti@0: 
kkonganti@0: 	return pdg_release.unique().join('').toString()
kkonganti@0: }