Skip to content

Commit

Permalink
update flu output
Browse files Browse the repository at this point in the history
  • Loading branch information
fraser-combe committed Nov 4, 2024
1 parent 979d636 commit 25a2da2
Showing 1 changed file with 5 additions and 9 deletions.
14 changes: 5 additions & 9 deletions workflows/utilities/wf_flu_track.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -112,16 +112,12 @@ workflow flu_track {
docker = assembly_metrics_docker
}
}
# Calculate the percentage of mapped reads for flu samples
if (defined(irma.seg_ha_bam) || defined(irma.seg_na_bam)) {
call assembly_metrics.stats_n_coverage as flu_stats_n_coverage {
input:
samplename = samplename,
bamfile = select_first([irma.seg_ha_bam, irma.seg_na_bam])
}
}
# combine HA & NA assembly coverages
String ha_na_assembly_coverage_string = "HA: " + select_first([ha_assembly_coverage.depth, ""]) + ", NA: " + select_first([na_assembly_coverage.depth, ""])

# combine HA & NA mapped reads percentages
String ha_na_mapped_reads_percentage_string = "HA: " + select_first([ha_assembly_coverage.percentage_mapped_reads, ""]) + ", NA: " + select_first([na_assembly_coverage.percentage_mapped_reads, ""])

# ABRICATE will run if assembly is provided, or was generated with IRMA
if (defined(irma.irma_assemblies) && defined(irma.irma_assembly_fasta)){
call abricate.abricate_flu {
Expand Down Expand Up @@ -264,7 +260,7 @@ workflow flu_track {
File? irma_na_bam = irma.seg_na_bam
String ha_na_assembly_coverage = ha_na_assembly_coverage_string
# calulate mapped reads percentage for flu samples
Float? percentage_mapped_reads = flu_stats_n_coverage.percentage_mapped_reads
String percentage_mapped_reads = flu_stats_n_coverage.percentage_mapped_reads
# GenoFLU outputs
String? genoflu_version = genoflu.genoflu_version
String? genoflu_genotype = genoflu.genoflu_genotype
Expand Down

0 comments on commit 25a2da2

Please sign in to comment.