某某电气

密封防水壳体成套配件解决方案提供商

全国免费服务热线

400-123-4567

当前位置:主页»新闻动态»常见问题解答»

科学网DeepARG的short

文章出处:威尼斯人平台 人气:发表时间:2019-12-19 17:36

近来想用DeepARG倒腾抗性基因的事,由于数据是测序的PE reads。不过安装的DeepARG版本提示输入文件为fasta格式或比对格式。

这两天看到其官网主页上提供了PE reads的网页版分析方案,在其代码放置主页上看到流程更新,https://bitbucket.org/gusphdproj/deeparg-ss/src/master/


第一:下载和安装

参考官网提供的介绍

https://bench.cs.vt.edu/deeparg


第二:尝试运行

运行脚本如下,其中数据是软件自带

export PATH="*/Python-2.7.14/bin:$PATH" export PYTHONPATH="*/Python-2.7.14/lib/python2.7/site-packages/:$PYTHONPATH" python */deeparg-ss/short_reads_pipeline/short_reads_pipeline.py --forward_pe_file ./reads/F.fq.gz --reverse_pe_file ./reads/R.fq.gz  --path_to_executables */short_reads_pipeline/bin/ --path_to_deeparg */deeparg-ss/ --output_file temp

前3步正常运行,第4步出现报错,报错信息如下:

Step 4: Quantification of ARG-like counts global name 'path_to_deeparg' is not defined

排查代码,在pairedEndPipelineClass.py中发现


调用了quantificationClass.py中的merge函数

查看其中的代码发现此间调用了deepargClass.py的dsize函数

def dsize(): return {i.split()[0].split("|")[-1].upper(): i.split()  for i in open(path_to_deeparg+'/database/v2/features.gene.length')}

将其中path_to_deeparg改成服务器安装DeepARG的路径

重新运行,第4步正常运行


第三:再次运行

运行脚本照旧,此时第1、2、3、4步正常运行

第5步出现报错提示dataset的路径出现错误,第5步调用了d16spipelineClass.py,查看源码,发现其中dataset的路径是

class d16sPipe():     def __init__(self, ggdata="/data/short_reads_pipeline/bin/gg13/dataset"):

于是将ggdata的路径设置为gg13的在服务器上的路径

重新运行,屏幕显示有很多乱码符号,查看报错原因是bam二进制文件打开出现问题,在d16spipelineClass.py查看源码发现

os.system(" ".join([     "samtools sort",     fi+".bam",     " -o", fi+".sorted.bam"     ]))

由于samtools版本的更新,此处运行不符合新的参数,

按照samtools的参数提示,修改后重新运行,结果顺利输出

2.png

此文关键字:DeepARG的short_reads_pipeline使用

同类文章排行

最新资讯文章

全国免费服务热线:13800000000 联系电话:400-123-4567

传真:+86-123-4567 地址:广东省广州市天河区88号

威尼斯人平台 版权所有 Power by DeDe58粤547800223-d3号

织梦模板