To run the Whole Genome Coverage Analysis tool:
Toolbox | Resequencing Analysis () | Whole Genome Coverage Analysis ()
In the first dialog, select a reads track or read mapping and click Next. This opens the dialog shown in figure 27.20.
Set the p-value and minimum length cutoff. Click Next and specify the result handling (figure 27.21).
Selecting the "Create report" will generate the report made of 2 tables (figure 27.22). The first one, called References, lists per chromosome the number of reads, their length, and how many signatures of unexpectedly low or high coverage was found in the mapping. The second table lists on 2 rows low and high coverage signatures found, as well as how many reads were used to calculate these signatures.
Selecting the "Create regions" will generate the annotation track carrying the name of the original file followed by (COV). This file can be visualized as an annotation track or as a table depending on the users choice. The annotation table contains a row for each detected low or high coverage region, with information describing the location, the type and the p-value of the detected region. The p-value of a region is defined as the average of the p-values calculated for each of the positions in the region.
An example of a track output of the Whole Genome Coverage Analysis tool is shown in figure 27.23.
The Whole Genome Coverage Analysis table includes the following columns (figure 27.23):
- The name is taken from the reference sequence used for mapping
- The start and end position of this region on the reference sequence
- The type of annotation: high or low coverage
- The calculated significance p-value for this region
- The length of the region
For the visual inspection and comparison to known gene/transcripts or other kind of annotations, all region are also annotated on the read mapping.