# feature-classifier error while using UNITE

**URL:** https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151
**Category:** Technical Support
**Tags:** feature-classifier, unite
**Created:** [February 11, 2022, 4:44pm UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151 "2022-02-11T16:44:02Z")
**Posts on this page:** 7
**Page:** 1

<div class="post-metadata">

### Author: ![April\_Oliver](https://forum.qiime2.org/user_avatar/forum.qiime2.org/april_oliver/32/12053_2.png) [@April\_Oliver](https://forum.qiime2.org/u/April_Oliver)
#### Post date: [February 11, 2022, 4:44pm UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151/1 "2022-02-11T16:44:03Z")

</div>

I'm struggling to assign taxonomy to ITS amplicons. I downloaded the UNITE v8 database and imported the reference sequences and taxonomy using the following:

```auto
qiime tools import \
   --type FeatureData[Sequence] \
   --input-path sh_refs_qiime_ver8_97_10.05.2021.fasta \
   --output-path unite-v8-97-seqs.qza

qiime tools import \
   --type FeatureData[Taxonomy] \
   --input-format HeaderlessTSVTaxonomyFormat \
   --input-path sh_taxonomy_qiime_ver8_97_10.05.2021.txt \
   --output-path unite-v8-97-tax.qza

```

I then used VSEARCH to assign taxonomy:

```auto
qiime feature-classifier classify-consensus-vsearch \
   --i-query dada2_output/dada2_rep_seqs.qza \
   --i-reference-reads ../ReferenceSequences/unite-v8-97-seqs.qza \
   --i-reference-taxonomy ../ReferenceSequences/unite-v8-97-tax.qza \
   --o-classification taxa/asv_classification.qza

```

And received this error:

> Plugin error from feature-classifier: 'Identifier 227 was reported in taxonomic search results, but was not present in the reference taxonomy.'

I compared the original FASTA and taxonomy files, but didn't find any discrepancy between the 227th sequence and taxonomic ID. Is this an error in the original database, or something with the way I imported these files as q2 artifacts? I also attempted this with the dynamic clustered database and received the same error.

---

<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: [February 11, 2022, 6:30pm UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151/2 "2022-02-11T18:30:52Z")

</div>

Hello April,

This tutorial about [importing the UNITE database into qiime2](https://john-quensen.com/tutorials/training-the-qiime2-classifier-with-unite-its-reference-sequences/) includes a section about reformatting the `fasta` headers.

```bash
awk '/^>/ {print($0)}; /^[^>]/ {print(toupper($0))}' sh_refs_qiime_ver8_99_04.02.2020_dev.fasta | tr -d ' ' > sh_refs_qiime_ver8_99_04.02.2020_dev_uppercase.fasta

```

Could you try that out and let me know if it helps?

---

<div class="post-metadata">

### Author: ![April\_Oliver](https://forum.qiime2.org/user_avatar/forum.qiime2.org/april_oliver/32/12053_2.png) [@April\_Oliver](https://forum.qiime2.org/u/April_Oliver)
#### Post date: [February 11, 2022, 7:25pm UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151/3 "2022-02-11T19:25:24Z")

</div>

Thanks for the link, Colin! I reformatted the fasta headers, but unfortunately I am getting the same error as before. I also used the developer files as the tutorial pointed out.

---

<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: [February 11, 2022, 8:14pm UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151/4 "2022-02-11T20:14:05Z")

</div>

Well I'm stumped! 🤔

Do you get this same error with other identities (like 99 or dynamic, which should both be better anyway)? What version of Qiime are you using?

---

<div class="post-metadata">

### Author: ![April\_Oliver](https://forum.qiime2.org/user_avatar/forum.qiime2.org/april_oliver/32/12053_2.png) [@April\_Oliver](https://forum.qiime2.org/u/April_Oliver)
#### Post date: [February 11, 2022, 8:53pm UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151/5 "2022-02-11T20:53:09Z")

</div>

I am using the latest Qiime 2 v. 2021.11. Yes, I did try using the 99% and dynamic identity set and got the same error.

Strangely enough, I was able to train a classifier with reference sequences and taxonomy, and then used that to assign taxonomy. I cannot fathom why VSEARCH won't work when training a classifier will.

Luckily this reference dataset is somewhat small compared to Silva, so I can easily use classify-sklearn where normally I would opt to use VSEARCH due to my hardware limitations.

---

<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: [February 11, 2022, 9:04pm UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151/6 "2022-02-11T21:04:19Z")

</div>

> [@April\_Oliver](#):
>
> Luckily this reference dataset is somewhat small compared to Silva, so I can easily use classify-sklearn where normally I would opt to use VSEARCH due to my hardware limitations.

That's great to hear! 👍

> [@April\_Oliver](#):
>
> VSEARCH won't work when training a classifier will

I'll look into this, more!

---

<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: [March 15, 2022, 3:04am UTC](https://forum.qiime2.org/t/feature-classifier-error-while-using-unite/22151/7 "2022-03-15T03:04:54Z")

</div>

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