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 )
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.)
