From dcce5692a9fe6c6fc76cadebbf960c4ea4f81f21 Mon Sep 17 00:00:00 2001 From: HNHalstead <47102761+HNHalstead@users.noreply.github.com> Date: Fri, 2 Feb 2024 11:43:21 -0800 Subject: [PATCH] Update task_qc_check_phb.wdl --- tasks/quality_control/task_qc_check_phb.wdl | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/tasks/quality_control/task_qc_check_phb.wdl b/tasks/quality_control/task_qc_check_phb.wdl index 66fbb3df9..717deb28f 100644 --- a/tasks/quality_control/task_qc_check_phb.wdl +++ b/tasks/quality_control/task_qc_check_phb.wdl @@ -476,7 +476,8 @@ task qc_check_phb_waphl { import pandas as pd import numpy as np - + + alert = "" qc_check = "no qc check for this taxa" if "~{expected_taxon}" and "~{predicted_taxon}": if "~{expected_taxon}" != "~{predicted_taxon}": @@ -505,7 +506,7 @@ task qc_check_phb_waphl { if "~{est_coverage_raw}" and qc_check == "PASS": if float("~{est_coverage_raw}") < float(qc_check_df.loc["~{predicted_taxon}", "est_coverage_raw"]): qc_check = "FAIL" - + if "~{number_N}" and "~{number_Total}" and qc_check == "PASS": perc_N = int("~{number_N}")/int("~{number_Total}") if perc_N >= float(qc_check_df.loc["~{predicted_taxon}", "perc_N"]):