近来想用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的参数提示，修改后重新运行，结果顺利输出