# Manually Adding Positive Controls to classifier or reference database

**URL:** https://forum.qiime2.org/t/manually-adding-positive-controls-to-classifier-or-reference-database/28733
**Category:** General Discussion
**Tags:** taxonomy, rescript, feature-table
**Created:** [December 27, 2023, 7:48pm UTC](https://forum.qiime2.org/t/manually-adding-positive-controls-to-classifier-or-reference-database/28733 "2023-12-27T19:48:51Z")
**Posts on this page:** 1
**Showing post:** 3

<div class="post-metadata">

### Author: ![alexkrohn](https://forum.qiime2.org/user_avatar/forum.qiime2.org/alexkrohn/32/11149_2.png) [@alexkrohn](https://forum.qiime2.org/u/alexkrohn)
#### Post date: [August 30, 2024, 3:26pm UTC](https://forum.qiime2.org/t/manually-adding-positive-controls-to-classifier-or-reference-database/28733/3 "2024-08-30T15:26:20Z")

</div>

Hi @Nicholas_Bokulich. Finally getting around to testing this.

I've imported a multi-fasta of reference samples sequenced with Sanger sequence using `qiime tools import --type 'FeatureData[Sequence'`.

Here's a screenshot of the first few sequences in the FASTA:

 ![Screenshot 2024-08-30 at 11.17.57 AM](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/e/7/e72c04acb0a0808e7de208e6e258ee8c09953364.png)

Based on [these instructions](https://forum.qiime2.org/t/issue-importing-tsv-taxonomy-file/1873), I created a headerless taxonomy file with just the accession number from my FASTA. As with NCBI, I assumed Qiime would consider everything left of the first space in the header to be the accession number.

 ![Screenshot 2024-08-30 at 11.18.53 AM](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/a/1/a1f3b37cccbe9269d609370e95e9572be8f6aa65.png)

Both the taxonomy and sequences import without a problem, but if I just merge them to my GenBank taxonomy and sequences, I don't really have a way of knowing whether the new taxonomy and sequences that I imported are properly connected.

Is there a way to verify that? Or can you affirm that this should work?

Alternatively, it seems you can use the [feature IDs of the sequences, or the accession number](https://forum.qiime2.org/t/how-do-i-import-a-taxonomy-table-what-format-is-required/22160) in the taxonomy file. Is there an easy way to find the feature IDs from a FASTA that I just imported?

Thanks as always for your help!

---

_[View the full topic](https://forum.qiime2.org/t/manually-adding-positive-controls-to-classifier-or-reference-database/28733)._
