Snippy 怎么用？一份从安装到多样本比对的完整实操笔记

Snippy 怎么用一份从安装到多样本比对的完整实操笔记【免费下载链接】snippy:scissors: :zap: Rapid haploid variant calling and core genome alignment项目地址: https://gitcode.com/gh_mirrors/sn/snippySnippy 是一套面向单倍体基因组的快速变异检测流水线它能在参考序列与测序数据之间一次性定位 SNP 和插入缺失并自动完成注释与核心基因组比对。如果你是做细菌、病毒或质粒测序的研究人员下面这份笔记会按选型 → 安装 → 首跑 → 批量的顺序带你走完从零到产出变异清单的全过程。它到底帮你做什么设想一个常见场景你手上有十几个样本的测序数据想快速知道它们相对某个参考基因组各自发生了哪些突变还想比较它们之间的亲缘关系。手工比对、逐位筛查显然不现实Snippy 就是把这条流水线压缩成一条命令的工具。适用对象细菌、病毒、质粒、线粒体等单倍体基因组输出内容变异位点清单VCF/TAB、比对文件BAM、修正版基因组consensus、核心 SNP 比对core.aln工作方式内部依次调用 bwa 做比对、freebayes 找变异、snpEff 做注释对外只暴露一个入口输入形态适用情况双端 FASTQ最常见的 Illumina 测序产物单端 FASTQIon Torrent 等平台或双端文件只剩其一contigs原始 reads 已丢失仅剩拼装结果⚠️ 注意Snippy 面向单倍体样本。人类这类二倍体需要区分杂合位点不在它的处理范围内。开工前先核对这几样✅ 参考基因组FASTA 或 GenBank 格式均可多 contig 不影响✅ 测序数据FASTQ 或 FASTA支持 gzip 压缩✅ 一个专门存放结果的文件夹✅ 硬件多核 CPU 能显著提速官方在 64 核机器上验证过✅ 外部依赖bwa、minimap2、samtools、bcftools、bedtools、freebayes、snpEff、samclip、seqtk 等最后一项最容易被忽视Snippy 本体只是一个调度脚本真正干重活的是背后那一串外部程序。任何一个缺失流程都会在中途停下。所以判断安装是否成功不能只看 Snippy 本身。三条安装路径怎么取舍路径适合谁优点要注意什么Conda大多数用户依赖自动装齐、环境隔离、卸载干净需要先配置好 Bioconda 源HomebrewmacOS 用户一条命令搞定与系统软件统一管理依赖按 brew 方式处理源码想追最新版的人代码最新、便于二次开发依赖全部自己装PATH 需手动配置Conda 路线一条命令把依赖一起装齐conda install -c conda-forge -c bioconda -c defaults snippyHomebrew 路线brew install brewsci/bio/snippy源码路线先把代码拿到手再补依赖git clone https://gitcode.com/gh_mirrors/sn/snippy.git export PATH$PWD/snippy/bin:$PATH源码方式只是拿到 Snippy 本体配套程序仍需另行安装例如conda install -c bioconda bwa samtools bcftools freebayes snpeff samclip seqtk连续实操装完 → 验证 → 跑第一个案例第 1 步确认版本号。snippy --version预期会打印出版本字符串当前仓库为 5.0.0-dev如果提示找不到命令多半是 PATH 没配对回头检查一下路径配置。第 2 步检查全部依赖。snippy --check命令会逐个探测比对、排序、变异识别、注释等环节用到的外部程序可用的显示 OK 或版本信息。若有缺失就针对性补装然后重新检查直到全部通过。这一步值得认真做——检查通过再上真实数据能避免跑到一半才发现缺工具。第 3 步准备一份练手数据。仓库的 test 目录里已经备好 example.fna参考序列、example.gbk带注释版本和 example.bed区域文件。真实 reads 文件比较大先用 wgsim 按参考序列模拟一对双端 readscd test wgsim -S 7 -h -r 0.005 -N 12000 -1 100 -2 100 -d 200 example.fna sim_1.fq sim_2.fq这行命令的作用以参考序列为模板随机撒入约 0.5% 的变异生成 1.2 万对、长度 100 bp 的虚拟双端 reads用来验证整条流程是否通畅。第 4 步正式检出变异。snippy --cpus 4 --outdir demo1 --ref example.fna --R1 sim_1.fq --R2 sim_2.fq运行过程中会依次看到比对、变异识别等阶段日志收尾时打印结果目录路径并显示 Done。第 5 步查看产出。ls demo1此时文件夹里已经躺着整套标准产物。结果文件夹里都有什么文件作用snps.vcf注释后的标准变异文件snps.tab / snps.csv便于阅读的表格汇总snps.bam比对结果含未比对与多比对 readssnps.consensus.fa把变异回贴到参考序列的修正版基因组snps.raw.vcf / snps.filt.vcf过滤前后的变异记录snps.html可在浏览器打开的网页版汇总想直接看变异清单执行head -5 demo1/snps.tab表格各列的含义依次是变异所在的序列名CHROM、位置POS、类型TYPE取值 snp/mnp/ins/del/complex、参考碱基REF、样本碱基ALT、支持各碱基的 reads 计数EVIDENCE。如果参考用的是 example.gbk 这类带注释的文件表格还会多出基因名、产物描述和影响效应列——变异落在哪个基因、是否改变氨基酸一眼就能看到这是 Snippy 很贴心的设计。从单样本扩展到批量分析样本一多逐个手敲命令就不现实了。Snippy 提供了批量入口先准备一个制表符分隔的清单文件 input.tab每行描述一个样本SampleA /path/to/A_1.fq.gz /path/to/A_2.fq.gz SampleB /path/to/B.fq.gz SampleC /path/to/C.contigs.fa生成并检查批量脚本snippy-multi input.tab --ref Reference.gbk --cpus 16 runme.sh less runme.sh确认脚本内容无误后放行sh runme.sh脚本会逐个样本输出结果文件夹并在末尾自动调用 snippy-core把所有样本都有覆盖的位点聚合成核心 SNP 比对产出 core.aln多序列比对和 core.vcf多样本 VCF。core.aln 可以直接交给 FastTree 等建树工具绘制系统发育树做亲缘关系分析。三个高频调参场景场景一深度太高跑得特别慢。有的样本深度高达上千倍而多数变异在 50~100 倍深度下就能可靠检出。用 --subsample 按比例随机抽读snippy --subsample 0.1 --outdir demo2 --ref ref.fna --R1 R1.fq.gz --R2 R2.fq.gz场景二只关心特定区域。比如只想筛耐药基因上的突变把目标区域写进 BED 文件用 --targets 限定范围能省下大量计算snippy --targets sites.bed --outdir demo3 --ref ref.fna --R1 R1.fq.gz --R2 R2.fq.gz场景三只有 contigs 没有原始 reads。改用 --ctgs 传入 contigs 文件Snippy 会把它撕成约 250 bp 的伪 reads 再比对产物与 reads 样本完全兼容可以混在一起参与批量分析snippy --outdir demo4 --ref ref.fna --ctgs sample.fasta另外三个核心过滤参数也值得记住--mincov位点最少覆盖数默认 10、--minfrac支持变异碱基的最低比例、--minqual最低质量值默认 100它们共同决定了最终哪些变异会被保留。出问题时的排查速查表症状常见原因处理办法提示 command not foundPATH 没配好或依赖缺失用 which 逐个定位缺失工具补进 PATH 或补装运行极慢数据深度过高用 --subsample 按比例抽读--check 报某项缺失外部依赖未安装用 conda 补装对应包再重跑检查结果缺注释列参考文件是 FASTA 而非 GenBank换成带注释的 .gbk 文件重跑收尾跑通之后继续做什么用测试数据完整走一遍验证 → 单样本 → 批量的流程把每个环节的输出记在脑子里再碰真实数据。正式样本开跑前备份好原始 reads并为每个样本建立清晰的命名规范方便日后追溯。记录 --version 的输出和关键参数——变异检测结果与版本强相关注明版本能让你的结果更可信、可复现。进阶方向结合 --report 生成逐位点可视化报告或把 core.aln 送入建树、重组过滤流程做群体分析。Snippy 的价值在于把比对—变异识别—注释—汇总这条长流水线封装成一条命令。当你把第一份 reads 顺利变成一张清晰的变异表格时后续的群体分析和系统发育研究也就有了可靠的数据地基。【免费下载链接】snippy:scissors: :zap: Rapid haploid variant calling and core genome alignment项目地址: https://gitcode.com/gh_mirrors/sn/snippy创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考