polishing coverage bugfix #620

mikolmogorov · Sep 8, 2023 · 7198b7d · 7198b7d
1 parent 911957b
commit 7198b7d
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/flye/__build__.py b/flye/__build__.py
@@ -1 +1 @@
-__build__ = 1794
+__build__ = 1795
diff --git a/flye/polishing/polish.py b/flye/polishing/polish.py
@@ -125,6 +125,8 @@ def polish(contig_seqs, read_seqs, work_dir, num_iters, num_threads, read_platfo
     with open(stats_file, "w") as f:
         f.write("#seq_name\tlength\tcoverage\n")
         for ctg_id in contig_lengths:
+            if ctg_id not in coverage_stats:
+                coverage_stats[ctg_id] = 0
             f.write("{0}\t{1}\t{2}\n".format(ctg_id,
                     contig_lengths[ctg_id], coverage_stats[ctg_id]))