Skip to content

Instantly share code, notes, and snippets.

@sinamajidian
Last active July 20, 2025 20:43
Show Gist options
  • Save sinamajidian/022066d28b2dcb5b0d2d112a1f987d17 to your computer and use it in GitHub Desktop.
Save sinamajidian/022066d28b2dcb5b0d2d112a1f987d17 to your computer and use it in GitHub Desktop.
How we ran metabuli for the Movi color study

In our Movi color study, we ran metabuli using the following script. We used two read datasets from CAMI (long and short read) for metagenomic classification.

Install

conda create -n metab python=3.12
conda activate metab

wget https://mmseqs.com/metabuli/metabuli-linux-avx2.tar.gz

Build

metabuli/bin/metabuli  build metabuli_db fasta.list accession2taxid --taxonomy-path . 

The file fasta.list includes 75166 fasta files (e.g. NC_000913.3.fa, NC_000914.2.fa, ).

classify

metabuli classify --seq-mode $mode cami_plant/reads/anonymous_reads_0.fq metabuli_db outdir_reads0_1thread jobid  --threads 1

mode is set to 1 for short read and 3 for long read.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment