kkonganti@0: // Help text for quast within CPIPES.
kkonganti@0: 
kkonganti@0: def quastHelp(params) {
kkonganti@0: 
kkonganti@0:     Map tool = [:]
kkonganti@0:     Map toolspecs = [:]
kkonganti@0:     tool.text = [:]
kkonganti@0:     tool.helpparams = [:]
kkonganti@0: 
kkonganti@0:     toolspecs = [
kkonganti@0:         'quast_run': [
kkonganti@0:             clihelp: 'Run quast tool. Default: ' +
kkonganti@0:                 (params.quast_run ?: false),
kkonganti@0:             cliflag: null,
kkonganti@0:             clivalue: null
kkonganti@0:         ],
kkonganti@0:         'quast_min_contig': [
kkonganti@0:             clihelp: 'Lower threshold for contig length. ' +
kkonganti@0:                 "Default: ${params.quast_min_contig}",
kkonganti@0:             cliflag: '-m',
kkonganti@0:             clivalue: (params.quast_min_contig ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_split_scaffolds': [
kkonganti@0:             clihelp: "Split assemblies by continuous fragments of N's" +
kkonganti@0:                 'and add such "contigs" to the comparison' +
kkonganti@0:                 "Default: ${params.quast_split_scaffolds}",
kkonganti@0:             cliflag: '-s',
kkonganti@0:             clivalue: (params.quast_split_scaffolds ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_euk': [
kkonganti@0:             clihelp: 'Genome is eukaryotic (primarily affects gene prediction). ' +
kkonganti@0:                 "Default: ${params.quast_euk}",
kkonganti@0:             cliflag: '-e',
kkonganti@0:             clivalue: (params.quast_euk ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_fungal': [
kkonganti@0:             clihelp: 'Genome is fungal (primarily affects gene prediction). ' +
kkonganti@0:                 "Default: ${params.quast_fungal}",
kkonganti@0:             cliflag: '--fungal',
kkonganti@0:             clivalue: (params.quast_fungal ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_large': [
kkonganti@0:             clihelp: 'Use optimal parameters for evaluation of large genomes' +
kkonganti@0:                 "In particular, imposes '-e -m 3000 -i 500 -x 7000'." +
kkonganti@0:                 "Default: ${params.quast_large}",
kkonganti@0:             cliflag: '--large',
kkonganti@0:             clivalue: (params.quast_large ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_k': [
kkonganti@0:             clihelp: 'Compute k-mer-based quality metrics (recommended for large genomes). ' +
kkonganti@0:                 "Default: ${params.quast_k}",
kkonganti@0:             cliflag: '-k',
kkonganti@0:             clivalue: (params.quast_k ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_kmer_size': [
kkonganti@0:             clihelp: 'Size of k used in --quast_k. ' +
kkonganti@0:                 "Default: ${params.quast_kmer_size}",
kkonganti@0:             cliflag: '--k-mer-size',
kkonganti@0:             clivalue: (params.quast_kmer_size ?: 101)
kkonganti@0:         ],
kkonganti@0:         'quast_circos': [
kkonganti@0:             clihelp: 'Draw circos plot. ' +
kkonganti@0:                 "Default: ${params.quast_circos}",
kkonganti@0:             cliflag: '--circos',
kkonganti@0:             clivalue: (params.quast_circos ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_glimmer': [
kkonganti@0:             clihelp: 'Use GlimmerHMM for gene prediction. ' +
kkonganti@0:                 "Default: ${params.quast_glimmer}",
kkonganti@0:             cliflag: '--glimmer',
kkonganti@0:             clivalue: (params.quast_glimmer ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_gene_thr': [
kkonganti@0:             clihelp: 'Comma-separated list of threshold lengths of genes to ' +
kkonganti@0:                 'search with Gene Finding module. ' +
kkonganti@0:                 "Default: ${params.quast_gene_thr}",
kkonganti@0:             cliflag: '--gene-thresholds',
kkonganti@0:             clivalue: (params.quast_gene_thr ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_rna_finding': [
kkonganti@0:             clihelp: 'Predict ribosomal RNA genes using Barrnap. ' +
kkonganti@0:                 "Default: ${params.quast_rna_finding}",
kkonganti@0:             cliflag: '--rna-finding',
kkonganti@0:             clivalue: (params.quast_rna_finding ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_ref_size': [
kkonganti@0:             clihelp: 'Estimated reference size in-case the reference genome is not supplied. ' +
kkonganti@0:                 "Default: ${params.quast_ref_size}",
kkonganti@0:             cliflag: '--est-ref-size',
kkonganti@0:             clivalue: (params.quast_ref_size ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_ctg_thr': [
kkonganti@0:             clihelp: 'Comma-separated list of contig length thresholds. ' +
kkonganti@0:                 "Default: ${params.quast_ctg_thr}",
kkonganti@0:             cliflag: '--contig-thresholds',
kkonganti@0:             clivalue: (params.quast_ctg_thr ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_x_for_nx': [
kkonganti@0:             clihelp: "Value of 'x' for Nx, Lx, etc metrics reported in addition to N50, L50, etc. " +
kkonganti@0:                 "Default: ${params.quast_x_for_nx}",
kkonganti@0:             cliflag: '--x-for-Nx',
kkonganti@0:             clivalue: (params.quast_x_for_nx ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_use_all_alns': [
kkonganti@0:             clihelp: 'Compute genome fraction, # genes, # operons in QUAST v1.* style. ' +
kkonganti@0:                 "By default, QUAST filters Minimap's alignments to keep only best ones. " +
kkonganti@0:                 "Default: ${params.quast_glimmer}",
kkonganti@0:             cliflag: '--use-all-alignments',
kkonganti@0:             clivalue: (params.quast_use_all_alns ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_min_alignment': [
kkonganti@0:             clihelp: 'The minimum alignment length. ' +
kkonganti@0:                 "Default: ${params.quast_min_alignment}",
kkonganti@0:             cliflag: '-i',
kkonganti@0:             clivalue: (params.quast_min_alignment ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_min_identity': [
kkonganti@0:             clihelp: 'The minimum alignment identity (80.0, 100.0). ' +
kkonganti@0:                 "Default: ${params.quast_min_identity}",
kkonganti@0:             cliflag: '--min-identity',
kkonganti@0:             clivalue: (params.quast_min_identity ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_ambig_usage': [
kkonganti@0:             clihelp: 'Use none, one, or all alignments of a contig when all of them. ' +
kkonganti@0:                 'are almost equally good (see --quast_ambig_score). ' +
kkonganti@0:                 "Default: ${params.quast_ambig_usage}",
kkonganti@0:             cliflag: '-a',
kkonganti@0:             clivalue: (params.quast_ambig_usage ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_ambig_score': [
kkonganti@0:             clihelp: 'Score S for defining equally good alignments of a single contig. ' +
kkonganti@0:                 'All alignments are sorted by decreasing LEN * IDY% value. ' +
kkonganti@0:                 'All alignments with LEN * IDY% < S * best(LEN * IDY%) are ' +
kkonganti@0:                 'discarded. S should be between 0.8 and 1.0. ' +
kkonganti@0:                 "Default: ${params.quast_ambig_score}",
kkonganti@0:             cliflag: '--ambiguity-score',
kkonganti@0:             clivalue: (params.quast_ambig_score ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_strict_na': [
kkonganti@0:             clihelp: 'Break contigs in any misassembly event when compute NAx and NGAx. '+
kkonganti@0:                 'By default, QUAST breaks contigs only by extensive misassemblies (not local ones). ' +
kkonganti@0:                 "Default: ${params.quast_strict_na}",
kkonganti@0:             cliflag: '--strict-NA',
kkonganti@0:             clivalue: (params.quast_strict_na ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_x': [
kkonganti@0:             clihelp: 'Lower threshold for extensive misassembly size. All relocations with inconsistency ' +
kkonganti@0:                 'less than extensive-mis-size are counted as local misassemblies. ' +
kkonganti@0:                 "Default: ${params.quast_x}",
kkonganti@0:             cliflag: '-x',
kkonganti@0:             clivalue: (params.quast_x ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_local_mis_size': [
kkonganti@0:             clihelp: 'Lower threshold on local misassembly size. Local misassemblies with inconsistency ' +
kkonganti@0:                 'less than local-mis-size are counted as (long) indels. ' +
kkonganti@0:                 "Default: ${params.quast_local_mis_size}",
kkonganti@0:             cliflag: '--local-mis-size',
kkonganti@0:             clivalue: (params.quast_local_mis_size ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_sca_gap_size': [
kkonganti@0:             clihelp: 'Max allowed scaffold gap length difference. All relocations with inconsistency ' +
kkonganti@0:                 'less than scaffold-gap-size are counted as scaffold gap misassemblies. ' +
kkonganti@0:                 "Default: ${params.quast_sca_gap_size}",
kkonganti@0:             cliflag: '--scaffold-gap-max-size',
kkonganti@0:             clivalue: (params.quast_sca_gap_size ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_unaln_part_size': [
kkonganti@0:             clihelp: 'Lower threshold for detecting partially unaligned contigs. Such contig should have ' +
kkonganti@0:                 'at least one unaligned fragment >= the threshold. ' +
kkonganti@0:                 "Default: ${params.quast_unaln_part_size}",
kkonganti@0:             cliflag: '--unaligned-part-size',
kkonganti@0:             clivalue: (params.quast_unaln_part_size ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_skip_unaln_mis_ctgs': [
kkonganti@0:             clihelp: 'Do not distinguish contigs with >= 50% unaligned bases as a separate group ' +
kkonganti@0:                 'By default, QUAST does not count misassemblies in them. ' +
kkonganti@0:                 "Default: ${params.quast_skip_unaln_mis_ctgs}",
kkonganti@0:             cliflag: '--skip-unaligned-mis-contigs',
kkonganti@0:             clivalue: (params.quast_skip_unaln_mis_ctgs ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_fragmented': [
kkonganti@0:             clihelp: 'Reference genome may be fragmented into small pieces (e.g. scaffolded reference). ' +
kkonganti@0:                 "Default: ${params.quast_fragmented}",
kkonganti@0:             cliflag: '--fragmented',
kkonganti@0:             clivalue: (params.quast_fragmented ? ' ' : '')
kkonganti@0:         ],
kkonganti@0:         'quast_frag_max_ident': [
kkonganti@0:             clihelp: 'Mark translocation as fake if both alignments are located no further than N bases ' +
kkonganti@0:                 'from the ends of the reference fragments. ' +
kkonganti@0:                 "Default: ${params.quast_frag_max_ident}",
kkonganti@0:             cliflag: '--fragmented-max-indent',
kkonganti@0:             clivalue: (params.quast_frag_max_ident ?: '')
kkonganti@0:         ],
kkonganti@0:         'quast_plots_format': [
kkonganti@0:             clihelp: 'Save plots in specified format [default: pdf]. ' +
kkonganti@0:                 'Supported formats: emf, eps, pdf, png, ps, raw, rgba, svg, svgz' +
kkonganti@0:                 "Default: ${params.quast_plots_format}",
kkonganti@0:             cliflag: '--plots-format',
kkonganti@0:             clivalue: (params.quast_plots_format ?: '')
kkonganti@0:         ],
kkonganti@0:     ]
kkonganti@0: 
kkonganti@0:     toolspecs.each {
kkonganti@0:         k, v -> tool.text['--' + k] = "${v.clihelp}"
kkonganti@0:         tool.helpparams[k] = [ cliflag: "${v.cliflag}", clivalue: v.clivalue ]
kkonganti@0:     }
kkonganti@0: 
kkonganti@0:     return tool
kkonganti@0: }