Detection and characterization of low and high genome coverage regions using an efficient running median and a double threshold approach
Abstract
Motivation
Next Generation Sequencing (NGS) provides researchers with powerful tools to investigate both prokaryotic and eukaryotic genetics. An accurate assessment of reads mapped to a specific genome consists of inspecting thegenome coverageas number of reads mapped to a specific genome location. Most current methods use the average of the genome coverage (sequencing depth) to summarize the overall coverage. This metric quickly assess the sequencing quality but ignores valuable biological information like the presence of repetitive regions or deleted genes. The detection of such information may be challenging due to a wide spectrum of heterogeneous coverage regions, a mixture of underlying models or the presence of a non-constant trend along the genome. Using robust statistics to systematically identify genomic regions with unusual coverage is needed to characterize these regions more precisely.
Results
We implemented an efficient running median algorithm to estimate the genome coverage trend. The distribution of the normalized genome coverage is then estimated using a Gaussian mixture model. A<monospace>z</monospace>-score statistics is then assigned to each base position and used to separate the central distribution from the regions of interest (ROI) (i.e., under and over-covered regions). Finally, a double threshold mechanism is used to cluster the genomic ROIs. HTML reports provide a summary with interactive visual representations of the genomic ROIs.
Availability
An implementation of the genome coverage characterization is available within the<monospace>Sequana</monospace>project. The standalone application is called<monospace>sequana_coverage</monospace>. The source code is available on<monospace>GitHub</monospace>(<ext-link xmlns:xlink="http://www.w3.org/1999/xlink" ext-link-type="uri" xlink:href="http://github.com/sequana/sequana">http://github.com/sequana/sequana</ext-link>), and documentation on<monospace>ReadTheDocs</monospace>(<ext-link xmlns:xlink="http://www.w3.org/1999/xlink" ext-link-type="uri" xlink:href="http://sequana.readtheodcs.org">http://sequana.readtheodcs.org</ext-link>). An example of HTML report is provided on<ext-link xmlns:xlink="http://www.w3.org/1999/xlink" ext-link-type="uri" xlink:href="http://sequana.github.io">http://sequana.github.io</ext-link>.
Contact
<email>dimitri.desvillechabrol@pasteur.fr</email>,<email>thomas.cokelaer@pasteur.fr</email>
Related articles
Related articles are currently not available for this article.