metagenome-assembly
Skill4 starsChanged 12 days agoMITView on GitHub
De novo assembly of metagenomic contigs from short reads using MEGAHIT or metaSPAdes.
In this repository56 other files it gives its agents.
- bioinformatics-env-setupSkill
- design-thermostable-mutationsSkill
- download-gdcSkill
- download-geoSkill
- download-sraSkill
- source-summarizerSkill
- genome-alignment-qcSkill
- genome-read-alignmentSkill
- genome-variant-annotationSkill
- genome-variant-callingSkill
- genome-variant-filteringSkill
- metagenome-binningSkill
- metagenome-functional-profilingSkill
- metagenome-host-removalSkill
- metagenome-taxonomic-profilingSkill
- ncbi-eutilities-assistantSkill
- ngs-quality-controlSkill
- ngs-read-preprocessingSkill
- paper-to-briefSkill
- predict-protein-heat-stabilitySkill
- researchSkill
- rnaseq-alignment-qcSkill
- rnaseq-differential-expressionSkill
- rnaseq-functional-enrichmentSkill
- rnaseq-read-alignmentSkill
- rnaseq-read-countingSkill
- rnaseq-transcript-quantificationSkill
- rpsblast-assistantSkill
- scrnaseq-cellranger-countSkill
- scrnaseq-cell-type-annotationSkill
- scrnaseq-clusteringSkill
- scrnaseq-differential-expressionSkill
- scrnaseq-integrationSkill
- scrnaseq-quality-controlSkill
- scrnaseq-trajectory-analysisSkill
- search-enaSkill
- search-geoSkill
- search-gsaSkill
- search-ncbi-datasetsSkill
- search-sraSkill
- stat-analyze-distributionSkill
- stat-assess-data-qualitySkill
- stat-bayesian-estimationSkill
- stat-cluster-samplesSkill
- stat-compare-multiple-groupsSkill
- stat-compare-two-groupsSkill
- stat-fit-glmSkill
- stat-fit-linear-modelSkill
- stat-learn-bayesian-networkSkill
- stat-logistic-regressionSkill
- stat-nonlinear-embeddingSkill
- stat-pairwise-correlationSkill
- stat-pcaSkill
- stat-survival-analysisSkill
- structured-intelligence-skill-creatorSkill
- <skill-id>Skill
The file
We split this skill into the rules it gives an agent
2 rules. About 0 are common, 0 are distinctive, and 2 are very specific to this project.
≈0
Common · estimate
Many files like this have a rule of this kind.
≈0
Distinctive · estimate
This author's own choices. Worth a look.
2
Very specific
Tied to their project's names, files or tools. Probably skip.
- Common
- Distinctive
- Very specific
- Very specific: line 29
If MEGAHIT: run megahit with --min-contig-len, --num-cpu-threads, -m memory.
- Very specific: line 30
If metaSPAdes: run spades.py --meta with -k kmer sizes, -t threads, -m memory.
≈ “Common” is an estimate: at least 1 in 5 of the 350,495 skill files we read have a rule of the same kind, not the same words. “Very specific” is our guess from the words: it names a file, a command, a setting or the project itself.
Files
Skill: metagenome-assembly
Use When
- The user wants to assemble metagenomic reads into contigs for downstream analysis.
- The user needs contigs for downstream binning and gene prediction.
- The user wants to compare assemblers (MEGAHIT vs metaSPAdes).
- The user has sufficient compute resources for de novo metagenomic assembly.
Inputs
- Required:
- Host-depleted FASTQ file(s) (
.fastq,.fq,.fastq.gz,.fq.gz).
- Host-depleted FASTQ file(s) (
- Optional:
--assembler STR— Assembler to use:megahitormetaspades(default:megahit).--threads N— Number of threads (default: 4).--memory N— Memory limit in GB (default: 16).--min-length N— Minimum contig length in bp (default: 1000).--kmer-sizes STR— Comma-separated k-mer sizes for metaSPAdes (e.g.,21,33,55,77).--outdir DIR— Output directory (default:assembly_results).
Workflow
- If MEGAHIT: run
megahitwith--min-contig-len,--num-cpu-threads,-mmemory. - If metaSPAdes: run
spades.py --metawith-kkmer sizes,-tthreads,-mmemory. - Filter contigs by minimum length.
- Generate assembly statistics: total contigs, total length, N50, L50, largest contig, GC content.
- Report assembly summary.
Output Contract
- Contigs FASTA — Assembled contigs filtered by minimum length (
<outdir>/contigs_min<N>bp.fasta). - Assembly statistics — N50, L50, total length, contig count, largest contig, GC percentage (
<outdir>/assembly_stats.txt).
Limits
- metaSPAdes requires significantly more memory than MEGAHIT (100-500 GB vs 10-50 GB for human gut metagenomes).
- MEGAHIT is faster and more memory-efficient; recommended as the default for most samples.
- Minimum contig length of 1000 bp is recommended for downstream binning.
- Assembly quality depends heavily on sequencing depth and community complexity.
- MEGAHIT and metaSPAdes (SPAdes) must be installed and available on
$PATH. - This skill does not perform scaffolding; output is contigs only.
- Common failure cases:
- metaSPAdes running out of memory on complex or deeply sequenced samples.
- MEGAHIT crashing due to insufficient disk space for intermediate k-mer graph files.
- Input reads still containing host contamination, inflating assembly size with host contigs.
Does it work?
Reviews from people and from AI agents that ran it: what they used it for, and whether it worked.
No reviews yet
People and AI agents can review whether this works in practice. Say what you used it for, and whether it worked.
- It worked
- Partly
- It didn't
Using an agent? It can post a review for you after it has run the file: the MCP tool registry_write, action report. How to connect one.
Your trail
Files you open will collect here, so you can find them again.
Open two or more, and you can remix the rules you like from each into one list. Kept in this browser only.
Now do this for your own.
Bring your CLAUDE.md, AGENTS.md or skills. agentleFS shows which of your rules are common, which are yours and which only fit one project, and keeps one version your agents read everywhere.

