Commit 564d06e4 authored by cmaffeo2's avatar cmaffeo2
Browse files

Added logic to skip ssDNA in _get_intrahelical_neighbor, fixing ENM

parent 13f6cbb3
......@@ -40,6 +40,7 @@ class Nucleotide(Group):
def get_intrahelical_nucleotide(self,offset=1):
""" Returns nucleotide in same strand-side of helix with some offset along 5'-to-3' """
from ..segmentmodel import SingleStrandedSegment
strand_piece = self.parent
is_fwd = strand_piece.is_fwd
......@@ -60,7 +61,10 @@ class Nucleotide(Group):
return None
else:
seg, nt_idx, is_fwd = ret
return seg._get_atomic_nucleotide(nt_idx, is_fwd)
if isinstance(seg, SingleStrandedSegment):
return None
else:
return seg._get_atomic_nucleotide(nt_idx, is_fwd)
else:
return seg._get_atomic_nucleotide(nt_idx, is_fwd)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment