Hello!
I encountered some errors when using dRep.
The version of dRep that I used was:
Check_dependencies:
I originally intended to perform redundancy removal on 38,530 viral genomes, but upon checking the output folder "dereplicated_genomes", I found that there were 35,259 viral genomes present. Here is the code I used:
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
I checked the logger.log file and I suspect the problem might lie in the second step:dRep dereplicate Step 2. Cluster
The following content is excerpted from 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
But I don't know how to solve this problem.
Thank you!
Hello!
I encountered some errors when using dRep.
The version of dRep that I used was:
Check_dependencies:
I originally intended to perform redundancy removal on 38,530 viral genomes, but upon checking the output folder "dereplicated_genomes", I found that there were 35,259 viral genomes present. Here is the code I used:
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
I checked the logger.log file and I suspect the problem might lie in the second step:dRep dereplicate Step 2. Cluster
The following content is excerpted from 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
But I don't know how to solve this problem.
Thank you!