Skip to content

Commit

Permalink
Fixes CDS start and end for alternative reading frame
Browse files Browse the repository at this point in the history
The last fix didn't correctly handle the negative strand.
  • Loading branch information
courcelm authored Feb 20, 2018
1 parent 4cdd6da commit ae794eb
Showing 1 changed file with 5 additions and 1 deletion.
6 changes: 5 additions & 1 deletion pyGeno/Transcript.py
Original file line number Diff line number Diff line change
Expand Up @@ -117,7 117,11 @@ def setV(k, v) :

if len(cDNA) == 0 and e.frame != 0 :
e.CDS = e.CDS[e.frame:]
e.CDS_start = e.frame

if e.strand == ' ':
e.CDS_start = e.frame
else:
e.CDS_end -= e.frame

if len(e.CDS):
cDNA.append(''.join(e.CDS))
Expand Down

0 comments on commit ae794eb

Please sign in to comment.