Journal
HUMAN MUTATION
Volume 37, Issue 12, Pages 1263-1271Publisher
WILEY
DOI: 10.1002/humu.23114
Keywords
whole genome sequencing; whole exome sequencing; NGS; NA12878; alignment; variant calling; bioinformatics; computing speed; benchmark
Categories
Funding
- Spanish Ministry of Economy and Competitiveness
- Generalitat de Catalunya
- European Regional Development Fund (ERDF)
- RD-Connect Project (EC) [305444]
- ELIXIR-EXCELERATE (EC) [676559]
- MICINN [TIN2014-53234- C2-1-R]
Ask authors/readers for more resources
As whole genome sequencing becomes cheaper and faster, it will progressively substitute targeted next-generation sequencing as standard practice in research and diagnostics. However, computing cost-performance ratio is not advancing at an equivalent rate. Therefore, it is essential to evaluate the robustness of the variant detection process taking into account the computing resources required. We have benchmarked six combinations of state-of-the-art read aligners (BWA-MEM and GEM3) and variant callers (FreeBayes, GATK Haplotype-Caller, SAMtools) on whole genome and whole exome sequencing data from the NA12878 human sample. Results have been compared between them and against the NIST Genome in a Bottle (GIAB) variants reference dataset. We report differences in speed of up to 20 times in some steps of the process and have observed that SNV, and to a lesser extent InDel, detection is highly consistent in 70% of the genome. SNV, and especially InDel, detection is less reliable in 20% of the genome, and almost unfeasible in the remaining 10%. These findings will aid in choosing the appropriate tools bearing in mind objectives, workload, and computing infrastructure available. Published 2016 Wiley Periodicals, Inc.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available