# Testing ITS1 Pipeline: Finding mock communities or correcting ITSxpress error

**URL:** https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481
**Category:** General Discussion
**Tags:** q2-itsxpress
**Created:** [February 7, 2025, 4:38pm UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481 "2025-02-07T16:38:20Z")
**Posts on this page:** 8
**Page:** 1

<div class="post-metadata">

### Author: ![mbode24](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/m/b782af/32.png) [@mbode24](https://forum.qiime2.org/u/mbode24)
#### Post date: [February 7, 2025, 4:38pm UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/1 "2025-02-07T16:38:20Z")

</div>

I am attempting to verify my pipeline for ITS1 analysis with a mock community from mockrobiota ([mockrobiota/data/mock-25 at master · caporaso-lab/mockrobiota · GitHub](https://github.com/caporaso-lab/mockrobiota/tree/master/data/mock-25)). Most of the pipeline seemed to work well but when I got to ITSxpress it removed all my reads. This seems odd because I didn't have this problem when I ran my own data.

The only thing I could think of is that my data is only the ITS1 region and the mock community contains the whole ITS region. However, when I've tried to modify the --p-region ( ALL, ITS1, ITS2) it still wasn't working properly.

I am wondering if anyone has any ideas how to either troubleshoot this or know of a different mock community specifically with the ITS1 region that I can use.

---

<div class="post-metadata">

### Author: ![salias](https://forum.qiime2.org/user_avatar/forum.qiime2.org/salias/32/18594_2.png) [@salias](https://forum.qiime2.org/u/salias)
#### Post date: [February 12, 2025, 8:52am UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/2 "2025-02-12T08:52:20Z")

</div>

Hi @mbode24 and welcome to the :qiime2: forum!

Could you share the full command you run for ITSxpress?

Best wishes,

Sergio

---

<div class="post-metadata">

### Author: ![mbode24](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/m/b782af/32.png) [@mbode24](https://forum.qiime2.org/u/mbode24)
#### Post date: [February 12, 2025, 2:39pm UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/3 "2025-02-12T14:39:08Z")

</div>

Hi Sergio,

This is the command I ran for ITSxpress. I've tried modifying the --p-region with each option (ALL and ITS2) but still no output.

qiime itsxpress trim-pair-output-unmerged  
--i-per-sample-sequences sequences.qza  
--p-region ITS1  
--p-taxa F  
--p-cluster-id 1  
--p-threads 15  
--o-trimmed trimmed\_exact.qza

Best wishes,  
Michelle

---

<div class="post-metadata">

### Author: ![salias](https://forum.qiime2.org/user_avatar/forum.qiime2.org/salias/32/18594_2.png) [@salias](https://forum.qiime2.org/u/salias)
#### Post date: [February 12, 2025, 3:10pm UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/4 "2025-02-12T15:10:14Z")

</div>

> [@mbode24](#):
>
> I've tried modifying the --p-region with each option (ALL and ITS2)

I think the region in this case is ITS2 (see the paper, [PMC5118932](https://pmc.ncbi.nlm.nih.gov/articles/PMC5118932/)).

> [@mbode24](#):
>
> --p-cluster-id 1

Have you tried playing with this option? I don't think this will solve the issue (I normally use 1 and never had any problem), but it's worth trying. According to the [documentation](https://github.com/USDA-ARS-GBRU/itsxpress), you can use values in the range 0.99-1.0

Best,

Sergio

---

<div class="post-metadata">

### Author: ![mbode24](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/m/b782af/32.png) [@mbode24](https://forum.qiime2.org/u/mbode24)
#### Post date: [February 12, 2025, 4:11pm UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/5 "2025-02-12T16:11:37Z")

</div>

Hi Sergio,

Thanks for the quick response. Yeah when I read it was ITS2 I thought it might work if I switched it in my pipeline but was surprised when it didn't. When I opened the fasta file example it described multiple regions ([mockrobiota/data/mock-25/source/expected-sequences.fasta at master · caporaso-lab/mockrobiota · GitHub](https://github.com/caporaso-lab/mockrobiota/blob/master/data/mock-25/source/expected-sequences.fasta)) so I thought that might be why.

Unfortunately, I did try to modify p-cluster id as well with no luck.

Would you have any idea why this might happen? If you don't think it's fixable, would you have any suggestions for an alternative mock dataset to use?

Thanks for your time,  
Michelle

---

<div class="post-metadata">

### Author: ![salias](https://forum.qiime2.org/user_avatar/forum.qiime2.org/salias/32/18594_2.png) [@salias](https://forum.qiime2.org/u/salias)
#### Post date: [February 13, 2025, 10:13am UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/6 "2025-02-13T10:13:50Z")

</div>

> [@mbode24](#):
>
> Would you have any idea why this might happen?

Not really, seems really strange to me! 🤔

Which steps did you perform prior to ITSxpress? How many reads are you losing in each step?

---

<div class="post-metadata">

### Author: ![mbode24](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/m/b782af/32.png) [@mbode24](https://forum.qiime2.org/u/mbode24)
#### Post date: [February 27, 2025, 7:21pm UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/7 "2025-02-27T19:21:13Z")

</div>

Hi Sergio,

Sorry for the late response, I just got back from a trip.

I followed this tutorial ([Q2-ITSxpress: A tutorial on a QIIME 2 plugin to trim ITS sequences](https://forum.qiime2.org/t/q2-itsxpress-a-tutorial-on-a-qiime-2-plugin-to-trim-its-sequences/5780)) but with initial fastp cleaning as well. There was no unexpected loss in reads until I got to ITSxpress step.

Thanks,  
Michelle

---

<div class="post-metadata">

### Author: ![salias](https://forum.qiime2.org/user_avatar/forum.qiime2.org/salias/32/18594_2.png) [@salias](https://forum.qiime2.org/u/salias)
#### Post date: [March 3, 2025, 10:24am UTC](https://forum.qiime2.org/t/testing-its1-pipeline-finding-mock-communities-or-correcting-itsxpress-error/32481/12 "2025-03-03T10:24:00Z")

</div>

Hi!

Maybe that specific community is reverse complement. I'm not sure if ITSxpress searches in both orientations, so maybe you could try to reverse complement the community before importing to QIIME 2 (with something like `seqtk` or Biopython or something else).

Best,

Sergio
