Hello! I am currently working on the DADA2 step in the qiime2 workflow. My data was already demultiplexed prior to qiime2 analysis (data was retrieved from SRA repository, so the details i have about the data is info provided on SRA as well as the paper). i successfully loaded the data as an artifact and was able to summarize it to look at the average quality over the reads to instruct my dada2 parameters. I then input the data with the following parameters into DADA2
qiime dada2 denoise-paired --i-demultiplexed-seqs pe-demux.qza --p-trunc-len-f 20 --p-trunc-len-r 20 --o-table table.qza --o-representative-sequences rep-seqs.qza --o-denoising-stats denoising-stats.qza --verbose
I have tried to run it a few times and am consistently getting this error:
Input must be a valid sequence table.
Calls: removeBimeraDenovo -> isBimeraDenovoTable
Execution halted
Traceback (most recent call last):
File "/data/apps/extern/anaconda/envs/qiime2/2022.2/lib/python3.8/site-packages/q2_dada2/_denoise.py", line 279, in denoise_paired
run_commands([cmd])
File "/data/apps/extern/anaconda/envs/qiime2/2022.2/lib/python3.8/site-packages/q2_dada2/_denoise.py", line 36, in run_commands
subprocess.run(cmd, check=True)
File "/data/apps/extern/anaconda/envs/qiime2/2022.2/lib/python3.8/subprocess.py", line 516, in run
raise CalledProcessError(retcode, process.args,
subprocess.CalledProcessError: Command '['run_dada_paired.R', '/tmp/tmplppuhq62/forward', '/tmp/tmplppuhq62/reverse', '/tmp/tmplppuhq62/output.tsv.biom', '/t
mp/tmplppuhq62/track.tsv', '/tmp/tmplppuhq62/filt_f', '/tmp/tmplppuhq62/filt_r', '20', '20', '0', '0', '2.0', '2.0', '2', '12', 'independent', 'consensus', '
1.0', '1', '1000000']' returned non-zero exit status 1.
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/data/apps/extern/anaconda/envs/qiime2/2022.2/lib/python3.8/site-packages/q2cli/commands.py", line 339, in call
results = action(**arguments)
File "", line 2, in denoise_paired
File "/data/apps/extern/anaconda/envs/qiime2/2022.2/lib/python3.8/site-packages/qiime2/sdk/action.py", line 245, in bound_callable
outputs = self.callable_executor(scope, callable_args,
File "/data/apps/extern/anaconda/envs/qiime2/2022.2/lib/python3.8/site-packages/qiime2/sdk/action.py", line 391, in callable_executor
output_views = self._callable(**view_args)
File "/data/apps/extern/anaconda/envs/qiime2/2022.2/lib/python3.8/site-packages/q2_dada2/_denoise.py", line 292, in denoise_paired
raise Exception("An error was encountered while running DADA2"
Exception: An error was encountered while running DADA2 in R (return code 1), please inspect stdout and stderr to learn more.
Plugin error from dada2:
An error was encountered while running DADA2 in R (return code 1), please inspect stdout and stderr to learn more.
See above for debug info.
Running external command line application(s). This may print messages to stdout and/or stderr.
The command(s) being run are below. These commands cannot be manually re-run as they will depend on temporary files that no longer exist.
Command: run_dada_paired.R /tmp/tmplppuhq62/forward /tmp/tmplppuhq62/reverse /tmp/tmplppuhq62/output.tsv.biom /tmp/tmplppuhq62/track.tsv /tmp/tmplppuhq62/filt_f /tmp/tmplppuhq62/filt_r 20 20 0 0 2.0 2.0 2 12 independent consensus 1.0 1 1000000
I found a forum question that pointed to the number of threads being used as a possible explanation, but i have reduced the threads from 24 to 6, and am still getting the same error.
I am running qiime2 on an HPC with Qiime natively installed (not in my own personal conda environment).
Thanks for helping!