Metagenome Binning

Workflow

1. Calculate contig coverage depth with MetaBAT2's jgi_summarize_bam_contig_depths.
2. Run MetaBAT2 for initial binning.
3. Optionally run additional binners (MaxBin2, CONCOCT) for DAS Tool refinement.
4. If multiple binners used: run DAS Tool to select the best non-redundant bin set.
5. Run CheckM2 predict on the final bin set for completeness and contamination estimates.
6. Classify bins by MIMAG standards: high-quality (>=90% complete, <5% contamination), medium-quality (>=50%, <10%), low-quality (below thresholds).
7. Report: number of bins, quality distribution, total high/medium/low quality MAGs.

Output Contract

• MAG FASTA files — One FASTA file per bin (<outdir>/bins/<bin_id>.fa).
• CheckM2 quality report — Completeness and contamination per bin (<outdir>/checkm2/quality_report.tsv).
• Bin quality summary — Per-bin completeness, contamination, N50, genome size, and MIMAG classification (<outdir>/bin_summary.txt).

Limits

• MetaBAT2 and CheckM2 must be installed and available on $PATH.
• CheckM2 requires a pre-downloaded database (approximately 3 GB).
• DAS Tool requires results from multiple binners as input; it is optional but recommended.
• Assembly must have sufficient depth for binning (>10x per genome is recommended).
• Contigs shorter than 1500 bp are typically excluded from binning.
• MIMAG quality standards: high-quality draft (>=90% completeness, <5% contamination), medium-quality draft (>=50% completeness, <10% contamination).
• Common failure cases:
  • BAM file not sorted or indexed, causing jgi_summarize_bam_contig_depths to fail.
  • CheckM2 database not downloaded or path not set, preventing quality assessment.
  • Insufficient sequencing depth (<10x per genome) producing fragmented or empty bins.