| 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 |
|
} |
