From 7198b7da2a5a244352204bb5c6bc68a4d50e08ce Mon Sep 17 00:00:00 2001 From: Mikhail Kolmogorov Date: Fri, 8 Sep 2023 12:21:59 -0700 Subject: [PATCH] polishing coverage bugfix #620 --- flye/__build__.py | 2 +- flye/polishing/polish.py | 2 ++ 2 files changed, 3 insertions(+), 1 deletion(-) diff --git a/flye/__build__.py b/flye/__build__.py index a3cc96e00..b839c4f53 100644 --- 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 index fd78ad90e..78060c170 100644 --- 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]))