# Minimum read for denoising process

**URL:** https://forum.qiime2.org/t/minimum-read-for-denoising-process/8210
**Category:** General Discussion
**Tags:** dada2, deblur
**Created:** [February 6, 2019, 4:26pm UTC](https://forum.qiime2.org/t/minimum-read-for-denoising-process/8210 "2019-02-06T16:26:58Z")
**Posts on this page:** 5
**Page:** 1

<div class="post-metadata">

### Author: ![devonorourke](https://forum.qiime2.org/user_avatar/forum.qiime2.org/devonorourke/32/3110_2.png) [@devonorourke](https://forum.qiime2.org/u/devonorourke)
#### Post date: [February 6, 2019, 4:26pm UTC](https://forum.qiime2.org/t/minimum-read-for-denoising-process/8210/1 "2019-02-06T16:26:58Z")

</div>

Hopefully a quick question:  
From what I can gather on posts @wasade has commented (ex. [here](https://forum.qiime2.org/t/figuring-out-why-sequences-are-lost-in-deblur/6718/9)) Deblur has parameters to define how many reads are required to be considered going into filtering (`--p-min-size`) and how many reads are required following filtering to be included in the resulting frequency table (`--p-min-reads`).  
I can't find any such information on DADA2; I'm specifically interested in the [paired-end approach](https://docs.qiime2.org/2019.1/plugins/available/dada2/denoise-paired/) if it makes a difference. I'm wondering (1) whether DADA2 considers reads of any abundance when building its error models, and (2) following the paired-end joining, if it removes any reads that pass its filtering but have less than N total merged read pairs. Perhaps @benjjneb can help.  
Thanks very much!

---

<div class="post-metadata">

### Author: ![benjjneb](https://forum.qiime2.org/user_avatar/forum.qiime2.org/benjjneb/32/1602_2.png) [@benjjneb](https://forum.qiime2.org/u/benjjneb)
#### Post date: [February 6, 2019, 9:50pm UTC](https://forum.qiime2.org/t/minimum-read-for-denoising-process/8210/2 "2019-02-06T21:50:37Z")

</div>

> [@devonorourke](#):
>
> (1) whether DADA2 considers reads of any abundance when building its error models

Yes.

> [@devonorourke](#):
>
> (2) following the paired-end joining, if it removes any reads that pass its filtering but have less than N total merged read pairs.

No. But note there is an effective threshold that removes all singletons due to the division criteria in the core algorithm. But there is no additional threshold enforced later, that is left up to the user if they so desire.

---

<div class="post-metadata">

### Author: ![devonorourke](https://forum.qiime2.org/user_avatar/forum.qiime2.org/devonorourke/32/3110_2.png) [@devonorourke](https://forum.qiime2.org/u/devonorourke)
#### Post date: [February 6, 2019, 10:18pm UTC](https://forum.qiime2.org/t/minimum-read-for-denoising-process/8210/3 "2019-02-06T22:18:51Z")

</div>

Thanks @benjjneb,

(1) I should have been more specific. Deblur has a parameter where it _ignores_ reads below a user defined (or default) value; I was wondering if DADA2 has a specific minimum value, or if _all_ read abundances are considered.

( 2) This threshold... it removes singletons? But it wouldn't remove doubletons then?

Cheers!

---

<div class="post-metadata">

### Author: ![benjjneb](https://forum.qiime2.org/user_avatar/forum.qiime2.org/benjjneb/32/1602_2.png) [@benjjneb](https://forum.qiime2.org/u/benjjneb)
#### Post date: [February 7, 2019, 2:53pm UTC](https://forum.qiime2.org/t/minimum-read-for-denoising-process/8210/4 "2019-02-07T14:53:52Z")

</div>

> [@devonorourke](#):
>
> (1) I should have been more specific. Deblur has a parameter where it _ignores_ reads below a user defined (or default) value; I was wondering if DADA2 has a specific minimum value, or if _all_ read abundances are considered.

The q2-dada2 plugin does not have such a parameter.

> [@devonorourke](#):
>
> ( 2) This threshold… it removes singletons? But it wouldn’t remove doubletons then?

Yes, although doubletons will still probably be diminished somewhat in paired end data, as they must be detected in both the F and R reads independently.

---

<div class="post-metadata">

### Author: ![devonorourke](https://forum.qiime2.org/user_avatar/forum.qiime2.org/devonorourke/32/3110_2.png) [@devonorourke](https://forum.qiime2.org/u/devonorourke)
#### Post date: [February 7, 2019, 3:16pm UTC](https://forum.qiime2.org/t/minimum-read-for-denoising-process/8210/5 "2019-02-07T15:16:21Z")

</div>

Thanks @benjjneb - very helpful. Appreciate the quick response!
