我在使用dRep的时候出现了一些错误。
我使用的dRep的版本是:v3.4.5
在使用依赖项检测时的输出:只有centrifuge是error,别的都是all good
我原本要对38530个病毒基因组进行去冗余操作,但是在输出的文件夹dereplicated_genomes中发现存在35259个病毒基因组。以下是我的运行代码:
time dRep dereplicate drep/ -g /data/ldd/data_for_processing/fna_list.txt -p 128 -l 2000 --S_algorithm fastANI --ignoreGenomeQuality -pa 0.8 -sa 0.97 -nc 0.85 -comW 0 -conW 0 -strW 0 -N50W 0 -sizeW 1 -centW 0 --skip_plots
我检查了logger.log文件,我猜测问题可能出现在第二步。dRep dereplicate Step 2. Cluster
以下内容截取自logger.log
09-09 16:09 DEBUG running cluster 221
09-09 16:09 DEBUG /home/user/miniconda3/envs/checkm1/bin/fastANI --ql /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList --rl /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList -o /data/ldd/data_for_processing/drep/data/fastANI_files/fastANI_out_xoimckskpn --matrix -t 128 --minFraction 0 xoimckskpn
09-09 16:09 DEBUG running cluster 5566
09-09 16:09 DEBUG /home/user/miniconda3/envs/checkm1/bin/fastANI --ql /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList --rl /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList -o /data/ldd/data_for_processing/drep/data/fastANI_files/fastANI_out_imxwmalabg --matrix -t 128 --minFraction 0 imxwmalabg
09-09 16:09 ERROR CRITICAL ERROR WITH SECONDARY CLUSTERING CODE imxwmalabg; SKIPPING
09-09 16:09 ERROR CRITICAL ERROR WITH PRIMARY CLUSTER 5566; TRYING AGAIN
09-09 16:09 DEBUG /home/user/miniconda3/envs/checkm1/bin/fastANI --ql /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList --rl /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList -o /data/ldd/data_for_processing/drep/data/fastANI_files/fastANI_out_mhnylcorro --matrix -t 128 --minFraction 0 mhnylcorro
09-09 16:09 ERROR CRITICAL ERROR WITH SECONDARY CLUSTERING CODE mhnylcorro; SKIPPING
09-09 16:09 ERROR DOUBLE CRITICAL ERROR AGAIN WITH PRIMARY CLUSTER 5566; SKIPPING
09-09 16:09 DEBUG running cluster 11328
09-09 16:09 DEBUG /home/user/miniconda3/envs/checkm1/bin/fastANI --ql /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList --rl /data/ldd/data_for_processing/drep/data/fastANI_files/tmp/genomeList -o /data/ldd/data_for_processing/drep/data/fastANI_files/fastANI_out_fnuwccsket --matrix -t 128 --minFraction 0 fnuwccsket
因为不知道应该加什么标签,所以把Python、java和c++都加上去了