项目场景:
提示:这里简述项目相关背景:
使用snippy4.6.0版本call snps,
问题描述
提示:这里描述项目中遇到的问题:
手动下载的安装包,参考序列使用fasta文件时可以运行,使用gbk文件时出现报错
File %public/homeb20213050430/wyf/sxm/gene NDM/Result/169 extract snippy/reference/.refgenesgff line 2622CDS 1478615076snippyWARNING TRANSCRIPT NOT FOUND:CAMt fd TCT TRANSCRIPT JFONICPM 006. C TSC TRANSCRIPT JFONICPM 026 a Cee 'GENE FOcodon start : 1
gene : VirB4id : JFONICPM 00026
inference : ab initio prediction:Prodigal:002006,similar to AA sequence:UniProtKB:O9RPY1locus tag : FONICPM 00026name : virB4
product : Type IV secretion system protein virB4
source : snippytransl table : 1
PASIAO ·:MSVEKGL IRPALIRCLGVPL YPFL GMCVCVLIGVWHEAMYALIPGWYAKRVTKIDERFDIYLRMOKGNPL ANKRENAVHYAGSSYDAVDISK VDNEM
type : CDSFile%publichome20213050430/wyf'sxm/gene NDMMResult/169 extract sippyreferencel/refigenes gff line 27 2 1 Sippy CDS15101WARNING TRANSCRIPT NOT FOUND: TOO may "WARNNG RANSCRIPT NOT FOUND' wmngs. no futher waings wil be shown.
WARNING FRAMES ZERO: All frames are zero! This seems rather odd. please check that 'frame' iformation in vour 'genes' file is accurate.
ERROR: CDs check file /xxx/xx/reference/./ref'cds.fa' not found.
ERROR: Protein check file %ublichome b20213050430wvf/sxm/gene NDM/Result/169 extract snippv/reference//ref/orotein.fa' not found
ERROR: Database check failed.
原因分析:
提示:这里填写问题的分析:
例如:查询后提示snippy和snpeff版本可能存在冲突,但根据snippy官方文件显示,版本符合要求,snippy版本为4.6.0,snpeff版本为5.1d。
进一步查询后发现,snpeff构建过程出现问题,此问题常出现于最新的5.1版本。
解决方案:
提示:这里填写该问题的具体解决方案:
参考其他人的报错解决方案(https://groups.google.com/g/linux.debian.bugs.dist/c/660LIVsRuVc)后,解决方法如下:
1、使用conda创建新的虚拟环境,安装snippy=4.6.0和snpeff=5.0,为提升安装速度,首先安装mamba
##创建虚拟环境
conda create -n snippy
##激活虚拟环境
conda activate snippy
conda install -c conda-forge mamba
##使用mamba安装snippy和snpeff的特定版本
mamba install -y -c conda-forge -c bioconda -c defaults snpeff=5.0 snippy=4.6.0
##检查是否安装成功,版本是否正确