From 4027b3ae2ee07011dc72a7b98db1ee28b8c1e9c1 Mon Sep 17 00:00:00 2001
From: Chris Maffeo <cmaffeo2@illinois.edu>
Date: Wed, 19 Sep 2018 12:40:02 -0500
Subject: [PATCH] clear_atomic now clears sequences from segments when they are
the wrong length
---
mrdna/segmentmodel.py | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/mrdna/segmentmodel.py b/mrdna/segmentmodel.py
index 7f3c015..8f1b305 100644
--- a/mrdna/segmentmodel.py
+++ b/mrdna/segmentmodel.py
@@ -1573,6 +1573,10 @@ class SegmentModel(ArbdModel):
for d in ('fwd','rev'):
seg.strand_pieces[d] = []
self._generate_strands()
+ ## Clear sequence if needed
+ for seg in self.segments:
+ if seg.sequence is not None and len(seg.sequence) != seg.num_nt:
+ seg.sequence = None
def clear_beads(self):
return self._clear_beads()
@@ -2428,7 +2432,7 @@ class SegmentModel(ArbdModel):
for d in ('fwd','rev'):
seg.strand_pieces[d] = []
except:
- ...
+ pass
self.strands = strands = []
""" Ensure unconnected ends have 5prime Location objects """
--
GitLab