SNP calling, genotype calling, and sample allele frequency estimation from new-generation sequencing data

Research output: Contribution to journalJournal articleResearchpeer-review


We present a statistical framework for estimation and application of sample allele frequency spectra from New-Generation Sequencing (NGS) data. In this method, we first estimate the allele frequency spectrum using maximum likelihood. In contrast to previous methods, the likelihood function is calculated using a dynamic programming algorithm and numerically optimized using analytical derivatives. We then use a bayesian method for estimating the sample allele frequency in a single site, and show how the method can be used for genotype calling and SNP calling. We also show how the method can be extended to various other cases including cases with deviations from Hardy-Weinberg equilibrium. We evaluate the statistical properties of the methods using simulations and by application to a real data set.
Original languageEnglish
Article numbere37558
JournalPLoS ONE
Issue number7
Number of pages11
Publication statusPublished - 2012

Number of downloads are based on statistics from Google Scholar and

No data available

ID: 44047265