genomics-assembly

name	genomics-assembly
description	Genome assembly quality assessment: N50/N90/L50/L90 (QUAST-compatible), GC content, contig length distribution, completeness estimation. Wraps SPAdes, Megahit, Flye, Canu.
version	0.2.0
author	OmicsClaw
license	MIT
tags	["genomics","assembly","SPAdes","Flye","Canu"]
metadata	{"omicsclaw":{"domain":"genomics","emoji":"🧬","trigger_keywords":["genome assembly","de novo","SPAdes","Megahit","Flye","Canu"],"allowed_extra_flags":["--method"],"legacy_aliases":["assemble"],"saves_h5ad":false}}

🧬 Genome Assembly

De novo genome assembly for short and long reads. Wraps SPAdes, Megahit, Flye, and Canu.

CLI Reference

python omicsclaw.py run genomics-assembly --demo
python omicsclaw.py run genomics-assembly --input <reads.fastq> --output <dir>

Why This Exists

• Without it: Assemblies require intense memory management and parameter orchestration per graph build
• With it: Automated contig building and K-mer tuning logic across read modalities
• Why OmicsClaw: Unified containerized or local graph assembler invocation

Workflow

1. Calculate: Prepare k-mer frequencies or long-read overlaps.
2. Execute: Build de Bruijn or string graphs.
3. Assess: Perform contig polishing and scaffolding.
4. Generate: Output structural FASTA representations.
5. Report: Synthesize N50 stats and completeness metrics.

Example Queries

• "Assemble my isolate using SPAdes"
• "De novo genome assembly using Flye"

Output Structure

output_directory/
├── report.md
├── result.json
├── assembled.fa
├── figures/
│   └── assembly_graph.png
├── tables/
│   └── quast_metrics.csv
└── reproducibility/
    ├── commands.sh
    ├── requirements.txt
    └── checksums.sha256

Safety

• Local-first: Strict offline processing without external upload.
• Disclaimer: Requires OmicsClaw reporting structures and disclaimers.
• Audit trail: Hyperparameters and operational flow states are logged fully.

Integration with Orchestrator

Trigger conditions:

• Automatically invoked dynamically based on tool metadata and user intent matching.

Chaining partners:

• genomics-qc — Upstream read trimming
• annotation — Downstream genome annotation

Citations

• SPAdes
• Flye