Skip to main content
Glama

discover_site_avoidance

Read-only

Sweep a set of sequences against every exact recognition site in the NEB catalog 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

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 catalog; 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

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. Changed1 schema field changed
    • changedInput schema / properties / enzymes / description
      Previous value: -"restrict 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."New value: +"restrict the sweep to these enzymes' sites (NEB names). Empty sweeps every exact site in the catalog; degenerate (IUPAC) sites are never swept and are listed in the result. Advanced."
  2. First observed

TDQS

A4.1/5.0
Behavior4/5

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

Annotations already declare readOnlyHint=true, destructiveHint=false, openWorldHint=false, so the safety profile is covered. The description goes beyond that by disclosing the statistical behavior: FDR is controlled 'across every (sequence, site) test,' and the three nulls are explained, including what each does and does not treat as signal (e.g. uShuffle holding CpG depletion in the null). It omits compute cost/runtime characteristics, keeping it out of the top band.

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 action is front-loaded in the first sentence with the enrichment/avoidance outcome and the FDR guarantee. The remaining sentences are dense but each carries information; it runs long and would benefit from breaking the null-model discussion out, so it is not maximally tight.

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, the description need not explain return values, and annotations carry the safety profile. It covers purpose, statistical meaning, and the three null models, which is enough for correct invocation; the main gap is the absence of explicit when-to-use routing against restriction-site siblings.

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 description coverage is 100%, so the baseline would be 3. The description adds genuine meaning on top, elaborating the null_model choices (order-0 Poisson, uShuffle k-let shuffles, R'MES maximal-order Markov judging a length-k site against the record's own (k−1)-mers) far beyond the enum labels, which is the semantically dominant parameter.

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?

States a specific verb and resource with precise scope: 'Sweep a set of sequences against every exact recognition site in the NEB catalog and report the sites each one carries far less (avoided) or far more (enriched) often than a stated null predicts.' The statistical framing (FDR-controlled avoidance/enrichment) clearly separates it from simple finders like find_restriction_sites or enzyme_cut_counts, so an agent can tell what it does without opening the schema.

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?

It gives motivating context ('Avoidance is evidence of the restriction–modification systems a lineage has lived with, and a site a design can rely on being rare') but never states when to choose this over alternatives or any exclusions. Usage is implied by the biological framing rather than spelled out, and no sibling is named as an alternative.

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