gapit

Mass screening of contigs and reads for known genes — AMR, virulence, serotype, plasmid replicons, species targets, or any custom database. An agent-first Python reimplementation of abricate: byte-compatible TSV plus first-class JSON and Markdown.

Why gapit

• Frozen abricate-compatible baseline. Contig screening keeps the same BLAST pipeline, hit rules, and abricate-format TSV as a frozen baseline surface for existing pipelines. From v0.5.0 gapit pursues its own contract; abricate parity is no longer a release gate — the parity harness (pixi run -e parity parity) remains a regression reference for that frozen surface.
• Machine-readable by design. Versioned output schemas (gapit.report/1, gapit.reads/1, gapit.summary/1), self-describing via gapit schema, typed JSON error envelopes on stderr, documented exit codes. An agent can discover the whole contract without reading docs.
• Reads, not just contigs. gapit screen --r1/--r2 screens FASTQ through minimap2, and --r1 accepts assembly FASTA directly (content-detected, map-ont forced) for a fast presence survey. abricate cannot screen raw reads.

Install

From PyPI (external BLAST+ and minimap2 binaries required on PATH — see docs/installation.md):

pip install gapit

Or from a git clone with pixi, which manages the environment including the external binaries:

git clone https://github.com/indexofire/gapit.git
cd gapit
pixi install
pixi run gapit --version

The package supports Python 3.11+; the dev environment pins 3.14.

Quick start

# Screen contigs (abricate-compatible TSV on stdout; --db is required)
gapit screen contigs.fa --db ncbi

# Agent- and human-readable outputs
gapit screen contigs.fa --db ncbi --format json
gapit screen contigs.fa --db ncbi --format md

# Two-stage designation: screen writes the gene table, typing designates from it
# (typed gene databases — the bundled ecoli_dec carries the dual-scheme DEC panel)
gapit screen isolates.fa --db ecoli_dec --output dec.tsv
gapit typing dec.tsv

# Screen FASTQ reads; --read-type picks the minimap2 preset
gapit screen --r1 sample_R1.fastq.gz --r2 sample_R2.fastq.gz --read-type sr --db card
gapit screen --r1 ont_reads.fastq.gz --db ncbi --read-type map-ont --format json

# Summarize report tables into a gene presence/absence matrix
gapit summary *.tsv

# Databases: six ship in-wheel and work instantly; the rest download from upstream on fetch
gapit db fetch all         # default set (card, vfdb)
gapit db fetch ncbi
gapit db list
gapit db outdated                # flag databases older than the staleness threshold
gapit db search "tet(M)"         # look up genes across every installed database
gapit db build mydb my_genes.fa --tsv my_meta.tsv   # custom db from any FASTA

# Introspection
gapit schema report     # JSON Schema of gapit.report/1

Contig screening follows abricate defaults (--minid 80, --mincov 80). Reads-mode presence defaults to 90% alignment breadth (--min-breadth 90).

Databases

A database is a directory under the datadir, resolved from $GAPIT_DATADIR, then ~/.local/share/gapit/db (override per call with --datadir). Twenty-two databases are built into the catalog — twelve gene providers, seven kaptive cluster providers, and three typed bundles (ecoli_dec, lm_doumith, and the cluster-kind vp_ok). Every provider builds on gapit db fetch, which downloads from upstream — and six content-provenance-audited ones also ship inside the package (ecoli_dec and lm_doumith plus the ncbi, resfinder, ecoh snapshots — public domain / Apache-2.0 / BSD-3-Clause / MIT — and the MIT typed-cluster snapshot vp_ok), materialized into the datadir on first use with zero network; the bundled ecoli_dec additionally carries per-gene identity floors (floors.json, pic at 90%) that reads mode applies before presence calls. Several upstream licenses (CARD's McMaster terms, VFDB's CC BY-NC, Kaptive's GPL-3.0) forbid redistribution inside an MIT-licensed distribution, so everything else is fetch-on-demand — and db fetch <name> refreshes the provider-backed bundled names (ecoh, ncbi, resfinder) from the latest upstream.

gapit also reads datadirs built by abricate itself (legacy ~~~ headers). The reverse does not hold: gapit-native databases use the gapit/v1 header format (see SPEC.md §11), which abricate cannot read. Protein databases such as bacmet2 screen through blastx.

Gene cluster databases are a second kind: gapit db build NAME loci.gbk|gff3 builds a kind: cluster db (locus calls via minimap2, optional --typing FILE phenotype scoring), and the seven Kaptive cluster databases (kpsc_k, kpsc_o, kosc_k, kosc_o, ab_k, ab_o, ecoli_kps) fetch the Kaptive antigen-locus references at install time (GPL-3.0 content, so nothing is bundled — cite Kaptive/Wyres et al. 2020 for results). The typed cluster bundle vp_ok merges the vpasero OAgc + CPSgc references (MIT) into one bundled db — a pure built-in in the ecoli_dec/lm_doumith pattern, no db fetch — where a single gapit screen --db vp_ok renders gapit.cluster/2 with o_group, k_group, and the composed serotype O:K call (RIMD 2210633 → O3:K6), working straight after pip install gapit with zero network. gapit db list --json surfaces each provider's content license where one is pinned.

Output contract

• stdout purity. Data on stdout, diagnostics on stderr, always. --quiet silences stderr only.
• Deterministic. Stable sort orders, fixed tool parameters, no wall-clock timestamps inside data payloads.
• Errors print a JSON envelope to stderr and exit nonzero:

{"schema": "gapit.error/1", "code": "...", "message": "...", "context": {...}}

MCP server

gapit ships an MCP (Model Context Protocol) stdio server so agent runtimes can screen assemblies without parsing CLI output: gapit mcp or the gapit-mcp console script speaks newline-delimited JSON-RPC 2.0 on stdin/stdout (no extra dependencies — the protocol is hand-rolled). It exposes nine tools: screen (gapit.report/1 by default; aligner minimap2 for a fast assembly survey), screen_reads (FASTQ via minimap2, gapit.reads/1), summary (gapit.summary/1), schema, db_list, plus the database tools db_fetch, db_build, db_search, and db_outdated so an agent can self-provision and inspect databases mid-session. Tool failures return isError: true with the gapit.error/1 envelope as text. Register it with an MCP client:

{"mcpServers": {"gapit": {"command": "gapit-mcp"}}}

Development

User documentation lives in docs/index.md, rendered at https://indexofire.github.io/gapit/.

SPEC.md is the parity contract, PLAN.md the roadmap, AGENTS.md the contributor guide, CHANGELOG.md the change history.

License

gapit is MIT-licensed. It is a behavioral reimplementation of abricate (GPL-2.0) and copies no Perl code; abricate itself remains GPL-2.0. The only database content inside the package is the six audited, permissively licensed bundles under src/gapit/data/dbs/ (public domain / Apache-2.0 / BSD-3-Clause / MIT); every other provider downloads from upstream at fetch time, under its own license (SPEC.md §9).