# Recommendation for processing forward and reverse reads with low quality issues

**URL:** https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730
**Category:** User Support
**Tags:** dada2, quality
**Created:** [April 11, 2018, 11:49am UTC](https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730 "2018-04-11T11:49:20Z")
**Posts on this page:** 6
**Page:** 1

<div class="post-metadata">

### Author: ![jorge\_rico](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/j/71c47a/32.png) [@jorge\_rico](https://forum.qiime2.org/u/jorge_rico)
#### Post date: [April 11, 2018, 11:49am UTC](https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730/1 "2018-04-11T11:49:20Z")

</div>

Hi, I have a 16S rRNA amplicon data set (EMP standard primers) that was sequenced at a facility and the returned reads aren't great quality, especially the reverse reads. See screenshot and uploaded .qzv. My questions are the following:

1. would you recommend merging these reads (as opposed to just using forward reads)
2. with the Dada2 options ([denoise-single: Denoise and dereplicate single-end sequences — QIIME 2 2018.2.0 documentation](https://docs.qiime2.org/2018.2/plugins/available/dada2/denoise-single/)), which are the best to use to ensure my final sequence reads are high quality? It seems like the most relevant options are --p-trunc-len, --p-trim-left, --p-max-ee, --p-trunc-q.

It seems like --p-max-ee would kick out most really bad sequences, is that correct?

With --p-trunc-q having a default of 2, will the default actually do anything given that I don't have reads with quality that dips below 2, which seems really low. Am I understanding this correctly?

Alternatively, would you recommend processing these sequences outside of qiime2 because the quality is poor.

Finally, do you think I should inquire to the sequencing facility about this quality? Does it seem unacceptable?

-Jorge

 ![10%20PM](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/2X/2/2171ad80bc564f8e4427d2781b06b524ad663256.jpg)

[demux-paired-end.qzv](https://cdck-file-uploads-global.s3.dualstack.us-west-2.amazonaws.com/flex002/uploads/qiime21/original/2X/1/13bbb76d1c215bd59e48652fd66a96edf83d56da.qzv) (267.5 KB)

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [April 11, 2018, 11:52am UTC](https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730/2 "2018-04-11T11:52:57Z")

</div>



---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [April 11, 2018, 9:28pm UTC](https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730/3 "2018-04-11T21:28:00Z")

</div>

Hi @jorge_rico,

> [@jorge\_rico](#):
>
> would you recommend merging these reads (as opposed to just using forward reads)

I would either just use the forward reads, or attempt to merge the first 220 nt of forward with 100 nt of reverse. Definitely trim left on each of these to remove the low-quality bases at the start of each read.

> [@jorge\_rico](#):
>
> which are the best to use to ensure my final sequence reads are high quality? It seems like the most relevant options are --p-trunc-len, --p-trim-left, --p-max-ee, --p-trunc-q.

dada2 is going to take care of denoising these reads. The greater issue is ensuring that these reads are trimmed so that low-quality bases at the 3' end of each read aren't causing those reads to be tossed. So yes, --p-trunc-len and --p-trim-left will be the main parameters to use. You could also try playing around with --p-max-ee to squeeze out some more reads (since it looks like errors are interspersed through your reads).

> [@jorge\_rico](#):
>
> It seems like --p-max-ee would kick out most really bad sequences, is that correct?

correct. it pre-filters the reads prior to denoising.

> [@jorge\_rico](#):
>
> With --p-trunc-q having a default of 2, will the default actually do anything given that I don’t have reads with quality that dips below 2, which seems really low. Am I understanding this correctly?

correct. By default, effectively no trimming is performed. That's the goal. (and personally I like trimming at specific sites rather than using --p-trunc-q which will result in reads/ASVs of different lengths).

> [@jorge\_rico](#):
>
> Alternatively, would you recommend processing these sequences outside of qiime2 because the quality is poor.

No. This is not a qiime2 issue, it is a problem with your reads, so going outside of qiime2 will not help anything necessarily. If you want to do qiime1-style quality filtering and otu picking (much more permissive than dada2/deblur), this can be done with q2-quality-filter and q2-vsearch.

> [@jorge\_rico](#):
>
> Finally, do you think I should inquire to the sequencing facility about this quality? Does it seem unacceptable?

The quality is disappointing but not the worst I've seen. You could chat with them to see if they would be willing to re-run, but I am not sure that this is necessarily "unacceptable". I'd give it a shot and see how much data you can squeeze out...

I hope that helps!

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [April 11, 2018, 9:28pm UTC](https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730/4 "2018-04-11T21:28:03Z")

</div>



---

<div class="post-metadata">

### Author: ![jorge\_rico](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/j/71c47a/32.png) [@jorge\_rico](https://forum.qiime2.org/u/jorge_rico)
#### Post date: [May 7, 2018, 1:03pm UTC](https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730/5 "2018-05-07T13:03:20Z")

</div>

Thank you very much for all your responses, I ended up using the forward reads as you suggested!

---

<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: [June 7, 2018, 7:03pm UTC](https://forum.qiime2.org/t/recommendation-for-processing-forward-and-reverse-reads-with-low-quality-issues/3730/6 "2018-06-07T19:03:26Z")

</div>

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