Hi good folks!
I've been working with QIIME 2 before, and everything has worked flawlessly (this forum is a goldmine, thanks!)
On my recent project, I'm stuck at deblur for some reason I can't seem to figure out. I work in QIIME2/2020.11
Heres my code, following this SOP:
qiime tools import \
--type SampleData[PairedEndSequencesWithQuality] \
--input-path /raw_data \
--output-path /treated_data/reads.qza \
--input-format CasavaOneEightSingleLanePerSampleDirFmt
qiime cutadapt trim-paired \
--i-demultiplexed-sequences /treated_data/reads.qza \
--p-cores $NCORES \
--p-front-f CCAGCASCYGCGGTAATTCC \
--p-front-r ACTTTCGTTCTTGATYRATGA \
--p-discard-untrimmed \
--p-no-indels \
--o-trimmed-sequences /treated_data/reads_trimmed.qza
qiime vsearch join-pairs \
--i-demultiplexed-seqs /treated_data/reads_trimmed.qza \
--o-joined-sequences /treated_data/reads_trimmed_joined.qza
qiime quality-filter q-score \
--i-demux /treated_data/reads_trimmed_joined.qza \
--o-filter-stats /treated_data/filt_stats.qza \
--o-filtered-sequences /treated_data/reads_trimmed_joined_filt.qza
qiime deblur denoise-other \
--i-demultiplexed-seqs /treated_data/reads_trimmed_joined.qza \
--i-reference-seqs /gb203_pr2_all_10_28_99p_clean_prob-rm.qza \
--p-trim-length 1 \
--p-sample-stats \
--p-jobs-to-start $NCORES \
--p-min-reads 1 \
--output-dir /treated_data/deblur_output
My error reads:
Traceback (most recent call last):
File "/cluster/software/QIIME2/2020.11/lib/python3.6/site-packages/q2cli/commands.py", line 329, in __call__
results = action(**arguments)
File "<decorator-gen-510>", line 2, in denoise_other
File "/cluster/software/QIIME2/2020.11/lib/python3.6/site-packages/qiime2/sdk/action.py", line 245, in bound_callable
output_types, provenance)
File "/cluster/software/QIIME2/2020.11/lib/python3.6/site-packages/qiime2/sdk/action.py", line 390, in _callable_executor_
output_views = self._callable(**view_args)
File "/cluster/software/QIIME2/2020.11/lib/python3.6/site-packages/q2_deblur/_denoise.py", line 130, in denoise_other
hashed_feature_ids=hashed_feature_ids)
File "/cluster/software/QIIME2/2020.11/lib/python3.6/site-packages/q2_deblur/_denoise.py", line 192, in _denoise_helper
"of the data being denoised." % trim_length)
ValueError: No sequences passed the filter. It is possible the trim_length (1) may exceed the longest sequence, that all of the sequences are artifacts like PhiX or adapter, or that the positive reference used is not representative of the data being denoised.
I have triple checked my primers, and they are correct. I have never worked with denoise-other
(all my former projects was 16S rRNA).
All the best,
Stian