# deblur denoise-other ITS error

**URL:** https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321
**Category:** Technical Support
**Tags:** deblur
**Created:** [August 5, 2023, 7:38am UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321 "2023-08-05T07:38:36Z")
**Posts on this page:** 9
**Page:** 1

<div class="post-metadata">

### Author: ![Rachel\_Burke](https://forum.qiime2.org/user_avatar/forum.qiime2.org/rachel_burke/32/16273_2.png) [@Rachel\_Burke](https://forum.qiime2.org/u/Rachel_Burke)
#### Post date: [August 5, 2023, 7:38am UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/1 "2023-08-05T07:38:37Z")

</div>

I am having issues with deblur denoise-other for assembling ITS. I used the UNITE database--specifically the "Includes singletons set as RefS (in dynamic files)" one, and only used the dynamic file (not 99 or 97) for the reference sequences. I've gotten different error messages by adjusting the trim length to 175, 250, 300, 350--either "no sequences passed the filter" or "max() is an empty argument." Does anyone have ideas of what might be wrong? Specific error message below.

- Version of QIIME 2 I am running: qiime2-2023.5
- What is the exact command or commands you ran:

qiime tools import   
--input-path manifest.tsv   
--type 'SampleData[PairedEndSequencesWithQuality]'   
--output-path paired-end-filtered-demux.qza   
--input-format PairedEndFastqManifestPhred33V2

qiime cutadapt trim-paired   
--i-demultiplexed-sequences paired-end-filtered-demux.qza   
--p-front-f GTGARTCATCGAATCTTTG   
--p-adapter-r CCTSCSCTTANTDATATGC   
--p-max-n 0   
--o-trimmed-sequences pef-demux.qza

qiime vsearch merge-pairs   
--i-demultiplexed-seqs pef-demux.qza   
--o-merged-sequences fil-demux-joined.qza

qiime tools import   
--type 'FeatureData[Sequence]'   
--input-path sh\_refs\_qiime\_ver9\_dynamic\_25.07.2023.fasta   
--output-path its\_dynamic\_ref\_seqs.qza

qiime deblur denoise-other   
--i-demultiplexed-seqs fil-demux-joined.qza   
--i-reference-seqs its\_dynamic\_ref\_seqs.qza   
--p-trim-length 500   
--p-sample-stats   
--p-jobs-to-start 20   
--o-representative-sequences rep-seqs.qza   
--o-table table.qza   
--o-stats deblur-stats.qza

