跳转至

CLI 参考

gside species

物种鉴定主命令。

用法

gside species <contigs.fasta> [选项]

参数

参数 类型 默认 说明
contigs 位置参数 无 组装后的 contigs FASTA 文件
-m,--mode 选项 marker 鉴定模式(见下)
-D,--db-dir 选项 自动 数据库根目录
-f,--format 选项 json 输出格式:json / tsv / md
-n,--no-header 选项 关 隐藏 TSV 表头行(仅 tsv;md/json 保持结构)
-a,--auto-setup 选项 关 鉴定前自动安装缺失的数据库 tier(进度走 stderr)
-j,--jobs N 选项 1 并发鉴定的输入文件数(输出顺序恒为输入顺序)
-t,--threads N 选项 4 每次比对调用的线程数(blastn -num_threads / skani -t / mash -p)
-o,--output PATH 选项 stdout 报告写入 PATH 而非 stdout(截断写入;tsv/md 逐输入流式落盘、逐行 flush——批次中途失败已写前缀仍保留;json 为单文档收尾一次写入;stdout 不再输出数据,提示信息走 stderr)

鉴定模式

共 6 种模式,默认 marker:

Mode 说明 数据库需求
marker 靶基因组合规则(blastn) mini(内置)
panel ANI 比对精选面板(skani) panel tier
mash_refseq MinHash 距离(mash) mash tier
sourmash sourmash GTDB gather sourmash GTDB 库
gtdbtk GTDB-Tk 系统发育 placement(classify_wf) 外部安装的 GTDB-Tk(GTDBTK_DATA_PATH)
all 运行全部 5 种方法 + 仲裁 按需

注意 gtdbtk 较慢(每基因组分钟级);all 模式也会运行它——如需跳过请显式指定 -m。

结构化错误字段(面向 agent)

当某方法因数据库 tier 缺失而失败时,该方法的 JSON 条目在 error 字符串 之外还携带结构化字段,无需正则解析:

{
  "method": "sourmash",
  "error": "sourmash GTDB database not found in ... (run: gside db setup --tier sourmash)",
  "suggested_command": "gside db setup --tier sourmash",
  "missing_tier": "sourmash",
  "estimated_size": "~3.9GB"
}

加上 --auto-setup 后,gside 会自行执行 suggested_command(进度走 stderr),对该方法重试一次,并重算 verdict。

示例

# 靶基因鉴定(最快,无需额外数据库)
gside species contigs.fna --mode marker

# ANI 鉴定(精确)
gside species contigs.fna --mode panel

# sourmash GTDB gather
gside species contigs.fna --mode sourmash

# GTDB-Tk 系统发育 placement(慢,需外部 gtdbtk + GTDBTK_DATA_PATH)
gside species contigs.fna --mode gtdbtk

# 多法仲裁(marker + panel + mash_refseq + sourmash + gtdbtk)
gside species contigs.fna --mode all

# 指定数据库目录
gside species contigs.fna --mode panel --db-dir /path/to/dbs

# 自动安装缺失数据库 tier 后鉴定
gside species contigs.fna --mode all --auto-setup

# 报告写入文件(stdout 保持为空,提示信息走 stderr)
gside species contigs.fna -m all -o report.json

gside db

数据库管理子命令。

子命令

子命令 说明
status 显示数据库就绪状态
list 列出可用 tier
setup 安装/更新数据库

status 与 list 支持 -f/--format [table|tsv|json](默认 table)。 tsv 与 json 为机器可读格式(如 db status --format tsv 输出 COMPONENT\tREADY\tTIER\tPATH 表,READY 列为 ready/missing)。

setup 选项

选项 默认 说明
-t,--tier panel mini / panel / mash / sourmash / meta / meta_sylph / all
-s,--source 无 从本地已有数据目录复制
-f,--format text 输出格式:text / json(json 模式下进度信息走 stderr,结果字典以单个 JSON 文档输出到 stdout;退出码契约不变)

all 仅安装 panel + mash(约 310MB)。sourmash GTDB 库(约 3.9GB)不含 在 all 内,需单独安装:

gside db setup --tier sourmash

宏基因组 tier meta(Kraken2 索引,约 5.5GB)与 meta_sylph(sylph GTDB 库,约 14GB)需手动安装——db setup 只打印上游下载指引而不下载。 放置方法见宏分类学分析。

示例

gside db status
gside db status --format tsv    # 机器可读
gside db setup --tier all
gside db setup --tier sourmash   # 不含在 --tier all 内(~3.9GB)
gside db setup --tier panel --source /path/to/existing/db

gside validate

标记基因组装验证,可选叠加 CheckM2 + GTDB-Tk。

用法

gside validate <contigs.fasta> [选项]

参数

参数 类型 默认 说明
contigs 位置参数 — 组装后的 contigs FASTA 文件
--mode 选项 simple 验证深度:simple(标记基因)或 standard(叠加 CheckM2 + GTDB-Tk)
--output-dir 选项 自动 validation.json 输出目录

standard 模式优雅降级:缺 CheckM2/GTDB-Tk 数据库或二进制时只告警并返回 部分结果,不崩溃。通过 CHECKM2DB / GTDBTK_DATA_PATH(或 GTDBDB) 环境变量配置。

示例

gside validate contigs.fna --mode simple
gside validate contigs.fna --mode standard --output-dir ./taxonomy

gside meta

宏基因组物种分析(Kraken2+Bracken 与 sylph),输入 FASTQ 读段或组装 contigs。详见宏分类学分析。

用法

gside meta <reads.fastq> [reads_2.fastq] [选项]

参数

参数 类型 默认 说明
reads 位置参数 — 1 个文件 = 单端读段或 contigs FASTA;2 个文件 = 双端 FASTQ
-m,--mode 选项 all 分析方法:kraken2 / sylph / all(all 合并两份物种列表)
-D,--db-dir 选项 自动 数据库根目录($GSIDE_DB_DIR 或 data/db)
-f,--format 选项 json 输出格式:json / tsv
-o,--output PATH 选项 stdout 报告写入 PATH 而非 stdout(json 为单文档收尾一次写入;stdout 不再输出数据,提示信息走 stderr)
-t,--threads N 选项 4 每次工具调用的线程数

示例

gside meta sample_1.fastq sample_2.fastq
gside meta sample_1.fastq sample_2.fastq -m kraken2 -f tsv -o profile.tsv
gside meta contigs.fna          # 仅 kraken2(sylph 只接受读段)
gside meta reads.fastq -m sylph

gside --version

gside --version
# gside, version 0.1.4

gside --version --json     # 或 -V -J
# {"tool": "gside", "version": "0.1.4"}