From 755db83f328fd7f5f766f38dc51c754f063fe205 Mon Sep 17 00:00:00 2001 From: Martin Hunt Date: Thu, 16 Apr 2015 08:13:02 +0100 Subject: [PATCH] SSAPCE sd cannot exceed 1 --- ariba/clusters.py | 1 + 1 file changed, 1 insertion(+) diff --git a/ariba/clusters.py b/ariba/clusters.py index 00912fb7..e6784bf1 100644 --- a/ariba/clusters.py +++ b/ariba/clusters.py @@ -246,6 +246,7 @@ def _bam_to_clusters_reads(self): def _set_insert_size_data(self): assert len(self.insert_hist) > 0 (x, self.insert_size, pc95, self.insert_sspace_sd) = self.insert_hist.stats() + self.insert_sspace_sd = min(1, self.insert_sspace_sd) self.insert_proper_pair_max = 1.1 * pc95 if self.verbose: print('\nInsert size information from reads mapped to reference genes:')