Opened 10 years ago

Last modified 10 years ago

#658 closed enhancement

Adjust fragment size mean and standard deviation — at Version 1

Reported by: Nicklas Nordborg Owned by: Nicklas Nordborg
Priority: critical Milestone: Reggie v2.16.1
Component: net.sf.basedb.reggie Keywords:
Cc:

Description (last modified by Jari Häkkinen)

In the demux step we introduced a call to bowtie that is used to calculate the average fragment size and standard deviation (see [2680]). The values are used as input to Tophat for -mate-inner-dist and --mate-std-dev (see [2684]).

But after alignment Tophat reports values which are typically 13 (inner-dist) and 10 (std-dev) higher. Adjusting the input parameters with this amount gives higher alignment rate.

We tested effect of the ad hoc addition to the bowtie numbers on the overall tophat alignment, and compared the tophat aligned fragment lengths and widths with the input values to tophat. The net result is that the alignment rate increases approximately 3 per cent over a sample set of 150, and the lengths and widths are retained within a standard deviation of 0.02. (Data not shown.)

Change History (1)

comment:1 by Jari Häkkinen, 10 years ago

Description: modified (diff)
Note: See TracTickets for help on using tickets.