# DADA2 denoise : No reads passed the filter (PacBio sequencing data)

**URL:** https://forum.qiime2.org/t/dada2-denoise-no-reads-passed-the-filter-pacbio-sequencing-data/21609
**Category:** Technical Support
**Tags:** dada2
**Created:** [December 9, 2021, 7:16am UTC](https://forum.qiime2.org/t/dada2-denoise-no-reads-passed-the-filter-pacbio-sequencing-data/21609 "2021-12-09T07:16:10Z")
**Posts on this page:** 4
**Page:** 1

<div class="post-metadata">

### Author: ![Kim\_Jihoon](https://forum.qiime2.org/user_avatar/forum.qiime2.org/kim_jihoon/32/12296_2.png) [@Kim\_Jihoon](https://forum.qiime2.org/u/Kim_Jihoon)
#### Post date: [December 9, 2021, 7:16am UTC](https://forum.qiime2.org/t/dada2-denoise-no-reads-passed-the-filter-pacbio-sequencing-data/21609/1 "2021-12-09T07:16:10Z")

</div>

Hi! I'm new to QIIME2 and I have some problem with denoising my PacBio sequencing data.

I am running qiime2 2021.11 with conda in Ubuntu. I'm dealing with my 18S rRNA PacBio sequencing data by PacBio Sequel I sequencing platform. I've searched many topics related to this issue but I can't solve this problem so far.  
When I run dada2 denoise-single with .qza artifact file from my sequenced data, there are No Reads Passed The Filter. Command I used is:  
**qiime dada2 denoise-single --i-demultiplexed-seqs trimmed\_reads.qza --output-dir dada2\_output --p-n-threads 0 --p-trunc-len 0 --p-trunc-q 3 --p-max-ee 3 --verbose**

I tried changing parameters of trunc-q from 3 to 30, --p-max-ee from 3 to 30 or turn on/off the chimera-method, but still no reads passed the filter like following:  
[debug\_info.txt](https://forum.qiime2.org/uploads/short-url/a5mboA1RDxeAGh0VGDyiBD9hU9F.txt) (12.3 KB)  
As PacBio sequel I platform generates data with quality score "!", which means lowest quality score, to all sequenced data (because PacBio said as quality is already high enough, so no need to generate quality score), I've changed all "!" quality score to "#" (which means quality score = 2) to avoid all data being removed. So my demux.qza file possess interactive quality plot like this:

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

Thank you for reading this issue. I really appreciate your comments on this issue!

---

<div class="post-metadata">

### Author: ![llenzi](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/l/9fc29f/32.png) [@llenzi](https://forum.qiime2.org/u/llenzi)
#### Post date: [December 9, 2021, 9:22am UTC](https://forum.qiime2.org/t/dada2-denoise-no-reads-passed-the-filter-pacbio-sequencing-data/21609/2 "2021-12-09T09:22:02Z")

</div>

Hi @Kim_Jihoon,  
welcome in the forum!

For what I know, DADA2 can be used to process PacBio sequences, but you need to work in R.  
The dada2 denoise single plug in is not tailored for that.  
Few tutorial you may want to look at, are linked in this ResearchGate topic:  
[https://www.researchgate.net/post/Pacbio\_sequences\_and\_DADA2\_pipeline\_compatible](https://www.researchgate.net/post/Pacbio_sequences_and_DADA2_pipeline_compatible)

Hope you manage to work with this.  
Best wishes  
Luca

---

<div class="post-metadata">

### Author: ![Kim\_Jihoon](https://forum.qiime2.org/user_avatar/forum.qiime2.org/kim_jihoon/32/12296_2.png) [@Kim\_Jihoon](https://forum.qiime2.org/u/Kim_Jihoon)
#### Post date: [December 9, 2021, 9:24am UTC](https://forum.qiime2.org/t/dada2-denoise-no-reads-passed-the-filter-pacbio-sequencing-data/21609/3 "2021-12-09T09:24:32Z")

</div>

Hi @llenzi, thank you for your reply!

I'll check what you linked and try to do with it.

Thank you!

---

<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: [January 9, 2022, 3:24pm UTC](https://forum.qiime2.org/t/dada2-denoise-no-reads-passed-the-filter-pacbio-sequencing-data/21609/4 "2022-01-09T15:24:50Z")

</div>

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