Merge pull request #25 from martinghunt/sspace_sd_bug

SSAPCE sd cannot exceed 1
sanger-pathogens · Apr 16, 2015 · f26abe5 · f26abe5
2 parents ef781a1 + 755db83
commit f26abe5
Showing 1 changed file with 1 addition and 0 deletions.
diff --git 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:')