fixed mean coverage calculation (#15)

jonas-fuchs · Dec 18, 2023 · 09bf627 · 09bf627
1 parent 6114f1c
commit 09bf627
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/bamdash/__init__.py b/bamdash/__init__.py
@@ -1,3 +1,3 @@
 """interactively visualize coverage and tracks"""
 _program = "bamdash"
-__version__ = "0.2.2"
+__version__ = "0.2.3"
diff --git a/bamdash/scripts/data.py b/bamdash/scripts/data.py
@@ -32,7 +32,7 @@ def get_coverage_stats(coverage_df, start, stop, min_cov):
         mean_coverage = 0
         recovery = 0
     else:
-        mean_coverage = statistics.mean(df_subset["coverage"])
+        mean_coverage = sum(df_subset["coverage"])/(stop-start)
         recovery = len(df_subset["coverage"])/(stop-start+1)*100
 
     return round(mean_coverage), round(recovery, 2)