Opened 3 years ago

Last modified 3 years ago

#1337 closed defect

Bug in Tophat when --mate-inner-dist + --mate-std-dev < 0 — at Version 1

Reported by: Nicklas Nordborg Owned by: Nicklas Nordborg
Priority: major Milestone: Reggie v4.33
Component: net.sf.basedb.reggie Keywords:
Cc:

Description (last modified by Nicklas Nordborg)

This was initially thought to be a problem with an incorrect value for the --mate-inner-dist parameter due to lot of clipping when sequencing 2x150 reads (see #1336 and #1335).

However, investigations shows that there is a bug in Tophat when marking aligned pairs as a proper pairs that is triggered when the --mate-inner-dist parameter is negative and --mate-std-dev is less than abs(--mate-inner-dist). For example:

--mate-inner-dist -80 --mate-std-dev 79 // will trigger the bug
--mate-inner-dist -79 --mate-std-dev 79 // works

To avoid triggering the bug we should make sure that the --mate-inner-dist parameter always is large enough.

Change History (2)

comment:1 by Nicklas Nordborg, 3 years ago

Description: modified (diff)

by Nicklas Nordborg, 3 years ago

Note: See TracTickets for help on using tickets.