163 |
|
fusion = true; |
164 |
|
else |
165 |
|
{ |
166 |
< |
int inter_dist = 0; |
166 |
> |
int inner_dist = 0; |
167 |
|
if (lh.antisense_align()) |
168 |
< |
inter_dist = lh.left() - rh.right(); |
168 |
> |
// used rh.left() instead of rh.right() for the cases, |
169 |
> |
// where reads overlap with each other or reads span introns |
170 |
> |
inner_dist = lh.left() - rh.left(); |
171 |
|
else |
172 |
< |
inter_dist = rh.left() - lh.right(); |
172 |
> |
inner_dist = rh.left() - lh.left(); |
173 |
|
|
174 |
< |
if (inter_dist < -(int)max_insertion_length || inter_dist > (int)fusion_min_dist) |
174 |
> |
if (inner_dist < 0 || inner_dist > (int)fusion_min_dist) |
175 |
|
fusion = true; |
176 |
|
} |
177 |
|
} |