Skip to content

Commit d5d3593

Browse files
Dom LaetschDom Laetsch
authored andcommitted
Fixed cas cov parsing
1 parent 40d370b commit d5d3593

File tree

2 files changed

+3
-2
lines changed

2 files changed

+3
-2
lines changed

lib/BtCore.py

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -300,7 +300,7 @@ def parseCovs(self, covLibObjs):
300300
for name, cov in cov_dict.items():
301301
covLib.cov_sum += cov
302302
self.dict_of_blobs[name].addCov(covLib.name, cov)
303-
self.dict_of_blobs[name].read_cov = {covLib.name : read_cov_dict[name]}
303+
self.dict_of_blobs[name].addReadCov(covLib.name, read_cov_dict[name])
304304

305305
elif covLib.fmt == 'cov':
306306
cov_dict = BtIO.readCov(covLib.f, set(self.dict_of_blobs))

lib/BtPlot.py

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -324,7 +324,8 @@ def compute_stats(self):
324324
for cov_lib in self.cov_libs:
325325
stats[label]['cov_mean'][cov_lib] = mean(array(stats[label]['covs'][cov_lib])) if stats[label]['count_visible'] > 0 else 0.0
326326
stats[label]['cov_std'][cov_lib] = std(array(stats[label]['covs'][cov_lib])) if stats[label]['count_visible'] > 0 else 0.0
327-
stats[label]['reads_mapped_perc'][cov_lib] = stats[label]['reads_mapped'][cov_lib]/self.cov_libs_total_reads_dict[cov_lib]
327+
if self.cov_libs_total_reads_dict[cov_lib]:
328+
stats[label]['reads_mapped_perc'][cov_lib] = stats[label]['reads_mapped'][cov_lib]/self.cov_libs_total_reads_dict[cov_lib]
328329

329330
self.stats = stats
330331

0 commit comments

Comments
 (0)