A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data

Abstract

Motivation: Most existing methods for DNA sequence analysis rely on accurate sequences or genotypes. However, in applications of the next-generation sequencing (NGS), accurate genotypes may not be easily obtained (e.g. multi-sample low-coverage sequencing or somatic mutation discovery). These applications press for the development of new methods for analyzing sequence data with uncertainty. Results: We present a statistical framework for calling SNPs, discovering somatic mutations, inferring population genetical parameters and performing association tests directly based on sequencing data without explicit genotyping or linkage-based imputation. On real data, we demonstrate that our method achieves comparable accuracy to alternative methods for estimating site allele count, for inferring allele frequency spectrum and for association mapping. We also highlight the necessity of using symmetric datasets for finding somatic mutations and confirm that for discovering rare events, mismapping is frequently the leading source of errors. Availability: http://samtools.sourceforge.net. Contact: [email protected].

The published version (open access now) with addition of a test for multi-allelic sites

Cited by in corpus (14)

• Haplotype-based variant detection from short-read sequencing
• Towards Better Understanding of Artifacts in Variant Calling from High-Coverage Samples
• MOSAIK: A hash-based algorithm for accurate next-generation sequencing read mapping
• Joint assembly and genetic mapping of the Atlantic horseshoe crab genome reveals ancient whole genome duplication
• grenepipe: A flexible, scalable, and reproducible pipeline to automate variant and frequency calling from sequence reads
• BioWorkbench: A High-Performance Framework for Managing and Analyzing Bioinformatics Experiments
• Dynamic read mapping and online consensus calling for better variant detection
• A Simple Data-Adaptive Probabilistic Variant Calling Model
• A frame-based representation of genomic sequences for removing errors and rare variant detection in NGS data
• Chromosomal rearrangements and transposable elements in locally adapted island Drosophila
• Computer Architecture-Aware Optimisation of DNA Analysis Systems
• Computational Performance of a Germline Variant Calling Pipeline for Next Generation Sequencing
• WIP: Identifying Tutorial Affordances for Interdisciplinary Learning Environments
• De novo genomic analyses for non-model organisms: an evaluation of methods across a multi-species data set