# Cutadapt primer removal questions

**URL:** https://forum.qiime2.org/t/cutadapt-primer-removal-questions/32762
**Category:** Other Bioinformatics Tools
**Created:** [March 15, 2025, 5:21pm UTC](https://forum.qiime2.org/t/cutadapt-primer-removal-questions/32762 "2025-03-15T17:21:16Z")
**Posts on this page:** 5
**Page:** 1

<div class="post-metadata">

### Author: ![Phoe](https://forum.qiime2.org/user_avatar/forum.qiime2.org/phoe/32/18852_2.png) [@Phoe](https://forum.qiime2.org/u/Phoe)
#### Post date: [March 15, 2025, 5:21pm UTC](https://forum.qiime2.org/t/cutadapt-primer-removal-questions/32762/1 "2025-03-15T17:21:16Z")

</div>

Hello everyone!  
I was trying to cut my adapters and primers sequences in my FastQs (V3-V4 region).

I usually use Fastp for adapter removal but this time i wanted to use cutadapt just to learn one more tool.

So i used 2 different commands that i think can both work but i'm not sure if they are 100% correct.

Considering that the adapter+primer sequence are:

16S Amplicon PCR Forward Primer 5'

TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG **CCTACGGGNGGCWGCAG**

16S Amplicon PCR Reverse Primer= 5'

GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG **GACTACHVGGGTATCTAATCC**

(Stronged sequence is the primer)

I used these 2 commands:  
cutadapt -g CCTACGGGNGGCWGCAG -G GACTACHVGGGTATCTAATCC -a CTGCWGCCNCCCGTAGG -A GGATTAGATACCCBDGTAGTC -Q 0 --discard-untrimmed --pair-filter=any -m 50 -j 30 -o cleaned\_4\_20/test1\_r1 -p cleaned\_4\_20/test1\_r2 r1.fastq.gz r2.fastq.gz

cutadapt -a ^CCTACGGGNGGCWGCAG...GGATTAGATACCCBDGTAGTC -A ^GACTACHVGGGTATCTAATCC...CTGCWGCCNCCCGTAGG -Q 0 --discard-untrimmed --pair-filter=any -m 50 -j 30 -o cleaned\_4\_20/r1\_trimmed\_prova2.fastq.gz -p cleaned\_4\_20/r2\_trimmed\_prova2.fastq.gz r1.fastq.gz r2.fastq.gz

My consideration are:

- Assuming that the read-through always happens, i think that the 2nd command is more correct, but this does not always happens; and for this reason, 1st command should also good. Am i missing something?  
What would you use and why?

Also, is primer removal mandatory? I've read on DADA2 that it can create problems if ambiguous characters are there, but in my case there shouldn't be any ambiguous character in my reads.  
Thank you for your time and suggestions.

---

<div class="post-metadata">

### Author: ![SoilRotifer](https://forum.qiime2.org/user_avatar/forum.qiime2.org/soilrotifer/32/21071_2.png) [@SoilRotifer](https://forum.qiime2.org/u/SoilRotifer)
#### Post date: [March 15, 2025, 5:33pm UTC](https://forum.qiime2.org/t/cutadapt-primer-removal-questions/32762/2 "2025-03-15T17:33:14Z")

</div>

Hi @Phoe,

> [@Phoe](#):
>
> Also, is primer removal mandatory?

Yes. It is not only about the presence of the ambiguous IUPAC bases, but also the fact that those sequences are still primer sequences that have been incorporated into the PCR / sequencing product. Thus, they are not actually the sequences from the target organism. Remember PCR primers are "leaky" and can still bind to and amplify despite a few mismatches, which are basically masked by the incorporated primer.

Also you normally do not need to add the adapter sequence into cutadapt unless you are getting substantial read-through. You should be able to get away with simply :

```auto
cutadapt -g CCTACGGGNGGCWGCAG -G GACTACHVGGGTATCTAATCC ...

```

Also, you can simply use [cutadapt](https://amplicon-docs.readthedocs.io/en/latest/references/plugin-reference/plugins/cutadapt/trim-paired.html#q2-action-cutadapt-trim-paired) via QIIME 2.

---

<div class="post-metadata">

### Author: ![Phoe](https://forum.qiime2.org/user_avatar/forum.qiime2.org/phoe/32/18852_2.png) [@Phoe](https://forum.qiime2.org/u/Phoe)
#### Post date: [March 15, 2025, 6:33pm UTC](https://forum.qiime2.org/t/cutadapt-primer-removal-questions/32762/3 "2025-03-15T18:33:57Z")

</div>

Thank you for your answer @SoilRotifer !  
So the final command would be:

cutadapt -g CCTACGGGNGGCWGCAG -G GACTACHVGGGTATCTAATCC -Q 0 --discard-untrimmed --pair-filter=any -m 50 -j 30 -o cleaned\_4\_20/test1\_r1 -p cleaned\_4\_20/test1\_r2 r1.fastq.gz r2.fastq.gz

Without the 3' primer sequences? If that's a yes, the answer is it in the read-through like we were discussing earlier?

Also I read that the best practice would be:

- Adapter removal
- Merging
- Quality trimming
- Import in qiime2 and analysis  
Does it change anything if I do everything inside qiime2?

Thank you again!

---

<div class="post-metadata">

### Author: ![SoilRotifer](https://forum.qiime2.org/user_avatar/forum.qiime2.org/soilrotifer/32/21071_2.png) [@SoilRotifer](https://forum.qiime2.org/u/SoilRotifer)
#### Post date: [March 15, 2025, 6:52pm UTC](https://forum.qiime2.org/t/cutadapt-primer-removal-questions/32762/4 "2025-03-15T18:52:01Z")

</div>

Looks okay to me. You can probably omit `-Q`. All the steps you've outlined can be done within QIIME 2. Check out the [tutorials](https://library.qiime2.org/books). Nothing should change within QIIME 2 as many of the steps are using generally accepted approaches and wrap the wools in question.

---

<div class="post-metadata">

### Author: ![Phoe](https://forum.qiime2.org/user_avatar/forum.qiime2.org/phoe/32/18852_2.png) [@Phoe](https://forum.qiime2.org/u/Phoe)
#### Post date: [March 15, 2025, 8:08pm UTC](https://forum.qiime2.org/t/cutadapt-primer-removal-questions/32762/5 "2025-03-15T20:08:59Z")

</div>

Thank you so much for you help!

Will look to do all in QIIME2 next time for sure!
