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.
R version 4.1.0 (2021-05-18)
Loading required package: Rcpp
DADA2: 1.22.0 / Rcpp: 1.0.8.3 / RcppParallel: 5.1.5
Filtering ..................
Learning Error Rates
266474 total bases in 1238 reads from 18 samples will be used for learning the error rates.
234067 total bases in 1079 reads from 18 samples will be used for learning the error rates.
Denoise samples ..................
Error in FUN(X[[i]], ...) :
Non-corresponding derep-class and dada-class objects.
Calls: mergePairs -> lapply -> FUN
Execution halted
Traceback (most recent call last):
File "/public1/home/sca0719/anaconda3/envs/qiime2-2022.2/lib/python3.8/site-packages/q2_dada2/_denoise.py", line 279, in denoise_paired
run_commands([cmd])
File "/public1/home/sca0719/anaconda3/envs/qiime2-2022.2/lib/python3.8/site-packages/q2_dada2/_denoise.py", line 36, in run_commands
subprocess.run(cmd, check=True)
File "/public1/home/sca0719/anaconda3/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/tmp7lf5_feg/forward', '/tmp/tmp7lf5_feg/reverse', '/tmp/tmp7lf5_feg/output.tsv.biom', '/tmp/tmp7lf5_feg/track.tsv', '/tmp/tmp7lf5_feg/filt_f', '/tmp/tmp7lf5_feg/filt_r', '0', '0', '29', '18', '2.0', '2.0', '2', '12', 'independent', 'consensus', '1.0', '12', '1000000']' returned non-zero exit status 1.
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/public1/home/sca0719/anaconda3/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 "/public1/home/sca0719/anaconda3/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 "/public1/home/sca0719/anaconda3/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 "/public1/home/sca0719/anaconda3/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.
Hey @Timmy21,
Your trim and truncation parameters looks kinda weird to me. Can you send a picture of your Demux.qzv quality score. Can you also send me your demux.qza file so that I can try to run this?
Also can you tell me a little bit more about your data? What did you do to get to this step and what type of sequencing data is this?
Hey @Timmy21 ,
Thank you for the information. Your qzv looks fine to me. Can you send me your demux.qza so that I can try to run dada2 and see if I can recreate your error.
Also can you tell me a little bit more about your data? What did you do to get to this step(i.e., did you import multiplexed or demuxplexed) and what type of sequencing data is this( i.e. is it illumina emp data or something else)?
Hi cherman2
My sequencing data is illumina paired-end data.
I got demux.qza after I run the comand :
qiime tools import --type 'SampleData[PairedEndSequencesWithQuality]' --input-path manifest.tsv --output-path demux.qza --input-format PairedEndFastqManifestPhred33V2
@Timmy21 ,
Looking back at the link I sent you, there might be something wrong with your manifest that is causing "Non-corresponding derep-class and dada-class objects".
Can you send your manifest.tsv? and again, I need your Demux.qza, I can not help you and unless you can provide me with all the files I request. If you dont want to share this publicly you can Direct Message me the Demux.qza
I think your issue might be a competing R installation.
Can you run echo $PATH in your command line and send me the output?
If your path has the installation of R before the QIIME 2 env that might be causing your computer to reference the R installation on your computer and not the one in the QIIME 2 env.
@Timmy21, I think it might be a path issue, it looks like your system-wide R installation is showing up before your conda environment in the path/you may have installed R in your system after creating your conda environment, the easiest thing might be to create a fresh QIIME2 environment using the instructions here(this will also get you the latest version). This is worth trying before doing anything more complicated.