@@ -142,13 +142,19 @@ def curate_alt_prime(event_list, options):
142142 elif event_list [i ].exons1 [1 , 1 ] < event_list [i ].exons2 [1 , 1 ]:
143143 event_list [i ].exons2 [1 , 1 ] = event_list [i ].exons1 [1 , 1 ]
144144 corr_count += 1
145+ ### check whether our isoform convention is still met - if not, correct
146+ if np .sum (event_list [i ].exons1 [:, 1 ] - event_list [i ].exons1 [:, 0 ]) > np .sum (event_list [i ].exons2 [:, 1 ] - event_list [i ].exons2 [:, 0 ]):
147+ event_list [i ].exons1 , event_list [i ].exons2 = event_list [i ].exons2 , event_list [i ].exons1
145148 elif np .all (event_list [i ].exons1 [1 , :] == event_list [i ].exons2 [1 , :]):
146149 if event_list [i ].exons1 [0 , 0 ] > event_list [i ].exons2 [0 , 0 ]:
147150 event_list [i ].exons2 [0 , 0 ] = event_list [i ].exons1 [0 , 0 ]
148151 corr_count += 1
149152 elif event_list [i ].exons1 [0 , 0 ] < event_list [i ].exons2 [0 , 0 ]:
150153 event_list [i ].exons1 [0 , 0 ] = event_list [i ].exons2 [0 , 0 ]
151154 corr_count += 1
155+ ### check whether our isoform convention is still met - if not, correct
156+ if np .sum (event_list [i ].exons1 [:, 1 ] - event_list [i ].exons1 [:, 0 ]) > np .sum (event_list [i ].exons2 [:, 1 ] - event_list [i ].exons2 [:, 0 ]):
157+ event_list [i ].exons1 , event_list [i ].exons2 = event_list [i ].exons2 , event_list [i ].exons1
152158
153159 ### remove events with non-overlapping alt_exons
154160 if len (rm_idx ) > 0 :
0 commit comments