Identify Graph Threshold Areas

The Identify Graph Threshold Areas tool uses graph tracks as input to identify graph regions that fall within certain limits or thresholds.

To run the tool go to the toolbox:

        Toolbox | Track Tools (Image track_tools) | Graphs (Image graphs_closed_16_h_p) | Identify Graph Threshold Areas

Both a lower and an upper threshold can be specified to create an annotation track for those regions of a graph track where the values are in the given range (see figure 23.33). The range chosen for the lower and upper thresholds will depend on the data (coverage, quality etc).

Image GraphThresholdAreaWindow
Figure 23.33: Specification of lower and upper thresholds.

The window-size parameter specifies the width of the window around every position that is used to calculate an average value for that position and hence "smoothes" the graph track beforehand. A window size of 1 will simply use the value present at every individual position and determine if it is within the upper and lower threshold. In contrast, a window size of 100 checks if the average value derived from the surrounding 100 positions falls between the minimum and maximum threshold. Such larger windows help to prevent "jumps" in the graph track from fragmenting the output intervals or help to detect over-represented regions in the track that are only visible when looked at in the context of larger intervals and lower resolution.

It is also possible to restrict the tool to certain regions with specifying a region track.

An example output is shown in figure 23.34 where the coverage graph has some local minima. However, by using the averaging window, the tool is able to produce a single unbroken annotation covering the entire region. Of course larger window sizes result in regions that are broader and hence their boundaries are less likely to exactly coincide with the borders of visually recognizable borders of regions in the track.

Image GraphThresholdAreaResultTrackList
Figure 23.34: Track list including a read coverage graph, a reads track, and two read graph threshold graph generated to annotate regions where the coverage was above 100. The top track graph threshold was generated with a window size of 1, while the one from below was generated with a window size of 150.