Skip to content

Commit

Permalink
apply the slice of _getSequence on sequence and not on data
Browse files Browse the repository at this point in the history
  • Loading branch information
ericloud committed Apr 10, 2017
1 parent 6078163 commit 3bb0518
Showing 1 changed file with 2 additions and 7 deletions.
9 changes: 2 additions & 7 deletions pyGeno/Chromosome.py
Original file line number Diff line number Diff line change
Expand Up @@ -64,25 +64,20 @@ def getSequenceData(self, slic) :
if sequenceModifier is not None :
if sequenceModifier.__class__ is SF.SequenceDel :
seqPos = seqPos + sequenceModifier.offset
#change the length of the sequence who can create some bug or side effect
#data = data[:seqPos] + data[seqPos + sequenceModifier.length:]
#To avoid to change the length of the sequence who can create some bug or side effect
data[seqPos:(seqPos + sequenceModifier.length)] = [''] * sequenceModifier.length
#raise TypeError("STOP!!!!!")

elif sequenceModifier.__class__ is SF.SequenceSNP :
data[seqPos] = sequenceModifier.alleles
elif sequenceModifier.__class__ is SF.SequenceInsert :
seqPos = seqPos + sequenceModifier.offset
data[seqPos] = "%s%s" % (data[seqPos], sequenceModifier.bases)
#raise TypeError("STOP!!!!!")
else :
raise TypeError("sequenceModifier on chromosome: %s starting at: %s is of unknown type: %s" % (self.chromosome.number, snp.start, sequenceModifier.__class__))

return data

def _getSequence(self, slic) :
# print data
return ''.join(self.getSequenceData(slic))
return ''.join(self.getSequenceData())[slic]

def __getitem__(self, i) :
return self._getSequence(i)
Expand Down

0 comments on commit 3bb0518

Please sign in to comment.