Hi there,
I am trying to use the decontam-identify method to find potential contaminants in ITS data, produced from Illumina sequencing.
I processed all of my data in Qiime2-2023.5, using an Ubuntu distribution via WSL. I had paired-end, demultiplexed reads which I imported, used the cutadapt plugin to remove primers and then de-noised using the DADA2 pipeline. I then assigned taxonomy using the feature-classifier plugin (naive bayes classifier trained on UNITE 99 developers ref seqs).
I'm planning on using all 3 methods (prevalence, frequency and combined), comparing what ASVs each identify as potential contamination, and then deciding which scores to use to filter out contamination.
I have tried to use decontam-identify with my feature table output from dada2 denoising, using the following code:
qiime quality-control decontam-identify --i-table table2.qza --m-metadata-file Metadata.tsv --p-method 'combined' --p-freq-concentration-column DNA-Conc --p-prev-control-column Sample-Type --p-prev-control-indicator "Control" --o-decontam-scores decontam-scores-comb-method.qza --verbose
but am receiving an error message:
Plugin error from quality-control:
An error was encountered while running Decontam in R (return code 1), please inspect stdout and stderr to learn more.
I will attach the debug info in a .txt file!
decontam-identify debug info.txt (3.6 KB)
I've had a bit of a read around but can't figure out what may be causing this - am a bit of a novice with QIIME and only slightly more advanced with R, but I thought it may be to do with the version of R on my PC?
I wasn't sure if I have to download the decontam package on R before running this command? I tried to do this
install.packages ("decontam")
but it returned
Warning in install.packages :
package ‘decontam’ is not available for this version of R
(which is R v4.3.1)
Any help would be great!
Thanks!