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/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
R version 4.1.0 (2021-05-18)
Loading required package: Rcpp
DADA2: 1.22.0 / Rcpp: 184.108.40.206 / 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
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
File "/public1/home/sca0719/anaconda3/envs/qiime2-2022.2/lib/python3.8/site-packages/q2_dada2/_denoise.py", line 36, in run_commands
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.
Can you give the command you ran to get this error and your demux file?
I have seen errors like this in the github for dada2: Non-corresponding derep-class and dada-class objects. · Issue #866 · benjjneb/dada2 · GitHub . That might be helpful to you.
My command as follow:
qiime dada2 denoise-paired
--p-trim-left-f 29 --p-trim-left-r 18
--p-trunc-len-f 0 --p-trunc-len-r 0
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?
per-sample-fastq-counts.tsv (344 Bytes)
demux.qzv (311.7 KB)
This is my Demux.qzv rusult. thanks for your continued attention
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)?
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
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
Hi again @Timmy21 ,
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.
Let me know!
run the command "echo PATH"
(qiime2-2022.2) [[email protected]%bscc-a5 qiime2] echo $PATH
An off-topic reply has been split into a new topic: Incompatible Classes when running q2-DADA2
Please keep replies on-topic in the future.
@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.
This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.