Samtools is a mature, production-grade command-line toolkit for processing high-throughput sequencing alignment data in SAM, BAM, and CRAM formats. It provides dozens of subcommands covering the full NGS workflow: format conversion (view), sorting, indexing, duplicate marking (markdup), mate-pair fixing (fixmate), depth/coverage calculation, variant calling preparation (mpileup), consensus calling, amplicon analysis, and comprehensive statistics (stats/flagstat). The architecture is a single C binary built on the htslib library, with each subcommand implemented as a separate source file dispatched from bamtk.c. The codebase is well-structured with reusable abstractions for pileup iteration, BED interval indexing, temporary file management, and global option handling. It is actively maintained with extensive regression tests, cross-platform CI (Linux, macOS, Windows, FreeBSD, OpenBSD, Solaris), and a large user community in bioinformatics.
The Genome Analysis Toolkit (GATK4) is a production-grade bioinformatics framework developed by the Broad Institute for variant discovery and genotyping in high-throughput sequencing data. It provides a comprehensive suite of command-line tools for germline and somatic variant calling (HaplotypeCaller, Mutect2), copy number variation analysis (GermlineCNVCaller), structural variant detection (FindBreakpointEvidenceSpark), functional annotation (Funcotator), base quality score recalibration (BQSR), and read processing. The architecture is built around a plugin-style engine with walker base classes (ReadWalker, VariantWalker, LocusWalker) that define traversal patterns over genomic data, with extensive use of Apache Spark for distributed execution on clusters and cloud platforms (GCS, HDFS). Key entry points include the `gatk` wrapper script and individual tool classes registered via the Barclay argument parser, with workflows orchestrated through WDL/Cromwell pipelines. External dependencies include htsjdk for BAM/VCF I/O, BWA for alignment, PyMC3 for Bayesian CNV modeling, and PyTorch for deep learning-based variant scoring. The codebase is heavily tested with unit and integration tests, and supports both local and cloud-based execution.