Skip to main content
Glama

jean — Molecular Biology & Genomics

population_fst

Read-only

Differentiation between populations of a call set: Weir and Cockerham's θ over two or more populations, or Hudson's or Patterson's Fst between pairs, per site and averaged with a block-jackknife standard error (scikit-allel). Every pair when more than two populations are compared. Validated on every site of 1000 Genomes phase 3 APOL1 for the 5 super-populations and 26 populations.

Input Schema

TableJSON Schema
NameRequiredDescriptionDefault
vcfYesthe call set as VCF text (a .vcf, or a .vcf.gz already decompressed); every record's GT is read, ploidy 2.
_filesNoThe project: record references shaped like _subject. Their records become library documents a tool names in sequences/parts/document; a trace or alignment is found by file name.
compareNothe populations to compare, comma-separated; empty is every group of the column.
groupByNowhich column of the sample table groups the samples; empty is its first grouping column.
_subjectNoThe record to act on: exactly one of sequence, accession (+source), sample (an id from list_samples) or file ({name, text} or {name, base64}); record_index picks one record of several.
estimatorNoweirCockerham: θ from the variance components over every compared population at once. hudson: 1 − within / between pairwise differences (Bhatia 2013). patterson: F2-based, biallelic sites only. Values: weirCockerham (Weir–Cockerham θ), hudson (Hudson Fst), patterson (Patterson Fst).hudson
blockLengthNosites per jackknife block for the average's standard error; a trailing partial block is dropped, as scikit-allel does.
populationsYesa sample table: a header row with a `sample` column and one or more grouping columns (tab-separated, as 1000 Genomes' .panel). Samples it does not list belong to no group.

Output Schema

TableJSON Schema
NameRequiredDescriptionDefault
sitesYesper site: pos, and per comparison the site's Fst (numerator / denominator), null where undefined
matrixYespairwise average Fst between every compared population (Hudson and Patterson)
averageYesper comparison: the populations, Fst over all sites, its jackknife standard error, blocks

Schema Changelog

Changes observed during successful MCP inspections.

  1. First observed

TDQS

A3.7/5.0
Behavior4/5

Does the description disclose side effects, auth requirements, rate limits, or destructive behavior?

Annotations already carry the safety profile (readOnly, non-destructive, closed-world). Beyond that, the description discloses genuinely non-obvious behavior: results are produced per site AND averaged with a block-jackknife standard error, all pairs are enumerated when more than two populations are given, and it states the tool was validated on 1000 Genomes phase 3 APOL1. That added context exceeds what the annotations provide.

Agents need to know what a tool does to the world before calling it. Descriptions should go beyond structured annotations to explain consequences.

Conciseness4/5

Is the description appropriately sized, front-loaded, and free of redundancy?

The core purpose is front-loaded in the first clause, with method and scope details following in order of relevance. It is dense and somewhat run-on, and the trailing validation sentence is ancillary, but overall it is appropriately sized with little waste.

Shorter descriptions cost fewer tokens and are easier for agents to parse. Every sentence should earn its place.

Completeness4/5

Given the tool's complexity, does the description cover enough for an agent to succeed on first attempt?

With an output schema present, return values need not be explained, and annotations cover the safety profile. The description supplies estimator behavior and scope adequately for correct invocation, though it omits practical prerequisites (e.g. what sample-table columns are expected, genotype ploidy handling) that the schema only partially conveys.

Complex tools with many parameters or behaviors need more documentation. Simple tools need less. This dimension scales expectations accordingly.

Parameters3/5

Does the description clarify parameter syntax, constraints, interactions, or defaults beyond what the schema provides?

Schema description coverage is 100%, so the schema already documents vcf, compare, groupBy, populations, estimator, and blockLength in detail. The description reinforces that every pair is compared when >2 populations are supplied, but adds little parameter-level syntax or format meaning beyond the schema. Baseline 3 applies when the schema does the heavy lifting.

Input schemas describe structure but not intent. Descriptions should explain non-obvious parameter relationships and valid value ranges.

Purpose4/5

Does the description clearly state what the tool does and how it differs from similar tools?

The description names a specific analytic operation ('Differentiation between populations of a call set') plus the concrete statistics computed (Weir-Cockerham θ, Hudson's/Patterson's Fst). An agent can tell this is a population-genetics Fst tool, not a generic stats or ancestry tool. It does not explicitly differentiate itself from close siblings such as population_site_stats, population_ancestry, or genotype_pca.

Agents choose between tools based on descriptions. A clear purpose with a specific verb and resource helps agents select the right tool.

Usage Guidelines3/5

Does the description explain when to use this tool, when not to, or what alternatives exist?

Usage is only implied: it notes θ applies 'over two or more populations' and Hudson/Patterson work 'between pairs', which hints at estimator selection, but there is no explicit when-to-use/when-not guidance, no prerequisites (e.g. genotype-only VCF, ploidy 2), and no named alternatives. An agent must infer context from the schema's minimum requirements.

Agents often have multiple tools that could apply. Explicit usage guidance like "use X instead of Y when Z" prevents misuse.

Try in Browser

Glama MCP Gateway

Add one secure layer between your agents and this server.

Resources