Skip to main content
Glama

jean — Molecular Biology & Genomics

discover_site_avoidance

Read-only

Sweep a set of sequences against every exact recognition site in the NEB catalogue and report the sites each one carries far less (avoided) or far more (enriched) often than a stated null predicts, with the false discovery rate controlled across every (sequence, site) test. Avoidance is evidence of the restriction–modification systems a lineage has lived with, and a site a design can rely on being rare. Three nulls: base composition (order-0, Poisson); k-let-preserving shuffles (uShuffle), which holds dinucleotide skews such as CpG depletion in the null rather than reporting them as findings; and R'MES's maximal-order Markov model, the standard in published avoidance analyses, which judges a site of length k against the record's own (k−1)-mers so only what those cannot explain is reported.

Input Schema

TableJSON Schema
NameRequiredDescriptionDefault
seedNoseed for the shuffles, so a run reproduces. Only applies when null_model is shuffle. Advanced.
alphaNofalse discovery rate at which a site is reported avoided or enriched.
_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.
enzymesNorestrict the sweep to these enzymes' sites (NEB names). Empty sweeps every exact site in the catalogue; degenerate (IUPAC) sites are never swept and are listed in the result. Advanced.
_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.
shufflesNoshuffles per record. The smallest p the null can give is 1 / (shuffles + 1). Only applies when null_model is shuffle.
positionsNoalso return every occurrence's coordinate and strand, for drawing on a map. Advanced.
sequencesNothe corpus: wire upstream records here, or list library document names. Each is one record of the sweep.
shuffle_kNok-let size the shuffles preserve. 2 holds every dinucleotide count, so CpG depletion is in the null. Only applies when null_model is shuffle.
correctionNoBenjamini–Hochberg (independent or positively dependent tests) or Benjamini–Yekutieli (any dependence, including overlapping sites on one sequence; less power). Values: BH (Benjamini–Hochberg), BY (Benjamini–Yekutieli).BH
null_modelNowhat a site is rare relative to: base composition, shuffles that keep the record's own k-let counts, or the record's own (k−1)-mers for a site of length k (R'MES maximal order, Gaussian score). Values: order0 (Base composition (order-0)), shuffle (k-let shuffles (uShuffle)), maximal (Maximal-order Markov (R'MES)).order0
include_subjectNoalso sweep the node's subject (the open or wired record).
min_site_lengthNoshortest recognition site swept. In bp.

Output Schema

TableJSON Schema
NameRequiredDescriptionDefault
sweepYesthe full result, including the records swept and the degenerate sites skipped
testsYesevery (record, site) test: observed, expected, ratio, both tails and both q-values; under the maximal-order null also the R'MES score z and its variance
avoidedYestests whose low-tail q is at or below alpha, most extreme ratio first
enrichedYestests whose high-tail q is at or below alpha, most extreme ratio first
nullModelYesthe null the p-values stand on, stated with every run
correctionYesthe multiple-testing correction that ran, and why it is valid
avoidedCountYeshow many (record, site) tests were called avoided
enrichedCountYeshow many (record, site) tests were called enriched

Schema Changelog

Changes observed during successful MCP inspections.

  1. First observed

TDQS

A4.6/5.0
Behavior5/5

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

Annotations only declare read-only, non-destructive, closed-world. The description adds what those annotations cannot: multiple-testing behavior (FDR controlled across every sequence/site pair, with BH vs BY choices described in the schema), that only exact sites are swept while degenerate IUPAC sites are excluded and merely listed, and the statistical meaning of each null. This is meaningful disclosure beyond structured hints.

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?

Purpose is front-loaded in the opening sentence and every subsequent sentence carries information: the biological interpretation, then the three nulls with their rationale. The sentences are long and clause-dense, and the null explanation partially overlaps the enum descriptions, but there is little outright filler.

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

Completeness5/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 formatting is already covered; all 13 parameters have full schema descriptions; and annotations cover the safety profile. The description adds the conceptual model (nulls, FDR, exact-vs-degenerate handling) an agent needs to call this statistical tool correctly, leaving no meaningful gap.

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

Parameters4/5

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

Schema coverage is 100%, so the baseline is 3, but the description goes further by explaining what the null_model alternatives actually mean statistically (order-0 vs k-let preserving vs maximal-order Markov) rather than just restating the enum labels. It does not add much for the other parameters, which the schema already documents.

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

Purpose5/5

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

The first sentence gives a specific verb (sweep) and resource (every exact recognition site in the NEB catalogue), plus the distinctive output (avoided/enriched sites vs a stated null, with FDR control). That combination is clearly separable from siblings like find_restriction_sites or enzyme_cut_counts, which merely report occurrences without a null model or significance testing.

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

Usage Guidelines4/5

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

The description supplies strong context for use — avoidance as evidence of restriction-modification systems, and rarity as a property a design can rely on — and explains why one would pick each null (shuffle keeps dinucleotide skews in the null rather than reporting them as findings; maximal follows published R'MES practice). It stops short of naming a sibling tool or stating explicit when-not conditions, which is the only gap.

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