diff --git a/tasks/taxon_id/task_nextclade.wdl b/tasks/taxon_id/task_nextclade.wdl index 5a54865a6..c435a6971 100644 --- a/tasks/taxon_id/task_nextclade.wdl +++ b/tasks/taxon_id/task_nextclade.wdl @@ -144,9 +144,19 @@ task nextclade_output_parser { nc_lineage = tsv_dict['lineage'] if nc_lineage is None: nc_lineage = "" + elif 'lineage' in tsv_dict: + nc_lineage = tsv_dict['Nextclade_pango'] + if nc_lineage is None: + nc_lineage = "" else: nc_lineage = "" Nextclade_Lineage.write(nc_lineage) + + with codecs.open("NEXTCLADE_QC", 'wt') as Nextclade_QC: + nc_qc = tsv_dict['qc.overallStatus'] + if nc_qc == '': + nc_qc = 'NA' + Nextclade_QC.write(nc_qc) CODE >>> runtime { @@ -164,6 +174,7 @@ task nextclade_output_parser { String nextclade_tamiflu_aa_subs = read_string("TAMIFLU_AASUBS") String nextclade_aa_dels = read_string("NEXTCLADE_AADELS") String nextclade_lineage = read_string("NEXTCLADE_LINEAGE") + String nextclade_qc = read_string("NEXTCLADE_QC") } } diff --git a/workflows/theiacov/wf_theiacov_illumina_pe.wdl b/workflows/theiacov/wf_theiacov_illumina_pe.wdl index 8dff1ee42..11fba29c9 100644 --- a/workflows/theiacov/wf_theiacov_illumina_pe.wdl +++ b/workflows/theiacov/wf_theiacov_illumina_pe.wdl @@ -352,6 +352,7 @@ workflow theiacov_illumina_pe { String nextclade_aa_dels = select_first([ha_na_nextclade_aa_dels, nextclade_output_parser.nextclade_aa_dels, ""]) String nextclade_clade = select_first([nextclade_output_parser.nextclade_clade, ""]) String? nextclade_lineage = nextclade_output_parser.nextclade_lineage + String? nextclade_qc = nextclade_output_parser.nextclade_qc # Nextclade Flu outputs - NA specific columns - tamiflu mutation String? nextclade_tamiflu_resistance_aa_subs = nextclade_output_parser_flu_na.nextclade_tamiflu_aa_subs # VADR Annotation QC