{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,12,19]],"date-time":"2025-12-19T09:32:37Z","timestamp":1766136757519,"version":"3.37.3"},"reference-count":40,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2016,10,3]],"date-time":"2016-10-03T00:00:00Z","timestamp":1475452800000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2016,10,3]],"date-time":"2016-10-03T00:00:00Z","timestamp":1475452800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"DOI":"10.13039\/100000002","name":"National Institutes of Health","doi-asserted-by":"publisher","award":["CA118444"],"award-info":[{"award-number":["CA118444"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]},{"name":"Mayo CCaTS grant","award":["UL1TR000135"],"award-info":[{"award-number":["UL1TR000135"]}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["BMC Bioinformatics"],"abstract":"<jats:title>Abstract<\/jats:title><jats:sec>\n                <jats:title>Background<\/jats:title>\n                <jats:p>GATK Best Practices workflows are widely used in large-scale sequencing projects and recommend post-alignment processing before variant calling. Two key post-processing steps include the computationally intensive local realignment around known INDELs and base quality score recalibration (BQSR). Both have been shown to reduce erroneous calls; however, the findings are mainly supported by the analytical pipeline that incorporates BWA and GATK UnifiedGenotyper. It is not known whether there is any benefit of post-processing and to what extent the benefit might be for pipelines implementing other methods, especially given that both mappers and callers are typically updated. Moreover, because sequencing platforms are upgraded regularly and the new platforms provide better estimations of read quality scores, the need for post-processing is also unknown. Finally, some regions in the human genome show high sequence divergence from the reference genome; it is unclear whether there is benefit from post-processing in these regions.<\/jats:p>\n              <\/jats:sec><jats:sec>\n                <jats:title>Results<\/jats:title>\n                <jats:p>We used both simulated and NA12878 exome data to comprehensively assess the impact of post-processing for five or six popular mappers together with five callers. Focusing on chromosome 6p21.3, which is a region of high sequence divergence harboring the human leukocyte antigen (HLA) system, we found that local realignment had little or no impact on SNP calling, but increased sensitivity was observed in INDEL calling for the Stampy\u2009+\u2009GATK UnifiedGenotyper pipeline. No or only a modest effect of local realignment was detected on the three haplotype-based callers and no evidence of effect on Novoalign. BQSR had virtually negligible effect on INDEL calling and generally reduced sensitivity for SNP calling that depended on caller, coverage and level of divergence. Specifically, for SAMtools and FreeBayes calling in the regions with low divergence, BQSR reduced the SNP calling sensitivity but improved the precision when the coverage is insufficient. However, in regions of high divergence (e.g., the HLA region), BQSR reduced the sensitivity of both callers with little gain in precision rate. For the other three callers, BQSR reduced the sensitivity without increasing the precision rate regardless of coverage and divergence level.<\/jats:p>\n              <\/jats:sec><jats:sec>\n                <jats:title>Conclusions<\/jats:title>\n                <jats:p>We demonstrated that the gain from post-processing is not universal; rather, it depends on mapper and caller combination, and the benefit is influenced further by sequencing depth and divergence level. Our analysis highlights the importance of considering these key factors in deciding to apply the computationally intensive post-processing to Illumina exome data.<\/jats:p>\n              <\/jats:sec>","DOI":"10.1186\/s12859-016-1279-z","type":"journal-article","created":{"date-parts":[[2016,10,3]],"date-time":"2016-10-03T11:25:22Z","timestamp":1475493922000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":33,"title":["Impact of post-alignment processing in variant discovery from whole exome data"],"prefix":"10.1186","volume":"17","author":[{"given":"Shulan","family":"Tian","sequence":"first","affiliation":[]},{"given":"Huihuang","family":"Yan","sequence":"additional","affiliation":[]},{"given":"Michael","family":"Kalmbach","sequence":"additional","affiliation":[]},{"given":"Susan L.","family":"Slager","sequence":"additional","affiliation":[]}],"member":"297","published-online":{"date-parts":[[2016,10,3]]},"reference":[{"issue":"7422","key":"1279_CR1","doi-asserted-by":"publisher","first-page":"56","DOI":"10.1038\/nature11632","volume":"491","author":"GR Abecasis","year":"2012","unstructured":"Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, Kang HM, Marth GT, McVean GA. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56\u201365.","journal-title":"Nature"},{"key":"1279_CR2","doi-asserted-by":"publisher","first-page":"274","DOI":"10.1186\/1471-2105-14-274","volume":"14","author":"X Yu","year":"2013","unstructured":"Yu X, Sun S. Comparing a few SNP calling algorithms using low-coverage sequencing data. BMC Bioinformatics. 2013;14:274.","journal-title":"BMC Bioinformatics"},{"issue":"4","key":"1279_CR3","doi-asserted-by":"publisher","first-page":"357","DOI":"10.1038\/tpj.2015.61","volume":"16","author":"Q Xu","year":"2016","unstructured":"Xu Q, Wu X, Li M, Huang H, Minica C, Yi Z, Wang G, Shen L, Xing Q, Shi Y, et al. Association studies of genomic variants with treatment response to risperidone, clozapine, quetiapine and chlorpromazine in the Chinese Han population. Pharmacogenomics J. 2016;16(4):357\u201365.","journal-title":"Pharmacogenomics J"},{"issue":"5","key":"1279_CR4","doi-asserted-by":"publisher","first-page":"491","DOI":"10.1038\/ng.806","volume":"43","author":"MA DePristo","year":"2011","unstructured":"DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491\u20138.","journal-title":"Nat Genet"},{"issue":"12","key":"1279_CR5","doi-asserted-by":"publisher","first-page":"1350","DOI":"10.1038\/ng.3121","volume":"46","author":"NI Weisenfeld","year":"2014","unstructured":"Weisenfeld NI, Yin S, Sharpe T, Lau B, Hegarty R, Holmes L, Sogoloff B, Tabbaa D, Williams L, Russ C, et al. Comprehensive variation discovery in single human genomes. Nat Genet. 2014;46(12):1350\u20135.","journal-title":"Nat Genet"},{"key":"1279_CR6","first-page":"160","volume":"4","author":"Z Wang","year":"2013","unstructured":"Wang Z, Liu X, Yang BZ, Gelernter J. The role and challenges of exome sequencing in studies of human diseases. Front Genet. 2013;4:160.","journal-title":"Front Genet"},{"issue":"1","key":"1279_CR7","doi-asserted-by":"publisher","first-page":"5","DOI":"10.1038\/jhg.2013.114","volume":"59","author":"B Rabbani","year":"2014","unstructured":"Rabbani B, Tekin M, Mahdieh N. The promise of whole-exome sequencing in medical genetics. J Hum Genet. 2014;59(1):5\u201315.","journal-title":"J Hum Genet"},{"issue":"8","key":"1279_CR8","doi-asserted-by":"publisher","first-page":"912","DOI":"10.1038\/ng.3036","volume":"46","author":"A Rimmer","year":"2014","unstructured":"Rimmer A, Phan H, Mathieson I, Iqbal Z, Twigg SR, Wilkie AO, McVean G, Lunter G. Integrating mapping-, assembly- and haplotype-based approaches for calling variants in clinical sequencing applications. Nat Genet. 2014;46(8):912\u20138.","journal-title":"Nat Genet"},{"issue":"14","key":"1279_CR9","doi-asserted-by":"publisher","first-page":"1838","DOI":"10.1093\/bioinformatics\/bts280","volume":"28","author":"H Li","year":"2012","unstructured":"Li H. Exploring single-sample SNP and INDEL calling with whole-genome de novo assembly. Bioinformatics. 2012;28(14):1838\u201344.","journal-title":"Bioinformatics"},{"issue":"9","key":"1279_CR10","doi-asserted-by":"publisher","first-page":"e75619","DOI":"10.1371\/journal.pone.0075619","volume":"8","author":"X Liu","year":"2013","unstructured":"Liu X, Han S, Wang Z, Gelernter J, Yang BZ. Variant callers for next-generation sequencing data: a comparison study. PLoS One. 2013;8(9):e75619.","journal-title":"PLoS One"},{"key":"1279_CR11","doi-asserted-by":"publisher","first-page":"247","DOI":"10.1186\/1471-2105-15-247","volume":"15","author":"AM Meynert","year":"2014","unstructured":"Meynert AM, Ansari M, FitzPatrick DR, Taylor MS. Variant detection sensitivity and biases in whole genome and exome sequencing. BMC Bioinformatics. 2014;15:247.","journal-title":"BMC Bioinformatics"},{"issue":"20","key":"1279_CR12","doi-asserted-by":"publisher","first-page":"2843","DOI":"10.1093\/bioinformatics\/btu356","volume":"30","author":"H Li","year":"2014","unstructured":"Li H. Toward better understanding of artifacts in variant calling from high-coverage samples. Bioinformatics. 2014;30(20):2843\u201351.","journal-title":"Bioinformatics"},{"issue":"14","key":"1279_CR13","doi-asserted-by":"publisher","first-page":"1754","DOI":"10.1093\/bioinformatics\/btp324","volume":"25","author":"H Li","year":"2009","unstructured":"Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754\u201360.","journal-title":"Bioinformatics"},{"issue":"9","key":"1279_CR14","doi-asserted-by":"publisher","first-page":"1297","DOI":"10.1101\/gr.107524.110","volume":"20","author":"A McKenna","year":"2010","unstructured":"McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297\u2013303.","journal-title":"Genome Res"},{"key":"1279_CR15","unstructured":"Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. http:\/\/arxiv.org\/abs\/1207.3907v2. 2012."},{"key":"1279_CR16","doi-asserted-by":"publisher","first-page":"8","DOI":"10.1186\/1471-2105-13-8","volume":"13","author":"D Challis","year":"2012","unstructured":"Challis D, Yu J, Evani US, Jackson AR, Paithankar S, Coarfa C, Milosavljevic A, Gibbs RA, Yu F. An integrative variant analysis suite for whole exome next-generation sequencing data. BMC Bioinformatics. 2012;13:8.","journal-title":"BMC Bioinformatics"},{"issue":"6","key":"1279_CR17","doi-asserted-by":"publisher","first-page":"1124","DOI":"10.1101\/gr.088013.108","volume":"19","author":"R Li","year":"2009","unstructured":"Li R, Li Y, Fang X, Yang H, Wang J, Kristiansen K, Wang J. SNP detection for massively parallel whole-genome resequencing. Genome Res. 2009;19(6):1124\u201332.","journal-title":"Genome Res"},{"key":"1279_CR18","first-page":"11.10.1","volume":"11","author":"GA Van der Auwera","year":"2013","unstructured":"Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J, et al. From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;11:11.10.1\u201311.10.33.","journal-title":"Curr Protoc Bioinformatics"},{"key":"1279_CR19","doi-asserted-by":"publisher","first-page":"17875","DOI":"10.1038\/srep17875","volume":"5","author":"S Hwang","year":"2015","unstructured":"Hwang S, Kim E, Lee I, Marcotte EM. Systematic comparison of variant calling pipelines using gold standard personal exome variants. Sci Rep. 2015;5:17875.","journal-title":"Sci Rep"},{"key":"1279_CR20","doi-asserted-by":"crossref","unstructured":"Liu Q, Guo Y, Li J, Long J, Zhang B, Shyr Y. Steps to ensure accuracy in genotype and SNP calling from Illumina sequencing data. BMC Genomics. 2012;13 Suppl 8:S8.","DOI":"10.1186\/1471-2164-13-S8-S8"},{"issue":"8","key":"1279_CR21","doi-asserted-by":"publisher","first-page":"1157","DOI":"10.1093\/bioinformatics\/btr076","volume":"27","author":"H Li","year":"2011","unstructured":"Li H. Improving SNP, discovery by base alignment quality. Bioinformatics. 2011;27(8):1157\u20138.","journal-title":"Bioinformatics"},{"issue":"6","key":"1279_CR22","doi-asserted-by":"publisher","first-page":"443","DOI":"10.1038\/nrg2986","volume":"12","author":"R Nielsen","year":"2011","unstructured":"Nielsen R, Paul JS, Albrechtsen A, Song YS. Genotype and SNP calling from next-generation sequencing data. Nat Rev Genet. 2011;12(6):443\u201351.","journal-title":"Nat Rev Genet"},{"issue":"3","key":"1279_CR23","doi-asserted-by":"publisher","first-page":"186","DOI":"10.1101\/gr.8.3.186","volume":"8","author":"B Ewing","year":"1998","unstructured":"Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998;8(3):186\u201394.","journal-title":"Genome Res"},{"issue":"6","key":"1279_CR24","doi-asserted-by":"publisher","first-page":"936","DOI":"10.1101\/gr.111120.110","volume":"21","author":"G Lunter","year":"2011","unstructured":"Lunter G, Goodson M. Stampy: a statistical algorithm for sensitive and fast mapping of Illumina sequence reads. Genome Res. 2011;21(6):936\u20139.","journal-title":"Genome Res"},{"issue":"1","key":"1279_CR25","doi-asserted-by":"publisher","first-page":"15","DOI":"10.1016\/j.ygeno.2013.02.011","volume":"102","author":"CA McLure","year":"2013","unstructured":"McLure CA, Hinchliffe P, Lester S, Williamson JF, Millman JA, Keating PJ, Stewart BJ, Dawkins RL. Genomic evolution and polymorphism: segmental duplications and haplotypes at 108 regions on 21 chromosomes. Genomics. 2013;102(1):15\u201326.","journal-title":"Genomics"},{"issue":"11","key":"1279_CR26","doi-asserted-by":"publisher","first-page":"627","DOI":"10.1038\/nrg3933","volume":"16","author":"MJ Chaisson","year":"2015","unstructured":"Chaisson MJ, Wilson RK, Eichler EE. Genetic variation and the de novo assembly of human genomes. Nat Rev Genet. 2015;16(11):627\u201340.","journal-title":"Nat Rev Genet"},{"issue":"4","key":"1279_CR27","doi-asserted-by":"publisher","first-page":"843","DOI":"10.1182\/blood-2012-03-413591","volume":"120","author":"SL Slager","year":"2012","unstructured":"Slager SL, Skibola CF, Di Bernardo MC, Conde L, Broderick P, McDonnell SK, Goldin LR, Croft N, Holroyd A, Harris S, et al. Common variation at 6p21.31 (BAK1) influences the risk of chronic lymphocytic leukemia. Blood. 2012;120(4):843\u20136.","journal-title":"Blood"},{"issue":"22","key":"1279_CR28","doi-asserted-by":"publisher","first-page":"8676","DOI":"10.1073\/pnas.1206614109","volume":"109","author":"C Wang","year":"2012","unstructured":"Wang C, Krishnakumar S, Wilhelmy J, Babrzadeh F, Stepanyan L, Su LF, Levinson D, Fernandez-Vina MA, Davis RW, Davis MM, et al. High-throughput, high-fidelity HLA genotyping with deep sequencing. Proc Natl Acad Sci U S A. 2012;109(22):8676\u201381.","journal-title":"Proc Natl Acad Sci U S A"},{"key":"1279_CR29","doi-asserted-by":"publisher","first-page":"355","DOI":"10.1186\/1471-2164-14-355","volume":"14","author":"K Hosomichi","year":"2013","unstructured":"Hosomichi K, Jinam TA, Mitsunaga S, Nakaoka H, Inoue I. Phase-defined complete sequencing of the HLA genes by next-generation sequencing. BMC Genomics. 2013;14:355.","journal-title":"BMC Genomics"},{"key":"1279_CR30","doi-asserted-by":"publisher","first-page":"6275","DOI":"10.1038\/ncomms7275","volume":"6","author":"G Highnam","year":"2015","unstructured":"Highnam G, Wang JJ, Kusler D, Zook J, Vijayan V, Leibovich N, Mittelman D. An analytical framework for optimizing variant discovery from personal genomes. Nat Commun. 2015;6:6275.","journal-title":"Nat Commun"},{"issue":"4","key":"1279_CR31","doi-asserted-by":"publisher","first-page":"593","DOI":"10.1093\/bioinformatics\/btr708","volume":"28","author":"W Huang","year":"2012","unstructured":"Huang W, Li L, Myers JR, Marth GT. ART: a next-generation sequencing read simulator. Bioinformatics. 2012;28(4):593\u20134.","journal-title":"Bioinformatics"},{"issue":"7","key":"1279_CR32","doi-asserted-by":"publisher","first-page":"873","DOI":"10.1093\/bioinformatics\/btq057","volume":"26","author":"TD Wu","year":"2010","unstructured":"Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010;26(7):873\u201381.","journal-title":"Bioinformatics"},{"issue":"21","key":"1279_CR33","doi-asserted-by":"publisher","first-page":"2790","DOI":"10.1093\/bioinformatics\/btt468","volume":"29","author":"FJ Sedlazeck","year":"2013","unstructured":"Sedlazeck FJ, Rescheneder P, von Haeseler A. NextGenMap: fast and accurate read mapping in highly polymorphic genomes. Bioinformatics. 2013;29(21):2790\u20131.","journal-title":"Bioinformatics"},{"issue":"16","key":"1279_CR34","doi-asserted-by":"publisher","first-page":"2078","DOI":"10.1093\/bioinformatics\/btp352","volume":"25","author":"H Li","year":"2009","unstructured":"Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The Sequence Alignment\/Map format and SAMtools. Bioinformatics. 2009;25(16):2078\u20139.","journal-title":"Bioinformatics"},{"issue":"3","key":"1279_CR35","doi-asserted-by":"publisher","first-page":"246","DOI":"10.1038\/nbt.2835","volume":"32","author":"JM Zook","year":"2014","unstructured":"Zook JM, Chapman B, Wang J, Mittelman D, Hofmann O, Hide W, Salit M. Integrating human sequence data sets provides a resource of benchmark SNP and indel genotype calls. Nat Biotechnol. 2014;32(3):246\u201351.","journal-title":"Nat Biotechnol"},{"issue":"2","key":"1279_CR36","doi-asserted-by":"publisher","first-page":"226","DOI":"10.1038\/ng.1028","volume":"44","author":"Z Iqbal","year":"2012","unstructured":"Iqbal Z, Caccamo M, Turner I, Flicek P, McVean G. De novo assembly and genotyping of variants using colored de Bruijn graphs. Nat Genet. 2012;44(2):226\u201332.","journal-title":"Nat Genet"},{"key":"1279_CR37","unstructured":"Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:13033997. 2013."},{"issue":"7","key":"1279_CR38","doi-asserted-by":"publisher","first-page":"1006","DOI":"10.1093\/bioinformatics\/btt730","volume":"30","author":"H Zhao","year":"2014","unstructured":"Zhao H, Sun Z, Wang J, Huang H, Kocher JP, Wang L. CrossMap: a versatile tool for coordinate conversion between genome assemblies. Bioinformatics. 2014;30(7):1006\u20137.","journal-title":"Bioinformatics"},{"key":"1279_CR39","doi-asserted-by":"publisher","first-page":"e600","DOI":"10.7717\/peerj.600","volume":"2","author":"CD Warden","year":"2014","unstructured":"Warden CD, Adamson AW, Neuhausen SL, Wu X. Detailed comparison of two popular variant calling packages for exome and targeted exon studies. PeerJ. 2014;2:e600.","journal-title":"PeerJ"},{"issue":"15","key":"1279_CR40","doi-asserted-by":"publisher","first-page":"7217","DOI":"10.1093\/nar\/gkv677","volume":"43","author":"Y Jiang","year":"2015","unstructured":"Jiang Y, Turinsky AL, Brudno M. The missing indels: an estimate of indel variation in a human genome and analysis of factors that impede detection. Nucleic Acids Res. 2015;43(15):7217\u201328.","journal-title":"Nucleic Acids Res"}],"container-title":["BMC Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-016-1279-z.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1186\/s12859-016-1279-z\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-016-1279-z.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,2,1]],"date-time":"2024-02-01T18:11:01Z","timestamp":1706811061000},"score":1,"resource":{"primary":{"URL":"https:\/\/bmcbioinformatics.biomedcentral.com\/articles\/10.1186\/s12859-016-1279-z"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2016,10,3]]},"references-count":40,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2016,12]]}},"alternative-id":["1279"],"URL":"https:\/\/doi.org\/10.1186\/s12859-016-1279-z","relation":{},"ISSN":["1471-2105"],"issn-type":[{"type":"electronic","value":"1471-2105"}],"subject":[],"published":{"date-parts":[[2016,10,3]]},"assertion":[{"value":"2 May 2016","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"26 September 2016","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"3 October 2016","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}}],"article-number":"403"}}