dont add 0 reads to summary

artic-network · Jan 18, 2024 · 73f0f77 · 73f0f77
1 parent 17925e9
commit 73f0f77
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/bin/extract_taxa_from_reads.py b/bin/extract_taxa_from_reads.py
@@ -381,6 +381,9 @@ def extract_taxa(
     sys.stderr.write("Write summary\n")
     summary = []
     for taxon in lists_to_extract:
+        if out_counts[taxon] == 0:
+            sys.stderr.write("No reads extracted  for taxid %s\n" %taxon)
+            continue
         if reads2:
             summary.append(
                 {