- What is the exact error message?  
/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/bin/deblur:534: DeprecationWarning: The 'warn' method is deprecated, use 'warning' instead  
logger.warn('deblur version %s workflow started on %s' %  
/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/bin/deblur:536: DeprecationWarning: The 'warn' method is deprecated, use 'warning' instead  
logger.warn('parameters: %s' % locals())  
Traceback (most recent call last):  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/q2cli/commands.py", line 468, in **call**  
results = action(\*\*arguments)  
File "", line 2, in denoise\_other  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 274, in bound\_callable  
outputs = self._callable\_executor_(  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 509, in _callable\_executor_  
output\_views = self.\_callable(\*\*view\_args)  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/q2\_deblur/\_denoise.py", line 118, in denoise\_other  
return \_denoise\_helper(  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/q2\_deblur/\_denoise.py", line 187, in \_denoise\_helper  
raise ValueError("No sequences passed the filter. It is possible "  
ValueError: No sequences passed the filter. It is possible the trim\_length (500) 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.

Plugin error from deblur:

No sequences passed the filter. It is possible the trim\_length (500) 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.

---

<div class="post-metadata">

### Author: ![crusher083](https://forum.qiime2.org/user_avatar/forum.qiime2.org/crusher083/32/9597_2.png) [@crusher083](https://forum.qiime2.org/u/crusher083)
#### Post date: [August 5, 2023, 1:59pm UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/2 "2023-08-05T13:59:00Z")

</div>

Hello,

as error says, no sequences passed the filter. You'll need to check if sequences in `fill-demux-joined.qza` (open it with archiver software, go to folder `./data`) are longer than your trim length.

Cheers,  
V

---

<div class="post-metadata">

### Author: ![Rachel\_Burke](https://forum.qiime2.org/user_avatar/forum.qiime2.org/rachel_burke/32/16273_2.png) [@Rachel\_Burke](https://forum.qiime2.org/u/Rachel_Burke)
#### Post date: [August 7, 2023, 5:28am UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/3 "2023-08-07T05:28:50Z")

</div>

I tried lowering my trim length, but now I'm getting this error:

/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/bin/deblur:534: DeprecationWarning: The 'warn' method is deprecated, use 'warning' instead  
logger.warn('deblur version %s workflow started on %s' %  
/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/bin/deblur:536: DeprecationWarning: The 'warn' method is deprecated, use 'warning' instead  
logger.warn('parameters: %s' % locals())  
/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/deblur/workflow.py:92: UserWarning: input file /gpfs/scratch/opt/tmp/tmp5v37v\_3v/all.seqs.fa.no\_artifacts does not appear to be FASTA or FASTQ  
warnings.warn(msg, UserWarning)  
Traceback (most recent call last):  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/q2cli/commands.py", line 468, in **call**  
results = action(\*\*arguments)  
File "", line 2, in denoise\_other  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 274, in bound\_callable  
outputs = self._callable\_executor_(  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 509, in _callable\_executor_  
output\_views = self.\_callable(\*\*view\_args)  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/q2\_deblur/\_denoise.py", line 118, in denoise\_other  
return \_denoise\_helper(  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/q2\_deblur/\_denoise.py", line 194, in \_denoise\_helper  
table = \_load\_table(tmp)  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/q2\_deblur/\_denoise.py", line 61, in \_load\_table  
table.update\_ids(sid\_map, axis='sample', inplace=True)  
File "/gpfs/shared/apps\_local/python/3.7-1.20/envs/qiime2-2023.5/lib/python3.8/site-packages/biom/table.py", line 1405, in update\_ids  
str\_dtype = 'U%d' % max([len(v) for v in id\_map.values()])  
ValueError: max() arg is an empty sequence

Do you know how to fix this?

---

<div class="post-metadata">

### Author: ![crusher083](https://forum.qiime2.org/user_avatar/forum.qiime2.org/crusher083/32/9597_2.png) [@crusher083](https://forum.qiime2.org/u/crusher083)
#### Post date: [August 7, 2023, 9:32am UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/4 "2023-08-07T09:32:58Z")

</div>

Pay attention to this line:

```auto
 UserWarning: input file /gpfs/scratch/opt/tmp/tmp5v37v_3v/all.seqs.fa.no_artifacts does not appear to be FASTA or FASTQ

```

It seems like something is wrong either with your data or with the import step. Could you please attach your input files?

Cheers,  
V

---

<div class="post-metadata">

### Author: ![Rachel\_Burke](https://forum.qiime2.org/user_avatar/forum.qiime2.org/rachel_burke/32/16273_2.png) [@Rachel\_Burke](https://forum.qiime2.org/u/Rachel_Burke)
#### Post date: [August 7, 2023, 7:24pm UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/5 "2023-08-07T19:24:36Z")

</div>

> **[Google Drive: Sign-in](https://accounts.google.com/v3/signin/identifier?continue=https:%2F%2Fdrive.google.com%2Fdrive%2Ffolders%2F1f5P1UeouSWV3hAqwuyx_ehRYAqATxwuQ?usp%3Dsharing&followup=https:%2F%2Fdrive.google.com%2Fdrive%2Ffolders%2F1f5P1UeouSWV3hAqwuyx_ehRYAqATxwuQ?usp%3Dsharing&osid=1&passive=1209600&service=wise&flowName=GlifWebSignIn&flowEntry=ServiceLogin&dsh=S-986882776:1788800915843418)**
>
> Access Google Drive with a Google account (for personal use) or Google Workspace account (for business use).

Here is the link to most of the data. I couldn't upload the paired-end-filtered-demux.qza, as it is too large. Also, the original sequencing data files are too large to include, so I just included one. They have been trimmed and filtered a little using CLC Genomics Workbench (base calling error probability of 0.05, maximum of 2 ambiguities per read, trimmed any adaptor sequences present).

 ![image](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/7/a/7a9aac1a24dfd5d29baead24a6969d70d170db39.png)

Here are the sizes of each of the files, if that helps.

Thank you so much for your help, I appreciate it greatly! I'm a junior in college, so I'm very new to bioinformatics 🙂

The sequencing data has both ITS and 16S, each with different primers... I'm running deblur denoise-16S using the same pipeline (except for different primer removal and not needing a reference sequences file) and the denoise-16S has worked on my data (or at least didn't give me any errors like the denoise-other for ITS did).

I've also tried dada2 denoise-paired but lost all my data.

---

<div class="post-metadata">

### Author: ![gregcaporaso](https://forum.qiime2.org/user_avatar/forum.qiime2.org/gregcaporaso/32/17769_2.png) [@gregcaporaso](https://forum.qiime2.org/u/gregcaporaso)
#### Post date: [August 8, 2023, 9:49pm UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/7 "2023-08-08T21:49:07Z")

</div>

> [@Rachel\_Burke](#):
>
> paired-end-filtered-demux.qza

@Rachel_Burke, have you run `qiime demux summarize` on that file? That will give you a summary of the demultiplexed sequences that may provide some clues here. Could you generate that, poke around to see if you can identify what's going wrong, and then share that file with us if you need more assistance?

---

<div class="post-metadata">

### Author: ![Rachel\_Burke](https://forum.qiime2.org/user_avatar/forum.qiime2.org/rachel_burke/32/16273_2.png) [@Rachel\_Burke](https://forum.qiime2.org/u/Rachel_Burke)
#### Post date: [August 9, 2023, 2:31am UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/8 "2023-08-09T02:31:13Z")

</div>

[visualization.qzv](https://forum.qiime2.org/uploads/short-url/bZJ8i2P0b4SU0GaK3CTOEZSfGsW.qzv) (315.5 KB)  
Here's what I've got for the demux summarize after merging with vsearch/before denoising with deblur. From what I found googling, it doesn't look too bad.

This means that the problem must be happening in the deblur denoise step, right? I tried running deblur denoise-16S on those reads (with different cutadapt primer removal), which worked but once I got to the OTU table, I only had 1 read left. I tried dada2 again but increased the expected error of the reverse reads to 5 (the reverse reads are lower quality), but also ended with 1 read in my OTU table.

---

<div class="post-metadata">

### Author: ![gregcaporaso](https://forum.qiime2.org/user_avatar/forum.qiime2.org/gregcaporaso/32/17769_2.png) [@gregcaporaso](https://forum.qiime2.org/u/gregcaporaso)
#### Post date: [August 9, 2023, 4:15pm UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/9 "2023-08-09T16:15:02Z")

</div>

@Rachel_Burke, that looks like the summary prior to merging with vsearch. Could you share this same visualization generated on the output of the merging step?

Another thought: have you tried comparing sequences in your reference data with your sequence reads, for example using command line BLAST? I am wondering if the sequences are in the same orientation (I'm not sure if deblur adjusts for orientation).

---

<div class="post-metadata">

### Author: ![system](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/2/1/21af5fe23cb6f4579467c66a9ed94e55274ca7bd.svg) [@system](https://forum.qiime2.org/u/system)
#### Post date: [September 9, 2023, 10:15pm UTC](https://forum.qiime2.org/t/deblur-denoise-other-its-error/27321/10 "2023-09-09T22:15:33Z")

</div>

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.
