# Dada2 Merging Issue

**URL:** https://forum.qiime2.org/t/dada2-merging-issue/27216
**Category:** User Support
**Tags:** dada2, merge
**Created:** [July 26, 2023, 9:13pm UTC](https://forum.qiime2.org/t/dada2-merging-issue/27216 "2023-07-26T21:13:50Z")
**Posts on this page:** 5
**Page:** 1

<div class="post-metadata">

### Author: ![elmo2420](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/e/c2a13f/32.png) [@elmo2420](https://forum.qiime2.org/u/elmo2420)
#### Post date: [July 26, 2023, 9:13pm UTC](https://forum.qiime2.org/t/dada2-merging-issue/27216/1 "2023-07-26T21:13:50Z")

</div>

My sequences aren't merging and I'm not sure where I am going wrong. For my parameters, I am trying to stay within 140-180 (forward) and 200-215(reverse) positions to be able to achieve an 80% for the input passed filter percentage. I also took into consideration that my quality score shouldn't be below 30 in the 25th percentile. I provided the command lines of my parameters:

```auto
qiime dada2 denoise-paired \
  --i-demultiplexed-seqs demultiplexed-sequences.qza \
  --p-trunc-len-f 237 \
  --p-trunc-len-r 175 \
  --p-trim-left-f 5 \
  --p-trim-left-r 5 \
  --p-max-ee-f 2.0 \
  --p-max-ee-r 2.0 \
  --o-representative-sequences rep-seqs-237-175.qza \
  --o-table feature-table-237-175.qza \
  --o-denoising-stats dada2-stats-237-175.qza

```

Here is my demultiplexed sequences interactive quality plot:

 ![Screenshot 2023-07-26 at 1.37.45 PM](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/7/5/75e3b3648d7e16e930835087115e92057460dc3d.jpeg)

And my dada2 summary results:

 ![Screenshot 2023-07-26 at 1.37.53 PM](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/f/6/f66093b01dfb337aba0341aad9697dbd51adfa4d.png)

One of my questions is if the --p-max-ee parameter could be causing this issue but I'm really sure about it. Any suggestion to fix this issue would be great! Thank you in advance! 🙂

---

<div class="post-metadata">

### Author: ![colinbrislawn](https://forum.qiime2.org/user_avatar/forum.qiime2.org/colinbrislawn/32/6221_2.png) [@colinbrislawn](https://forum.qiime2.org/u/colinbrislawn)
#### Post date: [July 27, 2023, 12:53am UTC](https://forum.qiime2.org/t/dada2-merging-issue/27216/2 "2023-07-27T00:53:46Z")

</div>

Hello Lei,

Welcome to the forums! :qiime2:

I like your approach. Let's get these reads to merge!

> [@elmo2420](#):
>
> One of my questions is if the --p-max-ee parameter could be causing this issue but I'm really sure about it.

Nope, that setting only changes the filtered reads before merging.

> [@elmo2420](#):
>
> to be able to achieve an 80% for the input passed filter percentage. I also took into consideration that my quality score shouldn't be below 30 in the 25th percentile.

Sounds good. But because the reads are not merging, get creative.

Because almost no reads merge, this makes me think the reads after trimming are not covering the amplicon you sequenced.

How long is your amplicon? Is it longer than 237+175==412 basepairs long?

---

<div class="post-metadata">

### Author: ![elmo2420](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/e/c2a13f/32.png) [@elmo2420](https://forum.qiime2.org/u/elmo2420)
#### Post date: [July 27, 2023, 3:58am UTC](https://forum.qiime2.org/t/dada2-merging-issue/27216/3 "2023-07-27T03:58:58Z")

</div>

Hi Colin,

Thank you for your response and the information you provided! The amplicon is 300 bp. I recently tried it again and I changed the `--p-trunc-len` to:

```auto
  --p-trunc-len-f 256 \
  --p-trunc-len-r 202 \

```

I used the equation `--p-trunc-len-f + --p-trunc-len-r - amplicon length` and my overlap was 150 bp, which got some of my reads to merge shown in the screenshot

 ![Screenshot 2023-07-26 at 7.27.47 PM](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/1/6/167ed9f5b8c387fb7ed4afc1fb655114704b1c07.jpeg)

---

<div class="post-metadata">

### Author: ![colinbrislawn](https://forum.qiime2.org/user_avatar/forum.qiime2.org/colinbrislawn/32/6221_2.png) [@colinbrislawn](https://forum.qiime2.org/u/colinbrislawn)
#### Post date: [July 27, 2023, 11:23am UTC](https://forum.qiime2.org/t/dada2-merging-issue/27216/4 "2023-07-27T11:23:06Z")

</div>

I'm glad to hear that!

Merging a 300 bp amplicon is hard. While keeping ~20% of the reads is unfortunate, looks like most of your samples have \>10k reads, which is pretty good.

> [@elmo2420](#):
>
> I used the equation `--p-trunc-len-f + --p-trunc-len-r - amplicon length` and my overlap was 150 bp

DADA2 tends to work better with smaller overlaps, so that may be worth trying to get more of your reads to merge.

---

<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: [August 27, 2023, 5:23pm UTC](https://forum.qiime2.org/t/dada2-merging-issue/27216/5 "2023-08-27T17:23:18Z")

</div>

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