diff --git a/pipes/WDL/workflows/assemble_refbased.wdl b/pipes/WDL/workflows/assemble_refbased.wdl index 3001f9dfb..f624252e4 100644 --- a/pipes/WDL/workflows/assemble_refbased.wdl +++ b/pipes/WDL/workflows/assemble_refbased.wdl @@ -67,18 +67,18 @@ workflow assemble_refbased { Float major_cutoff=0.75 Boolean skip_mark_dupes=false File? trim_coords_bed - } - Map[String,String] align_to_ref_options = { + Map[String,String] align_to_ref_options = { "novoalign": "-r Random -l 40 -g 40 -x 20 -t 501 -k", "bwa": "-k 12 -B 1", "minimap2": "" } - Map[String,String] align_to_self_options = { + Map[String,String] align_to_self_options = { "novoalign": "-r Random -l 40 -g 40 -x 20 -t 100", "bwa": "", "minimap2": "" } + } scatter(reads_unmapped_bam in reads_unmapped_bams) { call assembly.align_reads as align_to_ref { diff --git a/pipes/WDL/workflows/sarscov2_illumina_full.wdl b/pipes/WDL/workflows/sarscov2_illumina_full.wdl index b687d38c4..0c13d0d18 100644 --- a/pipes/WDL/workflows/sarscov2_illumina_full.wdl +++ b/pipes/WDL/workflows/sarscov2_illumina_full.wdl @@ -51,6 +51,7 @@ workflow sarscov2_illumina_full { Int min_genome_bases = 24000 Int max_vadr_alerts = 0 Int ntc_max_unambig = 3000 + Int? min_genome_coverage File? sample_rename_map @@ -121,11 +122,10 @@ workflow sarscov2_illumina_full { reads_unmapped_bams = name_reads.right, reference_fasta = reference_fasta, sample_name = name_reads.left, - aligner = "minimap2", skip_mark_dupes = ampseq, trim_coords_bed = bed_rename.outfile, major_cutoff = 0.75, - min_coverage = if ampseq then 50 else 3 + min_coverage = if defined(min_genome_coverage) then min_genome_coverage else (if ampseq then 50 else 3) } # log controls