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.
conda create -n metab python=3.12
conda activate metab
wget https://mmseqs.com/metabuli/metabuli-linux-avx2.tar.gz
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, ).
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.