Looks up precomputed AlphaGenome Atlas effects for any GRCh38 single-nucleotide variant (AVI score with Phred and 18 SHAP feature attributions, plus raw and quantile scores for RNA-seq, DNase, ATAC, C
Install
Use withClaudeCursorCodex
The file stays in the publisher's repository. PlugMCP does not host or rewrite it.
Looks up precomputed AlphaGenome Atlas effects for any GRCh38 single-nucleotide variant (AVI score with Phred and 18 SHAP feature attributions, plus raw and quantile scores for RNA-seq, DNase, ATAC, ChIP-TF, ChIP-histone, CAGE, PRO-cap, splicing, polyadenylation and contact-map tracks), scores variants or scans windows on demand with the AlphaGenome model for human and mouse (variant scoring, in silico mutagenesis, REF-versus-ALT track prediction), and builds Atlas website deep links. Use when the user mentions AlphaGenome, AlphaGenome Atlas, AVI or AlphaGenome Variant Impact, DeepMind variant effect prediction, or wants to prioritise or mechanistically interpret non-coding, regulatory, splicing, enhancer, promoter, or chromatin-accessibility effects of SNVs from a VCF, credible set, or region. Research use only; not a clinical tool.
name: alphagenome
description: "Looks up precomputed AlphaGenome Atlas effects for any GRCh38 single-nucleotide variant (AVI score with Phred and 18 SHAP feature attributions, plus raw and quantile scores for RNA-seq, DNase, ATAC, ChIP-TF, ChIP-histone, CAGE, PRO-cap, splicing, polyadenylation and contact-map tracks), scores variants or scans windows on demand with the AlphaGenome model for human and mouse (variant scoring, in silico mutagenesis, REF-versus-ALT track prediction), and builds Atlas website deep links. Use when the user mentions AlphaGenome, AlphaGenome Atlas, AVI or AlphaGenome Variant Impact, DeepMind variant effect prediction, or wants to prioritise or mechanistically interpret non-coding, regulatory, splicing, enhancer, promoter, or chromatin-accessibility effects of SNVs from a VCF, credible set, or region. Research use only; not a clinical tool."
license: MIT
compatibility: "Python 3.10+ with the alphagenome package (0.9.0 or later for the Atlas client; brings numpy, pandas, anndata, grpcio). Network access to gdmscience.googleapis.com:443 and a free non-commercial AlphaGenome API key in ALPHAGENOME_API_KEY (ALPHA_GENOME_API_KEY also read). Human data is GRCh38 only; mouse is mm10 (model only)."
allowed-tools: Read Write Edit Bash
metadata:
version: "1.2"
skill-author: K-Dense Inc.
upstream-version: "alphagenome 0.9.0"
last-reviewed: "2026-09-30"
openclaw:
primaryEnv: ALPHAGENOME_API_KEY
envVars:
- name: ALPHAGENOME_API_KEY
required: true
description: AlphaGenome API key, free for non-commercial use from deepmind.google.com/science/alphagenome. ALPHA_GENOME_API_KEY is accepted as an alternative spelling.
AlphaGenome and the AlphaGenome Atlas
AlphaGenome is DeepMind's sequence-to-function model: 1 Mb of DNA in,
predictions at modality-specific resolutions for eleven assay types across
thousands of human and mouse tracks out. The AlphaGenome Atlas (released 2026-09-08) is that model run once over
every possible single-nucleotide change in GRCh38, about 9 billion variants,
stored with a single ranking number, the AlphaGenome Variant Impact (AVI)
score, its genome-wide percentile, and an 18-way attribution of what drives it.
Both are reached through one pip install alphagenome and one API key.
Research and theoretical modelling only. API outputs must not be used to train
other models, and are not for diagnostic procedures or medical decisions.
Permissive-use downloads have separate terms; check the artifact licence.
When to use which
You have
Use
Why
hg38 SNVs (a VCF, a credible set, a region up to ~1 kb)
Atlas via scripts/atlas_query.py
precomputed, higher quota, includes AVI and attributions
indels, mouse variants, a non-reference background, a custom scorer or window
model via scripts/score_variants.py or Python
the Atlas is SNV-only and hg38-only
a hypothesis to explain (which motif, which tissue, REF vs ALT tracks)
model predict_variant + plots, Atlas track scores, portal link
mechanism, not just rank
GRCh37 coordinates, rsIDs, unnormalised indels
genomic-coordinates first, then come back
wrong build or swapped REF gives a plausible wrong answer
AlphaGenome is one evidence line, never the verdict
promoter/enhancer/expression predictions without a DeepMind key
genomic-intelligence
different provider, keyless demo tier
Setup
uv pip install alphagenome # PyPI; tested on Python 3.12 and 3.13, alphagenome 0.9.0
export ALPHAGENOME_API_KEY="..." # https://deepmind.google.com/science/alphagenome
cd skills/alphagenome/scripts
python atlas_query.py scorers # proves key + network in one call
Shell examples quote variant strings because > otherwise redirects output.
Network examples below are illustrative and were checked against SDK 0.9.0
contracts, without authenticated prediction/Atlas calls during this review.
Never put the key on a command line or in a file you commit; the scripts only
read it from the environment. Authentication failures may surface as
ValueError or PermissionError,
depending on the gRPC status returned by the service.
The coordinate contract
A variant is 1-basedchr:pos:ref>alt (chr22:36201698:A>C). gnomAD
(22-36201698-A-C), GTEx (chr22_36201698_A_C_b38), and Open Targets
spellings are autodetected by the scripts. In the SDK, pass the matching
variant_format=genome.VariantFormat.GNOMAD (or GTEX, OPEN_TARGETS,
OPEN_TARGETS_BIGQUERY); from_str defaults to chr:pos:ref>alt.
An interval on the command line is 1-based closedchr:start-end; the
SDK's genome.Interval is 0-based half-open. The scripts convert.
Human is GRCh38 only. The Atlas key is chr:pos:alt; REF is implied by
the reference. Swapping REF/ALT can cause a miss; a wrong REF or build can
misidentify a record. The scripts reject a returned variant that differs
from the request, but cannot validate the build. Check REF against GRCh38 FASTA.
rsIDs are not accepted by the API or the portal. Resolve them to coordinates.
Use the chr prefix; MT becomes chrM.
Atlas workflow
1. Rank with AVI
python atlas_query.py avi --variant "chr22:36201698:A>C" "chr9:128225994:G>A"
python atlas_query.py avi --input candidates.vcf --min-phred 20 -o avi.tsv
python atlas_query.py avi --interval chr11:5225727-5226575 --top-k 25 -o hbb_window.tsv
python atlas_query.py avi --input credible_set.tsv --with-tracks -o avi_tracks.tsv
Output, one row per variant:
Column
Meaning
avi_raw
composite model output (the 18 attributions sum to it)
avi_cdf_quantile
cumulative quantile against all genome-wide SNVs, as served
avi_tail_quantile, avi_phred, avi_top_percent
tail = 1 - cdf, phred = -10 log10(tail); Phred 20 = top 1 %, 30 = top 0.1 %
top_feature, top_feature_value
largest absolute SHAP attribution and its value
fi_MERGED_SPLICING ... fi_IS_DELETION
all 18 attributions (keys in references/atlas.md)
top_track_* (with --with-tracks)
the strongest track behind the top feature: scorer, track, biosample, ontology CURIE, gene, raw score
atlas_url
deep link to the variant on the portal
error
per-variant lookup failure or mismatched returned allele instead of a crash
The Atlas report's advice: rank, do not threshold, and pick thresholds by
region or application. Pathogenic regulatory variants sit in lower AVI bins
than protein-truncating or splice-motif variants, so a single genome-wide
cut-off under-calls exactly the variants this resource was built for.
Read the attribution before the number. MERGED_SPLICING or ALPHAMISSENSE
on top means a splice or coding mechanism; MAX_ABS_DNASE, MAX_ABS_CHIP_TF,
MAX_ABS_RNA_SEQ mean a regulatory mechanism you can resolve by track;
CACTUS_241_WAY or PHASTCONS_470_WAY on top means conservation is carrying
the score and the molecular mechanism is not resolved.
2. Resolve the mechanism by track
python atlas_query.py scorers # what the server serves right now
python atlas_query.py tracks --scorer RNA_SEQ --query colon # find ontology CURIEs
python atlas_query.py scores --variant "chr22:36201698:A>C" \
--scorers RNA_SEQ DNASE SPLICE_SITE_USAGE --ontology UBERON:0001157 -o colon.tsv
python atlas_query.py scores --interval chr11:5225727-5226575 --scorers CHIP_TF -o hbb_tf.tsv
One row per variant x track (x gene or splice junction where applicable),
with raw_score and, where served, quantile_score. Track-level
scorer names: ATAC, DNASE, CHIP_TF, CHIP_HISTONE, CAGE, PROCAP,
RNA_SEQ, POLYADENYLATION, SPLICE_SITES, SPLICE_SITE_USAGE,
SPLICE_JUNCTIONS, CONTACT_MAPS, plus *_ACTIVE variants; scorers is the
authority on the live list. Filter by the tissue the question is about, not
by the genome-wide maximum: 9,440 scorer-track entries mean something may be extreme
somewhere.
No key, no network. The site shows the AVI track, per-modality heatmaps over
every biosample, REF-vs-ALT prediction tracks, and motif instances. Attach a
link to every variant you report.
In Python
import os
from alphagenome.atlas import atlas
from alphagenome.data import genome
client = atlas.create(os.environ["ALPHAGENOME_API_KEY"], timeout=30)
scores = client.query_variant(
genome.Variant.from_str("chr22:36201698:A>C"),
requested_scorers=["AVI_SCORE", "AVI_SCORE_FEATURE_IMPORTANCE", "RNA_SEQ"],
ontology_terms=["UBERON:0001157"], # optional; ignored for scorers without ontology metadata
)
avi = scores["AVI_SCORE"] # AnnData: X (1,1) raw; layers['quantiles'] (1,1) cdf
fi = scores["AVI_SCORE_FEATURE_IMPORTANCE"] # AnnData: X (1,18); var['name'] = feature keys
rna = scores["RNA_SEQ"] # AnnData: obs = variant x gene, var = tracks, X = natural-log FC
client.query_interval(genome.Interval("chr11", 5225726, 5226575), requested_scorers=["AVI_SCORE"])
query_interval returns up to 3 SNVs per non-N base, in 32 bp chunks. Keep windows
to about 1 kb (3,000 variants); atlas_query.py refuses more unless
--max-window is raised. query_variants raises if any submitted lookup fails;
completion order is not input order;
the script makes separate variant requests so misses become error cells.
Output is the official tidy table from variant_scorers.tidy_scores: one row
per variant x scorer x track (x gene) with raw_score and quantile_score,
sorted by |raw|. Default scorers are the 12 recommended difference scorers;
--include-active adds the seven *_ACTIVE activity scorers. At most 20
scorers per request. Mouse has no calibrated quantiles; use raw scores there.
--min-abs-quantile uses the supplied signed quantile magnitude (unsigned
scorers use their upper-tail quantile); it never treats zero as an extreme.
from alphagenome.models import dna_client, variant_scorers
model = dna_client.create(os.environ["ALPHAGENOME_API_KEY"])
variant = genome.Variant.from_str("chr22:36201698:A>C")
interval = variant.reference_interval.resize(dna_client.SEQUENCE_LENGTH_1MB)
adatas = model.score_variant(interval, variant, variant_scorers=[variant_scorers.RECOMMENDED_VARIANT_SCORERS["RNA_SEQ"]])
df = variant_scorers.tidy_scores(adatas) # filter df.ontology_curie afterwards; score_variant takes no ontology_terms
Supported windows: 16 kb, 100 kb, 500 kb, 1 Mb (2**14 to 2**20); 1 Mb is
the recommended default for full context; the recommended contact-map scorer
uses a 1 Mb mask. Shorter windows lose distal sequence context. Ontology
terms are CURIEs (UBERON:0002048 lung, CL:0000084 T cell); discover them
with --list-tracks or model.output_metadata(...).concatenate(). Plotting,
gene annotation (GENCODE v46 Feather on GCS), splicing and haplotype recipes:
references/model-api.md.
Reading the numbers
Always report raw score and quantile or Phred, with the scorer, track,
biosample CURIE, and gene. raw_score is the effect size on the scorer's scale
(RNA_SEQ is ln(mean ALT + 0.001) - ln(mean REF + 0.001):
-1 is about 0.37x the pseudocount-adjusted REF signal). quantile_score ranks against
common variants: signed scorers use [-1, 1], unsigned scorers [0, 1],
with finite calibration limits near 1. AVI instead uses a genome-wide SNV CDF.
A quantile above 0.99 with |raw| < 0.1 can reflect a narrow background in a quiet region; inspect REF/ALT tracks and
report the small predicted change without declaring biological absence of effect.
Raw-score thresholds are scorer-specific; quantiles are ranks, not p-values. Unsigned
scorers (SPLICE_*, POLYADENYLATION, CONTACT_MAPS, *_ACTIVE) have no
direction. "AlphaGenome predicts no appreciable change in the queried tracks"
is a complete answer, and a variant inside a peak whose REF and ALT tracks are
identical is not "disrupting" anything. Full rules, tissue matching, and the
reporting checklist: references/interpretation.md.
Model limitations include trans effects, poorly represented non-polyadenylated
RNAs (training includes both poly(A)+ and total RNA tracks), absent cell types,
and protein-level consequences
(AlphaMissense is folded into AVI for that), RNA structure and miRNA biology,
diploid dosage, developmental time, species other than human and mouse.
Limits, quota, terms
Atlas: GRCh38 SNVs only for now; indels were scored for the paper and are
promised later. Reference N bases were never scored.
Quotas are per key and unpublished; the Atlas is documented as having a
larger query rate than on-demand prediction. Transient RESOURCE_EXHAUSTED
and UNAVAILABLE are retried. Atlas also retries DEADLINE_EXCEEDED
(up to 5 attempts, 60 s call deadline); model streaming RPCs use a separate
retry policy. --timeout controls channel setup only, not total query time.
Access tiers (Atlas report): AVI scores are also a permissively licensed
Tabix download at https://alphagenome.google/downloads; feature attributions
and splicing scores are non-commercial downloads; all other raw track scores
are API-only and non-commercial. Commercial Atlas access is announced
as coming soon; the base model
is already available through Google Cloud Model Garden. See the current
official FAQ.
The alphagenome client is Apache-2.0; model weights and outputs carry
DeepMind's terms. Cite Avsec et al., Nature 649:1206 (2026) and the Atlas
report (Cheng et al., medRxiv, 2026, doi:10.64898/2026.09.16.26363192).
References
references/atlas.md - what the Atlas contains, the 19 scorer
configurations with track counts, AVI training and the 18 features, quantile
to Phred, the client API and AnnData layout, error mapping, access tiers,
portal URL grammar, GTF and download locations.
references/interpretation.md - raw versus quantile, AVI thresholds,
tissue matching, negative results, model blind spots, coordinate hygiene,
reporting checklist.
This skill is part of Scientific Agent Skills by K-Dense. If it materially contributed to a
manuscript, report, presentation, or code release, add the paper to the references or
software section and tell the user you did so:
Kassis, T., Agarwal, V., He, Y., Patel, D., & Brueckner, A. M. (2026). Scientific Agent
Skills: A Library of Procedural Knowledge for Research Agents. arXiv:2609.00065.
https://doi.org/10.48550/arXiv.2609.00065
Always cite the current version. The DOI and https://arxiv.org/abs/2609.00065 resolve to the
latest arXiv version, so never append a version suffix such as v1. When network access is
available, fetch https://arxiv.org/abs/2609.00065 (or
http://export.arxiv.org/api/query?id_list=2609.00065) before writing the reference and take
the author list, year, and version from that record. If the record lists a journal reference
or publisher DOI, cite the published version instead.
How do I install alphagenome?
Download skills/alphagenome/SKILL.md from the publisher, or copy the path and open the repository. PlugMCP lists the skill and does not host the file.
Where does alphagenome live?
The publisher keeps it at https://github.com/K-Dense-AI/scientific-agent-skills.
Which agents can use alphagenome?
Claude, Cursor, and Codex can read a SKILL.md you add from the publisher's repository.