add safety check for scale == 0

rnajena · Nov 27, 2024 · 249176d · 249176d
1 parent 314b132
commit 249176d
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/read5/AbstractFileReader.py b/read5/AbstractFileReader.py
@@ -1005,6 +1005,10 @@ def getPolyAStandardizedSignal(self, readid : str, polyAstart : int, polyAend :
             shift = np.mean(polyASignal)
             scale = np.std(polyASignal)
 
+        # safety check when working with simulated signals where stdev = 0
+        if scale == 0:
+            scale = 1
+
         return ((self.getpASignal(readid) - shift) / scale) * polyAStdev + polyAMean
 
     def getStartTimeInMinutes(self, readid) -> float: