Create Mapping Graph
Create Mapping Graph can create a range of different graphs from a read mapping track. To run the tool, go to:
Toolbox | Utility Tools (
) | Tracks (
) | Graphs (
) | Create Mapping Graph (
)
Select a read mapping as input (figure 27.33).
Figure 27.33: Mapping graph tracks containing different types of information can be created using Create Mapping Graph.
Select the graph track types to create (figure 27.34). One graph track is created for each type selected.
Figure 27.34: Choose the types of graph tracks to be created.
Each position in a graph track contains the values of the track type selected. These are:
- Read coverage The number of reads contributing to the alignment at the position. A more detailed definition is in Reference sequence statistics.
- Non-specific read coverage The number of reads mapped at the position that would map equally well to other places in the reference sequence.
- Specific read coverage The number of reads that map uniquely at the position, i.e. they do not map equally well to other places in the reference sequence.
- Unaligned ends coverage The number of reads with unaligned ends at the position. Unaligned ends arise when a read has been locally aligned and the there are mismatches or gaps relative to the reference sequence at the end of the read. Unaligned regions do not contribute to coverage in other graph track types.
- Non-perfect read coverage The number of reads at the position with one or more mismatches or gaps relative to the reference sequence.
- Paired read coverage The number of intact read pairs mapped to the position. Coverage is counted as one in positions where the reads of a pair overlap.
- Paired read specific coverage The number of intact paired reads that map uniquely at the position, i.e. they do not map equally well to other places in the reference sequence.
- Paired end distance The average distance between the forward and reverse reads of pairs mapped to the position.
- Broken pair coverage The number of broken paired reads mapped to the position. A pair is marked as broken if only one read in the pair matches the reference, if the reads map to different chromosomes, or if the distance or relative orientation between the reads is not within the expected values.
- Reads start coverage The number of reads with their start mapped to the position.
- Forward read coverage The number of reads mapping in forward direction. First and second read of a pair will be counted separately.
- Reverse read coverage The number of reads mapping in reverse direction. First and second read of a pair will be counted separately.
An example of three different outputs is shown in figure 27.35.
Figure 27.35: Three types of graph tracks
The option "Fix graph bounds" found under Track layout in the Side Panel is useful to manually adjust the numbers on the y-axis.
When zoomed out, the graph tracks are composed of three curves showing the maximum, mean, and minimum value observed in a given region (figure 27.38). When zoomed in to single base resolution only one curve is shown, reflecting the exact observation at each individual position (figure 27.36).
Figure 27.36: Zoomed in on a graph track