Twelve Years of SAMtools and BCFtools

Overview

Journal Gigascience

Publisher Oxford University Press

Specialties Biology
Genetics

Date 2021 Feb 16

PMID 33590861

Citations 4014

Authors

Petr Danecek

James K Bonfield

Jennifer Liddle

John Marshall

Valeriu Ohan

Martin O Pollard

Andrew Whitwham

Thomas Keane

Shane A McCarthy

Robert M Davies

Heng Li

Affiliations

Soon will be listed here.

Abstract

Background: SAMtools and BCFtools are widely used programs for processing and analysing high-throughput sequencing data. They include tools for file format conversion and manipulation, sorting, querying, statistics, variant calling, and effect analysis amongst other methods.

Findings: The first version appeared online 12 years ago and has been maintained and further developed ever since, with many new features and improvements added over the years. The SAMtools and BCFtools packages represent a unique collection of tools that have been used in numerous other software projects and countless genomic pipelines.

Conclusion: Both SAMtools and BCFtools are freely available on GitHub under the permissive MIT licence, free for both non-commercial and commercial use. Both packages have been installed >1 million times via Bioconda. The source code and documentation are available from https://www.htslib.org.

Citing Articles

Next-generation sequencing-based population genetics unravels the evolutionary history of Rhodomyrtus tomentosa in China.

Xu X, Liao B, Liao S, Qin Q, He C, Ding X BMC Plant Biol. 2025; 25(1):338.

PMID: 40089704 DOI: 10.1186/s12870-025-06364-6.

Direct measurement of the male germline mutation rate in individuals using sequential sperm samples.

Shoag J, Srinivasa A, Loh C, Liu M, Lassen E, Melanaphy S Nat Commun. 2025; 16(1):2546.

PMID: 40089484 DOI: 10.1038/s41467-025-57507-0.

Benchmarking long-read structural variant calling tools and combinations for detecting somatic variants in cancer genomes.

Aydin S, Yilmaz K, Acar A Sci Rep. 2025; 15(1):8707.

PMID: 40082509 PMC: 11906795. DOI: 10.1038/s41598-025-92750-x.

Integration of therapeutic cargo into the human genome with programmable type V-K CAST.

Liu J, Aliaga Goltsman D, Alexander L, Khayi K, Hong J, Dunham D Nat Commun. 2025; 16(1):2427.

PMID: 40082411 PMC: 11906591. DOI: 10.1038/s41467-025-57416-2.

The genome sequence of a flea beetle, (Marsham, 1802).

Geiser M, Sims I Wellcome Open Res. 2025; 10:62.

PMID: 40078959 PMC: 11897695. DOI: 10.12688/wellcomeopenres.23697.1.

References

Danecek P, McCarthy S, Durbin R . A Method for Checking Genomic Integrity in Cultured Cell Lines from SNP Genotyping Data. PLoS One. 2016; 11(5):e0155014. PMC: 4866717. DOI: 10.1371/journal.pone.0155014. View

Schilbert H, Rempel A, Pucker B . Comparison of Read Mapping and Variant Calling Tools for the Analysis of Plant NGS Data. Plants (Basel). 2020; 9(4). PMC: 7238416. DOI: 10.3390/plants9040439. View

Danecek P, Bonfield J, Liddle J, Marshall J, Ohan V, Pollard M . Twelve years of SAMtools and BCFtools. Gigascience. 2021; 10(2). PMC: 7931819. DOI: 10.1093/gigascience/giab008. View

DePristo M, Banks E, Poplin R, Garimella K, Maguire J, Hartl C . A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011; 43(5):491-8. PMC: 3083463. DOI: 10.1038/ng.806. View

Pightling A, Petronella N, Pagotto F . Choice of reference-guided sequence assembler and SNP caller for analysis of Listeria monocytogenes short-read sequence data greatly influences rates of error. BMC Res Notes. 2015; 8:748. PMC: 4672502. DOI: 10.1186/s13104-015-1689-4. View

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N . The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25(16):2078-9. PMC: 2723002. DOI: 10.1093/bioinformatics/btp352. View

Liu F, Zhang Y, Zhang L, Li Z, Fang Q, Gao R . Systematic comparative analysis of single-nucleotide variant detection methods from single-cell RNA sequencing data. Genome Biol. 2019; 20(1):242. PMC: 6862814. DOI: 10.1186/s13059-019-1863-4. View

Bonfield J, Marshall J, Danecek P, Li H, Ohan V, Whitwham A . HTSlib: C library for reading/writing high-throughput sequencing data. Gigascience. 2021; 10(2). PMC: 7931820. DOI: 10.1093/gigascience/giab007. View

McLaren W, Gil L, Hunt S, Riat H, Ritchie G, Thormann A . The Ensembl Variant Effect Predictor. Genome Biol. 2016; 17(1):122. PMC: 4893825. DOI: 10.1186/s13059-016-0974-4. View

10.

Narasimhan V, Danecek P, Scally A, Xue Y, Tyler-Smith C, Durbin R . BCFtools/RoH: a hidden Markov model approach for detecting autozygosity from next-generation sequencing data. Bioinformatics. 2016; 32(11):1749-51. PMC: 4892413. DOI: 10.1093/bioinformatics/btw044. View

11.

Danecek P, McCarthy S . BCFtools/csq: haplotype-aware variant consequences. Bioinformatics. 2017; 33(13):2037-2039. PMC: 5870570. DOI: 10.1093/bioinformatics/btx100. View

12.

Li H . A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011; 27(21):2987-93. PMC: 3198575. DOI: 10.1093/bioinformatics/btr509. View

13.

Danecek P, Auton A, Abecasis G, Albers C, Banks E, DePristo M . The variant call format and VCFtools. Bioinformatics. 2011; 27(15):2156-8. PMC: 3137218. DOI: 10.1093/bioinformatics/btr330. View

14.

Bonfield J, McCarthy S, Durbin R . Crumble: reference free lossy compression of sequence quality values. Bioinformatics. 2018; 35(2):337-339. PMC: 6330002. DOI: 10.1093/bioinformatics/bty608. View

15.

Hsi-Yang Fritz M, Leinonen R, Cochrane G, Birney E . Efficient storage of high throughput DNA sequencing data using reference-based compression. Genome Res. 2011; 21(5):734-40. PMC: 3083090. DOI: 10.1101/gr.114819.110. View

16.

Gruning B, Dale R, Sjodin A, Chapman B, Rowe J, Tomkins-Tinch C . Bioconda: sustainable and comprehensive software distribution for the life sciences. Nat Methods. 2018; 15(7):475-476. PMC: 11070151. DOI: 10.1038/s41592-018-0046-7. View

17.

Auton A, Brooks L, Durbin R, Garrison E, Kang H, Korbel J . A global reference for human genetic variation. Nature. 2015; 526(7571):68-74. PMC: 4750478. DOI: 10.1038/nature15393. View

18.

Li H . Improving SNP discovery by base alignment quality. Bioinformatics. 2011; 27(8):1157-8. PMC: 3072548. DOI: 10.1093/bioinformatics/btr076. View

19.

Tarasov A, Vilella A, Cuppen E, Nijman I, Prins P . Sambamba: fast processing of NGS alignment formats. Bioinformatics. 2015; 31(12):2032-4. PMC: 4765878. DOI: 10.1093/bioinformatics/btv098. View