Accurate Variant Calling Research Articles

Founder reconstruction enables scalable and seamless pangenomic analysis.

MotivationVariant calling workflows that utilize a single reference sequence are the de facto standard elementary genomic analysis routine for resequencing projects. Various ways to enhance the reference with pangenomic information have been proposed, but scalability combined with seamless integration to existing workflows remains a challenge.ResultsWe present PanVC with founder sequences, a scalable and accurate variant calling workflow based on a multiple alignment of reference sequences. Scalability is achieved by removing duplicate parts up to a limit into a founder multiple alignment, that is then indexed using a hybrid scheme that exploits general purpose read aligners. Our implemented workflow uses GATK or BCFtools for variant calling, but the various steps of our workflow (e.g. vcf2multialign tool, founder reconstruction) can be of independent interest as a basis for creating novel pangenome analysis workflows beyond variant calling.Availability and implementationOur open access tools and instructions how to reproduce our experiments are available at the following address: https://github.com/algbio/panvc-founders.Supplementary information Supplementary data are available at Bioinformatics online.

Abstract 2282: Rapid and accurate variant calling of FFPE samples with the genexus system

Abstract Next Generation Sequencing technology has enhanced oncology research by enabling the detection of relevant variants for clinical and drug discovery programs. The Oncomine Comprehensive Assay v3, a pan-cancer panel, used with the Ion Torrent Genexus System allows for formalin-fixed paraffin embedded (FFPE) samples to be examined across 161 unique genes in an automated workflow from sample to result in 30hrs. This study demonstrates ≥ 90% accuracy for detecting SNV/Indel/CNV and Fusion variants. The Genexus System consists of two software linked instruments, the Genexus Purification System and the Genexus Integrated Sequencer. The purification system is used first to extract and quantify 6 FFPE sections from human colon and lung tissue for matching DNA & RNA pairs. The nucleic acid output plate from the purification instrument is directly transferred to the sequencer for library preparation and sequencing using the Oncomine Comprehensive Assay v3. Variant calling analysis of 6 DNA & RNA pairs are completed the following day. Results show that the variants were detected with ≥ 90% Sensitivity when compared to the Ion Torrent GeneStudio S5 System as an orthogonal method. All variants report p-values ≤ 10^-5 signifying more confidence that the variant call is correct. The Genexus System is a seamless solution for sample to variant calling results. Citation Format: Iris Casuga, Jason Gioia, Janice Au-Young, Frances Chan, Milton Huynh, Michael Crowe. Rapid and accurate variant calling of FFPE samples with the genexus system [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2282.

Association of mitochondrial variants and haplogroups identified by whole exome sequencing with Alzheimer's disease.

IntroductionFindings regarding the association between mitochondrial DNA (mtDNA) variants and Alzheimer's disease (AD) are inconsistent.MethodsWe developed a pipeline for accurate assembly and variant calling in mitochondrial genomes embedded within whole exome sequences (WES) from 10,831 participants from the Alzheimer's Disease Sequencing Project (ADSP). Association of AD risk was evaluated with each mtDNA variant and variants located in 1158 nuclear genes related to mitochondrial function using the SCORE test. Gene‐based tests were performed using SKAT‐O.ResultsAnalysis of 4220 mtDNA variants revealed study‐wide significant association of AD with a rare MT‐ND4L variant (rs28709356 C>T; minor allele frequency = 0.002; P = 7.3 × 10−5) as well as with MT‐ND4L in a gene‐based test (P = 6.71 × 10−5). Significant association was also observed with a MT‐related nuclear gene, TAMM41, in a gene‐based test (P = 2.7 × 10−5). The expression of TAMM41 was lower in AD cases than controls (P = .00046) or mild cognitive impairment cases (P = .03).DiscussionSignificant findings in MT‐ND4L and TAMM41 provide evidence for a role of mitochondria in AD.

ScSNV: accurate dscRNA-seq SNV co-expression analysis using duplicate tag collapsing

Identifying single nucleotide variants has become common practice for droplet-based single-cell RNA-seq experiments; however, presently, a pipeline does not exist to maximize variant calling accuracy. Furthermore, molecular duplicates generated in these experiments have not been utilized to optimally detect variant co-expression. Herein, we introduce scSNV designed from the ground up to “collapse” molecular duplicates and accurately identify variants and their co-expression. We demonstrate that scSNV is fast, with a reduced false-positive variant call rate, and enables the co-detection of genetic variants and A>G RNA edits across twenty-two samples.

Moss enables high sensitivity single-nucleotide variant calling from multiple bulk DNA tumor samples

Intra-tumor heterogeneity renders the identification of somatic single-nucleotide variants (SNVs) a challenging problem. In particular, low-frequency SNVs are hard to distinguish from sequencing artifacts. While the increasing availability of multi-sample tumor DNA sequencing data holds the potential for more accurate variant calling, there is a lack of high-sensitivity multi-sample SNV callers that utilize these data. Here we report Moss, a method to identify low-frequency SNVs that recur in multiple sequencing samples from the same tumor. Moss provides any existing single-sample SNV caller the ability to support multiple samples with little additional time overhead. We demonstrate that Moss improves recall while maintaining high precision in a simulated dataset. On multi-sample hepatocellular carcinoma, acute myeloid leukemia and colorectal cancer datasets, Moss identifies new low-frequency variants that meet manual review criteria and are consistent with the tumor’s mutational signature profile. In addition, Moss detects the presence of variants in more samples of the same tumor than reported by the single-sample caller. Moss’ improved sensitivity in SNV calling will enable more detailed downstream analyses in cancer genomics.

A unified haplotype-based method for accurate and comprehensive variant calling.

Almost all haplotype-based variant callers were designed specifically for detecting common germline variation in diploid populations, and give suboptimal results in other scenarios. Here we present Octopus, a variant caller that uses a polymorphic Bayesian genotyping model capable of modeling sequencing data from a range of experimental designs within a unified haplotype-aware framework. Octopus combines sequencing reads and prior information to phase-called genotypes of arbitrary ploidy, including those with somatic mutations. We show that Octopus accurately calls germline variants in individuals, including single nucleotide variants, indels and small complex replacements such as microinversions. Using a synthetic tumor data set derived from clean sequencing data from a sample with known germline haplotypes and observed mutations in a large cohort of tumor samples, we show that Octopus is more sensitive to low-frequency somatic variation, yet calls considerably fewer false positives than other methods. Octopus also outputs realigned evidence BAM files to aid validation and interpretation.

Single-cell DNA amplicon sequencing reveals clonal heterogeneity and evolution in T-cell acute lymphoblastic leukemia

AbstractT-cell acute lymphoblastic leukemia (T-ALL) is an aggressive leukemia that is most frequent in children and is characterized by the presence of few chromosomal rearrangements and 10 to 20 somatic mutations in protein-coding regions at diagnosis. The majority of T-ALL cases harbor activating mutations in NOTCH1 together with mutations in genes implicated in kinase signaling, transcriptional regulation, or protein translation. To obtain more insight in the level of clonal heterogeneity at diagnosis and during treatment, we used single-cell targeted DNA sequencing with the Tapestri platform. We designed a custom ALL panel and obtained accurate single-nucleotide variant and small insertion-deletion mutation calling for 305 amplicons covering 110 genes in about 4400 cells per sample and time point. A total of 108 188 cells were analyzed for 25 samples of 8 T-ALL patients. We typically observed a major clone at diagnosis (>35% of the cells) accompanied by several minor clones of which some were less than 1% of the total number of cells. Four patients had >2 NOTCH1 mutations, some of which present in minor clones, indicating a strong pressure to acquire NOTCH1 mutations in developing T-ALL cells. By analyzing longitudinal samples, we detected the presence and clonal nature of residual leukemic cells and clones with a minor presence at diagnosis that evolved to clinically relevant major clones at later disease stages. Therefore, single-cell DNA amplicon sequencing is a sensitive assay to detect clonal architecture and evolution in T-ALL.

Ratatosk: hybrid error correction of long reads enables accurate variant calling and assembly

A major challenge to long read sequencing data is their high error rate of up to 15%. We present Ratatosk, a method to correct long reads with short read data. We demonstrate on 5 human genome trios that Ratatosk reduces the error rate of long reads 6-fold on average with a median error rate as low as 0.22 %. SNP calls in Ratatosk corrected reads are nearly 99 % accurate and indel calls accuracy is increased by up to 37 %. An assembly of Ratatosk corrected reads from an Ashkenazi individual yields a contig N50 of 45 Mbp and less misassemblies than a PacBio HiFi reads assembly.

Parliament2: Accurate structural variant calling at scale.

BackgroundStructural variants (SVs) are critical contributors to genetic diversity and genomic disease. To predict the phenotypic impact of SVs, there is a need for better estimates of both the occurrence and frequency of SVs, preferably from large, ethnically diverse cohorts. Thus, the current standard approach requires the use of short paired-end reads, which remain challenging to detect, especially at the scale of hundreds to thousands of samples.FindingsWe present Parliament2, a consensus SV framework that leverages multiple best-in-class methods to identify high-quality SVs from short-read DNA sequence data at scale. Parliament2 incorporates pre-installed SV callers that are optimized for efficient execution in parallel to reduce the overall runtime and costs. We demonstrate the accuracy of Parliament2 when applied to data from NovaSeq and HiSeq X platforms with the Genome in a Bottle (GIAB) SV call set across all size classes. The reported quality score per SV is calibrated across different SV types and size classes. Parliament2 has the highest F1 score (74.27%) measured across the independent gold standard from GIAB. We illustrate the compute performance by processing all 1000 Genomes samples (2,691 samples) in <1 day on GRCH38. Parliament2 improves the runtime performance of individual methods and is open source (https://github.com/slzarate/parliament2), and a Docker image, as well as a WDL implementation, is available.ConclusionParliament2 provides both a highly accurate single-sample SV call set from short-read DNA sequence data and enables cost-efficient application over cloud or cluster environments, processing thousands of samples.

Read trimming has minimal effect on bacterial SNP-calling accuracy.

Read alignment is the central step of many analytic pipelines that perform variant calling. To reduce error, it is common practice to pre-process raw sequencing reads to remove low-quality bases and residual adapter contamination, a procedure collectively known as ‘trimming’. Trimming is widely assumed to increase the accuracy of variant calling, although there are relatively few systematic evaluations of its effects and no clear consensus on its efficacy. As sequencing datasets increase both in number and size, it is worthwhile reappraising computational operations of ambiguous benefit, particularly when the scope of many analyses now routinely incorporates thousands of samples, increasing the time and cost required. Using a curated set of 17 Gram-negative bacterial genomes, this study initially evaluated the impact of four read-trimming utilities (Atropos, fastp, Trim Galore and Trimmomatic), each used with a range of stringencies, on the accuracy and completeness of three bacterial SNP-calling pipelines. It was found that read trimming made only small, and statistically insignificant, increases in SNP-calling accuracy even when using the highest-performing pre-processor in this study, fastp. To extend these findings, >6500 publicly archived sequencing datasets from Escherichia coli , Mycobacterium tuberculosis and Staphylococcus aureus were re-analysed using a common analytic pipeline. Of the approximately 125 million SNPs and 1.25 million indels called across all samples, the same bases were called in 98.8 and 91.9 % of cases, respectively, irrespective of whether raw reads or trimmed reads were used. Nevertheless, the proportion of mixed calls (i.e. calls where <100 % of the reads support the variant allele; considered a proxy of false positives) was significantly reduced after trimming, which suggests that while trimming rarely alters the set of variant bases, it can affect the proportion of reads supporting each call. It was concluded that read quality- and adapter-trimming add relatively little value to a SNP-calling pipeline and may only be necessary if small differences in the absolute number of SNP calls, or the false call rate, are critical. Broadly similar conclusions can be drawn about the utility of trimming to an indel-calling pipeline. Read trimming remains routinely performed prior to variant calling likely out of concern that doing otherwise would typically have negative consequences. While historically this may have been the case, the data in this study suggests that read trimming is not always a practical necessity.

CROMqs: An infinitesimal successive refinement lossy compressor for the quality scores.

The amount of sequencing data is growing at a fast pace due to a rapid revolution in sequencing technologies. Quality scores, which indicate the reliability of each of the called nucleotides, take a significant portion of the sequencing data. In addition, quality scores are more challenging to compress than nucleotides, and they are often noisy. Hence, a natural solution to further decrease the size of the sequencing data is to apply lossy compression to the quality scores. Lossy compression may result in a loss in precision, however, it has been shown that when operating at some specific rates, lossy compression can achieve performance on variant calling similar to that achieved with the losslessly compressed data (i.e. the original data). We propose Coding with Random Orthogonal Matrices for quality scores (CROMqs), the first lossy compressor designed for the quality scores with the "infinitesimal successive refinability" property. With this property, the encoder needs to compress the data only once, at a high rate, while the decoder can decompress it iteratively. The decoder can reconstruct the set of quality scores at each step with reduced distortion each time. This characteristic is specifically useful in sequencing data compression, since the encoder does not generally know what the most appropriate rate of compression is, e.g. for not degrading variant calling accuracy. CROMqs avoids the need of having to compress the data at multiple rates, hence incurring time savings. In addition to this property, we show that CROMqs obtains a comparable rate-distortion performance to the state-of-the-art lossy compressors. Moreover, we also show that it achieves a comparable performance on variant calling to that of the lossless compressed data while achieving more than 50% reduction in size.

A customized scaffolds approach for the detection and phasing of complex variants by next-generation sequencing

Next-generation sequencing (NGS) is widely used in genetic testing for the highly sensitive detection of single nucleotide changes and small insertions or deletions. However, detection and phasing of structural variants, especially in repetitive or homologous regions, can be problematic due to uneven read coverage or genome reference bias, resulting in false calls. To circumvent this challenge, a computational approach utilizing customized scaffolds as supplementary reference sequences for read alignment was developed, and its effectiveness demonstrated with two CBS gene variants: NM_000071.2:c.833T>C and NM_000071.2:c.[833T>C; 844_845ins68]. Variant c.833T>C is a known causative mutation for homocystinuria, but is not pathogenic when in cis with the insertion, c.844_845ins68, because of alternative splicing. Using simulated reads, the custom scaffolds method resolved all possible combinations with 100% accuracy and, based on > 60,000 clinical specimens, exceeded the performance of current approaches that only align reads to GRCh37/hg19 for the detection of c.833T>C alone or in cis with c.844_845ins68. Furthermore, analysis of two 1000 Genomes Project trios revealed that the c.[833T>C; 844_845ins68] complex variant had previously been undetected in these datasets, likely due to the alignment method used. This approach can be configured for existing workflows to detect other challenging and potentially underrepresented variants, thereby augmenting accurate variant calling in clinical NGS testing.

Abstract 2505: The Oncomine BRCA expanded next-generation sequencing assay: Development and analytical validation of a new panel for detection of SNV, insertions, deletions and copy number variants in a panel of 15 genes involved in homologous recombination repair of double-strand break DNA damage

Abstract Introduction: Many tumors display phenotypic similarity across tissue types, characterized by defects in biological pathways leading to susceptibility to particular therapies. Deficiencies in one such pathway, homologous recombination DNA damage repair (HR DDR), affect the response of breast and ovarian cancers to agents that inhibit PARP function. The Oncomine BRCA Expanded panel targets 15 genes (BRCA1 and 2, ATM, BARD1, BRIP1, CDK12, CHEK2, FANCD2, MRE11, NBN, PALB2, PPP2R2A, RAD51B, RAD54L, and TP53) that play various roles in the HR DDR pathway. The assay includes 1011 amplicons covering both hotspot variants and coding regions of these genes. An additional set of pre-qualified genes is available for assay customization, including a recommended set of 25 genes of particular relevance to HR DDR pathway. Methods: Targeted regions were selected, and amplicons were designed, optimized, and experimentally verified to perform well with semi-conductor based sequencing. Measures used included gene-base and hotspot-base coverage uniformity, variant calling accuracy, and detection of known CNV gains in both FFPE and cell line DNA samples. Large gene rearrangements in BRCA1 and BRCA2, including exon deletions and duplications, could be detected by using a custom exon deletion pipeline. Analysis was performed on an integrated bioinformatics platform, Ion Reporter, which includes variant calling, identification and annotation of the likely deleterious mutations, and production of configurable reports, including visualization of the extent of copy number changes. Results: The assay's performance was evaluated using a variety of sample types, including commercially available controls and FFPE samples for breast, ovarian, prostate, and other tumor types. Sensitivity and positive predictive value for small variant detection was &amp;gt;95% in Acrometrix Hotspot Controls, and &amp;gt;90% of CNV gains greater than 6 were detected in cell line samples. Examples of all detected variant classes were displayed based on the Ion Reporter user interface. Conclusions: We introduce an expanded BRCA panel for researching genetic changes influencing homologous recombination DNA damage repair. The option to customize the gene content of the panel is enabled via a web-based interface to inventoried, pre-qualified targets. This approach allows identifying a broad set of mutation types in a single assay with high accuracy, enabling interrogation of HR DDR genes in FFPE preparations. For research use only. Not for use in diagnostic procedures. Citation Format: Charles Scafe, Yun Zhu, Chenchen Yang, Yu-Ting Tseng, Brooke McKnight, Fernando Farfan, Santoshi Bandla, Seth Sadis, Steven Roman, Fiona C. Hyland. The Oncomine BRCA expanded next-generation sequencing assay: Development and analytical validation of a new panel for detection of SNV, insertions, deletions and copy number variants in a panel of 15 genes involved in homologous recombination repair of double-strand break DNA damage [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 2505.

Accurate fetal variant calling in the presence of maternal cell contamination.

High-throughput sequencing of fetal DNA is a promising and increasingly common method for the discovery of all (or all coding) genetic variants in the fetus, either as part of prenatal screening or diagnosis, or for genetic diagnosis of spontaneous abortions. In many cases, the fetal DNA (from chorionic villi, amniotic fluid, or abortive tissue) can be contaminated with maternal cells, resulting in the mixture of fetal and maternal DNA. This maternal cell contamination (MCC) undermines the assumption, made by traditional variant callers, that each allele in a heterozygous site is covered, on average, by 50% of the reads, and therefore can lead to erroneous genotype calls. We present a panel of methods for reducing the genotyping error in the presence of MCC. All methods start with the output of GATK HaplotypeCaller on the sequencing data for the (contaminated) fetal sample and both of its parents, and additionally rely on information about the MCC fraction (which itself is readily estimated from the high-throughput sequencing data). The first of these methods uses a Bayesian probabilistic model to correct the fetal genotype calls produced by MCC-unaware HaplotypeCaller. The other two methods "learn" the genotype-correction model from examples. We use simulated contaminated fetal data to train and test the models. Using the test sets, we show that all three methods lead to substantially improved accuracy when compared with the original MCC-unaware HaplotypeCaller calls. We then apply the best-performing method to three chorionic villus samples from spontaneously terminated pregnancies.

Numt identification and removal with RtN!

Assays in mitochondrial genomics rely on accurate read mapping and variant calling. However, there are known and unknown nuclear paralogs that have fundamentally different genetic properties than that of the mitochondrial genome. Such paralogs complicate the interpretation of mitochondrial genome data and confound variant calling. Remove the Numts! (RtN!) was developed to categorize reads from massively parallel sequencing data not based on the expected properties and sequence identities of paralogous nuclear encoded mitochondrial sequences, but instead using sequence similarity to a large database of publicly available mitochondrial genomes. RtN! removes low-level sequencing noise and mitochondrial paralogs while not impacting variant calling, while competing methods were shown to remove true variants from mitochondrial mixtures. https://github.com/Ahhgust/RtN. Supplementary data are available at Bioinformatics online.

Dv-trio: a family-based variant calling pipeline using DeepVariant.

In 2018, Google published an innovative variant caller, DeepVariant, which converts pileups of sequence reads into images and uses a deep neural network to identify single-nucleotide variants and small insertion/deletions from next-generation sequencing data. This approach outperforms existing state-of-the-art tools. However, DeepVariant was designed to call variants within a single sample. In disease sequencing studies, the ability to examine a family trio (father-mother-affected child) provides greater power for disease mutation discovery. To further improve DeepVariant's variant calling accuracy in family-based sequencing studies, we have developed a family-based variant calling pipeline, dv-trio, which incorporates the trio information from the Mendelian genetic model into variant calling based on DeepVariant. dv-trio is available via an open source BSD3 license at GitHub (https://github.com/VCCRI/dv-trio/). e.giannoulatou@victorchang.edu.au. Supplementary data are available at Bioinformatics online.

A novel virtual barcode strategy for accurate panel-wide variant calling in circulating tumor DNA

BackgroundHybrid capture-based next-generation sequencing of DNA has been widely applied in the detection of circulating tumor DNA (ctDNA). Various methods have been proposed for ctDNA detection, but low-allelic-fraction (AF) variants are still a great challenge. In addition, no panel-wide calling algorithm is available, which hiders the full usage of ctDNA based ‘liquid biopsy’. Thus, we developed the VBCALAVD (Virtual Barcode-based Calling Algorithm for Low Allelic Variant Detection) in silico to overcome these limitations.ResultsBased on the understanding of the nature of ctDNA fragmentation, a novel platform-independent virtual barcode strategy was established to eliminate random sequencing errors by clustering sequencing reads into virtual families. Stereotypical mutant-family-level background artifacts were polished by constructing AF distributions. Three additional robust fine-tuning filters were obtained to eliminate stochastic mutant-family-level noises. The performance of our algorithm was validated using cell-free DNA reference standard samples (cfDNA RSDs) and normal healthy cfDNA samples (cfDNA controls). For the RSDs with AFs of 0.1, 0.2, 0.5, 1 and 5%, the mean F1 scores were 0.43 (0.25~0.56), 0.77, 0.92, 0.926 (0.86~1.0) and 0.89 (0.75~1.0), respectively, which indicates that the proposed approach significantly outperforms the published algorithms. Among controls, no false positives were detected. Meanwhile, characteristics of mutant-family-level noise and quantitative determinants of divergence between mutant-family-level noises from controls and RSDs were clearly depicted.ConclusionsDue to its good performance in the detection of low-AF variants, our algorithm will greatly facilitate the noninvasive panel-wide detection of ctDNA in research and clinical settings. The whole pipeline is available at https://github.com/zhaodalv/VBCALAVD.

Exploring the limit of using a deep neural network on pileup data for germline variant calling

Single-molecule sequencing technologies have emerged in recent years and revolutionized structural variant calling, complex genome assembly and epigenetic mark detection. However, the lack of a highly accurate small variant caller has limited these technologies from being more widely used. Here, we present Clair, the successor to Clairvoyante, a program for fast and accurate germline small variant calling, using single-molecule sequencing data. For Oxford Nanopore Technology data, Clair achieves better precision, recall and speed than several competing programs, including Clairvoyante, Longshot and Medaka. Through studying the missed variants and benchmarking intentionally overfitted models, we found that Clair may be approaching the limit of possible accuracy for germline small variant calling using pileup data and deep neural networks. Clair requires only a conventional central processing unit (CPU) for variant calling and is an open-source project available at https://github.com/HKU-BAL/Clair. A lack of accurate and efficient variant calling methods has held back single-molecule sequencing technologies from clinical applications. The authors present a deep-learning method for fast and accurate germline small variant calling, using single-molecule sequencing data.

NanoVar: accurate characterization of patients\u2019 genomic structural variants using low-depth nanopore sequencing

The recent advent of third-generation sequencing technologies brings promise for better characterization of genomic structural variants by virtue of having longer reads. However, long-read applications are still constrained by their high sequencing error rates and low sequencing throughput. Here, we present NanoVar, an optimized structural variant caller utilizing low-depth (8X) whole-genome sequencing data generated by Oxford Nanopore Technologies. NanoVar exhibits higher structural variant calling accuracy when benchmarked against current tools using low-depth simulated datasets. In patient samples, we successfully validate structural variants characterized by NanoVar and uncover normal alternative sequences or alleles which are present in healthy individuals.

Contaminant DNA in bacterial sequencing experiments is a major source of false genetic variability

BackgroundContaminant DNA is a well-known confounding factor in molecular biology and in genomic repositories. Strikingly, analysis workflows for whole-genome sequencing (WGS) data commonly do not account for errors potentially introduced by contamination, which could lead to the wrong assessment of allele frequency both in basic and clinical research.ResultsWe used a taxonomic filter to remove contaminant reads from more than 4000 bacterial samples from 20 different studies and performed a comprehensive evaluation of the extent and impact of contaminant DNA in WGS. We found that contamination is pervasive and can introduce large biases in variant analysis. We showed that these biases can result in hundreds of false positive and negative SNPs, even for samples with slight contamination. Studies investigating complex biological traits from sequencing data can be completely biased if contamination is neglected during the bioinformatic analysis, and we demonstrate that removing contaminant reads with a taxonomic classifier permits more accurate variant calling. We used both real and simulated data to evaluate and implement reliable, contamination-aware analysis pipelines.ConclusionAs sequencing technologies consolidate as precision tools that are increasingly adopted in the research and clinical context, our results urge for the implementation of contamination-aware analysis pipelines. Taxonomic classifiers are a powerful tool to implement such pipelines.