注意事项
1 资源版本一致性
所有分析场景都应遵循以下原则:
- 参考基因组 FASTA、
.fai、.dict必须与下游分析使用的资源完全一致。 - dbSNP、known sites、interval/BED、germline resource、PoN 等资源不可混用不同参考版本。
- hg19 / GRCh37、hg38 / GRCh38 等不同版本资源混用,可能导致比对、重校正、变异检测和结果解释全部产生系统性偏差;请特别避免“版本混用”。
2 WGS / WES 分析
- 比对时需要参考序列索引文件,请先使用 aligner 的
--build-index构建。索引类型用--index-type选择,可选bwaFast(默认,标准索引,结果与 bwa 一致)、bwaUltra(高内存索引,更快、占用更多内存,结果与 bwa 一致)、minibwa(minibwa 索引)。比对时的--index-type必须与建索引时一致。 - bqsr 默认只输出 recal.table;如果需要碱基质量校正后的 BAM 文件,需要设置
--out-BAM。 - 对于 WES / panel 项目,variantCaller 阶段应显式设置 interval/BED;做 BAM 统计时也应显式使用同一 BED。
- 如果不在 WES / panel 项目中指定 BED,默认统计口径会退化为全基因组范围,导致覆盖度和深度等指标不准确。
3 联合变异检测
- 输入的 GVCF 文件需要带
.tbi索引。 - 参考基因组 FASTA 需要同时具备
.fai与.dict。 - 样本量很大时,应先小范围按染色体或区间测试资源占用,再决定 jobs 与 memory 的配置。
- 队列项目建议统一保留单样本 GVCF,而不是只保留单样本 VCF。
4 体细胞变异检测
- 优先使用 tumor-normal 模式;若缺少 matched normal,再使用 tumor-only。
- Mutect2 的原始输出通常不应直接作为最终交付结果。
- 推荐继续使用 FilterMutectCalls 等后过滤步骤,得到更适合解释的结果。
- tumor-only 场景更依赖 germline resource、PoN 及后过滤策略,解释时应更保守。
5 压缩与解压
- 压缩输入的 FASTQ 文件名后缀需为
*.fq、*.FASTQ、*.fq.gz或*.FASTQ.gz。 - 输入参考序列文件应为文本 FASTA,文件名后缀需为
*.fa或*.fasta。 - 解压输入的 arc 文件必须是本程序压缩得到的结果文件。
- 解压时使用的参考序列必须与压缩时保持一致,否则会出现
[ref load md5 error]。 - FASTQ 不完整时,可能出现
[read FASTQ error]。 - read 的 seq 和 qual 不等长时,可能出现
[seqlen isnot equal quallen]。 - 文件读写异常时,可能出现
[read file err]或[write file err]。 - 若数据块校验失败,可能出现
[block check fail],通常说明压缩过程异常或压缩文件未完整保存。 - 大规模压缩场景建议提前设置
ulimit -n。
6 license 相关
若工具 DCS Tools 或 lickit 无法连接 licproxy,常见报错包括:
- connection refused
- couldn't connect to server
建议排查顺序:
6.1 检查与 license 代理是否联通
curl -v http://$PROXY_HOST/api/licensor/licenses/anon/pool
6.2 在部署 licproxy 的机器上检查进程
ps aux | grep licproxy
若进程不存在,请先按 快速开始 中的说明运行代理程序。
6.3 检查 PROXY_HOST 是否设置正确
例如:
export PROXY_HOST=172.20.9.149:8909
6.4 常见认证错误提示
- Permission denied:无权限或验证失败
- Size must be greater than 0:提交线程数必须大于零
- Pool size exceeded:提交线程数大于当前线程数
- License expired:许可证已过期
- Unknown error:未知错误
6.5 认证环境变量优先级
认证相关环境变量优先级为:
DCS_DONGLE > DCS_DONGLE_SERVER > PROXY_HOST
若同时设置多个环境变量,程序将按上述优先级顺序选取一种认证模式。
离线单机模式:
export DCS_DONGLE=1
离线集群模式:
unset DCS_DONGLE
export DCS_DONGLE_SERVER=127.0.0.1:8910
在线认证模式:
unset DCS_DONGLE
unset DCS_DONGLE_SERVER
export PROXY_HOST=127.0.0.1:8909