Merge pull request #111 from jts/fix/metadata-na

Fix/metadata na
jts · Apr 25, 2024 · 04c7654 · 04c7654
2 parents 7a19778 + 3028cc6
commit 04c7654
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/workflow/scripts/format_pileup.py b/workflow/scripts/format_pileup.py
@@ -12,7 +12,7 @@
 
 print("\t".join(["contig", "position", "depth", "ref_base", "count_A", "count_C", "count_G", "count_T", "count_del", "alt_frequency"]))
 for pc in samfile.pileup(ignore_orphans=False):
-    freqs = {'A': 0, 'T': 0, 'G': 0, 'C': 0, '-': 0, 'R': 0 }
+    freqs = {'A': 0, 'T': 0, 'G': 0, 'C': 0, '-': 0, 'R': 0, 'N': 0 }
     for pr in pc.pileups:
         if pr.is_del:
             freqs['-'] += 1

diff --git a/workflow/scripts/plot/plot_qc_sequencing.R b/workflow/scripts/plot/plot_qc_sequencing.R
@@ -122,6 +122,9 @@ if(!interactive()) {
 
     if(args$metadata != "") {
         metadata.raw <- read.table(args$metadata, header=T, sep="\t")
+        if(all(is.na(metadata.raw$ct))) {
+          metadata.raw$ct <- 0
+        }
 
         # clean metadata of unknown Cts
         metadata = subset(metadata.raw, ct != "unknown")