Hi, I am trying to work on the sequence quality control using few of my sequences using DADA2. Also, I have no prior experience in this.
The command I used is as follows:
qiime dada2 denoise-paired --i-demultiplexed-seqs demux.qza --p-trim-left-f 20 --p-trim-left-r 20 --p-trunc-len-f 100 --p-trunc-len-r 100 --p-n-threads 12 --o-representative-sequences rep-seqs.qza --o-table table.qza
But I am getting an error that says:
Plugin error from dada2:
An error was encountered while running DADA2 in R (return code 1), please inspect stdout and stderr to learn more.
Debug info has been saved to /tmp/qiime2-q2cli-err-uzeztcgz.log
On checking the log, it shows the following:
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/tmp0fmgeqwf/forward /tmp/tmp0fmgeqwf/reverse /tmp/tmp0fmgeqwf/output.tsv.biom /tmp/tmp0fmgeqwf/filt_f /tmp/tmp0fmgeqwf/filt_r 100 100 20 20 2.0 2 consensus 1.0 12 1000000
R version 3.4.1 (2017-06-30)
Loading required package: Rcpp
DADA2 R package version: 1.6.0
- Filtering Error in filterAndTrim(unfiltsF, filtsF, unfiltsR, filtsR, truncLen = c(truncLenF, :
These are the errors (up to 5) encountered in individual cores...
Error : cannot allocate vector of size 95.4 Mb
Error : cannot allocate vector of size 95.4 Mb
Execution halted
Traceback (most recent call last):
File "/home/qiime2/miniconda/envs/qiime2-2018.2/lib/python3.5/site-packages/q2_dada2/_denoise.py", line 179, in denoise_paired
run_commands([cmd])
File "/home/qiime2/miniconda/envs/qiime2-2018.2/lib/python3.5/site-packages/q2_dada2/_denoise.py", line 35, in run_commands
subprocess.run(cmd, check=True)
File "/home/qiime2/miniconda/envs/qiime2-2018.2/lib/python3.5/subprocess.py", line 398, in run
output=stdout, stderr=stderr)
subprocess.CalledProcessError: Command '['run_dada_paired.R', '/tmp/tmp0fmgeqwf/forward', '/tmp/tmp0fmgeqwf/reverse', '/tmp/tmp0fmgeqwf/output.tsv.biom', '/tmp/tmp0fmgeqwf/filt_f', '/tmp/tmp0fmgeqwf/filt_r', '100', '100', '20', '20', '2.0', '2', '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 "/home/qiime2/miniconda/envs/qiime2-2018.2/lib/python3.5/site-packages/q2cli/commands.py", line 246, in call
results = action(**arguments)
File "", line 2, in denoise_paired
File "/home/qiime2/miniconda/envs/qiime2-2018.2/lib/python3.5/site-packages/qiime2/sdk/action.py", line 228, in bound_callable
output_types, provenance)
File "/home/qiime2/miniconda/envs/qiime2-2018.2/lib/python3.5/site-packages/qiime2/sdk/action.py", line 363, in callable_executor
output_views = self._callable(**view_args)
File "/home/qiime2/miniconda/envs/qiime2-2018.2/lib/python3.5/site-packages/q2_dada2/_denoise.py", line 194, in denoise_paired
" and stderr to learn more." % e.returncode)
Exception: An error was encountered while running DADA2 in R (return code 1), please inspect stdout and stderr to learn more.
Before this I had imported the metadata using the command:
qiime tools import \
--type 'SampleData[PairedEndSequencesWithQuality]'
--input-path pe-64-manifest
--output-path paired-end-demux.qza
--source-format PairedEndFastqManifestPhred33