# How to generate an ASV table from DADA2 output?

**URL:** https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973
**Category:** User Support
**Created:** [May 28, 2019, 9:20am UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973 "2019-05-28T09:20:39Z")
**Posts on this page:** 10
**Page:** 1

<div class="post-metadata">

### Author: ![Clara](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/c/d26b3c/32.png) [@Clara](https://forum.qiime2.org/u/Clara)
#### Post date: [May 28, 2019, 9:20am UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/1 "2019-05-28T09:20:39Z")

</div>

Hello fellow QIIME2 experts,

I have a few questions regarding ASV table over here. I have generated my OTU table in the taxonomy analysis step and recently get to know very little about ASV (if someone can explain it to me in a simple words that will be much appreciated). Since I am using DADA2 in QIIME2, I wonder I have already generated my ASV outputs in terms of rep-seqs-dada2.qza, stats-dada2.qza and table-dada2.qza, am I right?

If I am right, how can I "use" this ASV table? Or how can I transform them into a readable table form? Is it by converting them from a .qza format into a .biom file? Can they be treated just like an OTU table?

Thanks.

Clara

---

<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: [May 28, 2019, 12:18pm UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/2 "2019-05-28T12:18:16Z")

</div>

> [@Clara](#):
>
> I have generated my OTU table in the taxonomy analysis step and recently get to know very little about ASV (if someone can explain it to me in a simple words that will be much appreciated).

You can think of ASVs as being 100% OTUs that have gone through some sort of denoising process. So instead of being operational taxonomic units clustered as some % threshold as a crude method for reducing noise by collapsing into clusters, they are _amplicon sequence variants_ that recover true signal from the sample. I recommend [this paper](https://www.nature.com/articles/ismej2017119) for a better idea of the transition from ASVs to OTUs.

> [@Clara](#):
>
> Since I am using DADA2 in QIIME2, I wonder I have already generated my ASV outputs in terms of rep-seqs-dada2.qza, stats-dada2.qza and table-dada2.qza, am I right?

Correct!

> [@Clara](#):
>
> how can I “use” this ASV table?

Same exact way as an OTU table. Both ASV and OTU tables are types of "feature tables" ([feature being a generic term since this type of observation matrix can really contain all sorts of information](https://forum.qiime2.org/t/what-is-a-feature-exactly/2201/3)). See the QIIME 2 tutorials for many examples of what one can do with a feature table.

> [@Clara](#):
>
> Or how can I transform them into a readable table form? Is it by converting them from a .qza format into a .biom file?

Correct

> [@Clara](#):
>
> Can they be treated just like an OTU table?

Yes!

Good luck!

---

<div class="post-metadata">

### Author: ![ben](https://forum.qiime2.org/user_avatar/forum.qiime2.org/ben/32/1790_2.png) [@ben](https://forum.qiime2.org/u/ben)
#### Post date: [May 28, 2019, 1:25pm UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/3 "2019-05-28T13:25:16Z")

</div>

Hi Clara, I replied to you directly in DMs, but I want to add one weird quirk I am not sure @Nicholas_Bokulich mentioned.

The taxonomic phylogeny is not tied to the ASV/OTU biom file. You will need to manually add them together. I think there's a couple of issues, but if you're importing into Phyloseq you need to get the taxonomic file with the OTU/ASV biom file and concatenate them according to ASV/OTU ID.

Ben

---

<div class="post-metadata">

### Author: ![lca123](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/l/47e85d/32.png) [@lca123](https://forum.qiime2.org/u/lca123)
#### Post date: [May 29, 2019, 1:32pm UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/4 "2019-05-29T13:32:14Z")

</div>

Hey there,  
I am a bit confused here... I usually get the output from DADA2 (rep-seqs) and classify them against the silva\_132 99% (available here =\> [Archive](https://www.arb-silva.de/download/archive/qiime/)) using qiime feature-classifier classify-consensus-vsearch to get taxonomy assignment for each feature. Am I doing it right?

---

<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: [May 29, 2019, 1:35pm UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/5 "2019-05-29T13:35:17Z")

</div>

> [@lca123](#):
>
> I usually get the output from DADA2 (rep-seqs) and classify them against the silva\_132 99%

You are doing it right, no worries @lca123

---

<div class="post-metadata">

### Author: ![Clara](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/c/d26b3c/32.png) [@Clara](https://forum.qiime2.org/u/Clara)
#### Post date: [May 29, 2019, 9:46pm UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/6 "2019-05-29T21:46:59Z")

</div>

Thanks for all the replies, it is very helpful.  
But I wonder which step makes the OTU table an OTU table? Because let's say we use DADA2 to filter and trim our raw sequences read, the outcome will actually be an ASV table. Next, if we use this ASV output in the feature-classifier step, the generated artifact taxonomy.qza will actually be an ASV table?

---

<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: [May 29, 2019, 11:58pm UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/7 "2019-05-29T23:58:46Z")

</div>

> [@Clara](#):
>
> But I wonder which step makes the OTU table an OTU table?

using one of the OTU clustering methods in the q2-vsearch plugin. I recommend reading the [overview tutorial](https://docs.qiime2.org/2019.4/tutorials/overview/) for a better idea.

> [@Clara](#):
>
> Because let’s say we use DADA2 to filter and trim our raw sequences read, the outcome will actually be an ASV table.

correct

> [@Clara](#):
>
> Next, if we use this ASV output in the feature-classifier step, the generated artifact taxonomy.qza will actually be an ASV table?

No, because the output of a taxonomy classifier is not a feature table at all, it is only a `FeatureData[Taxonomy]` artifact, which is effectively a list of sequence IDs and their taxonomic identities. It does not contain abundance information.

---

<div class="post-metadata">

### Author: ![Clara](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/c/d26b3c/32.png) [@Clara](https://forum.qiime2.org/u/Clara)
#### Post date: [May 31, 2019, 5:29am UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/8 "2019-05-31T05:29:52Z")

</div>

> [@Nicholas\_Bokulich](#):
>
> using one of the OTU clustering methods in the q2-vsearch plugin. I recommend reading the [overview tutorial](https://docs.qiime2.org/2019.4/tutorials/overview/) for a better idea.

Thanks for the link. By looking at the vsearch plugin, I have 2 questions at the moment and answers for it, but need some validation if my analysis is right:

1. What is the different if I use a **qiime vsearch cluster-features-open-reference** and **qiime feature-classifier classify-sklearn**?  
Since both of these functions started with a FeatureData[Sequence], I am thinking is that means that if I use qiime vsearch cluster-features-open-reference (taken from the tutorial),  
qiime vsearch cluster-features-open-reference   
--i-table table.qza   
--i-sequences rep-seqs.qza   
--i-reference-sequences 85\_otus.qza   
--p-perc-identity 0.85   
--o-clustered-table table-or-85.qza   
--o-clustered-sequences rep-seqs-or-85.qza   
--o-new-reference-sequences new-ref-seqs-or-85.qza

**an OTU table will be generated** (since this is an OTU clustering step) and I will then use this rep-seqs-or-85.qza file in the **qiime feature-classifier classify-sklearn** to generate a taxonomy.qzv file and subsequently in a qiime taxa barplot to obtain the final OTU table?

qiime feature-classifier classify-sklearn   
--i-classifier gg-13-8-99-515-806-nb-classifier.qza   
--i-reads rep-seqs.qza   
--o-classification taxonomy.qza

qiime metadata tabulate   
--m-input-file taxonomy.qza   
--o-visualization taxonomy.qzv

qiime taxa barplot   
--i-table table.qza   
--i-taxonomy taxonomy.qza   
--m-metadata-file sample-metadata.tsv   
--o-visualization taxa-bar-plots.qzv

1. If I do not use qiime vsearch cluster-features-open-reference (or de novo or closed reference), and straight away classify using qiime feature-classifier classify-sklearn followed by qiime metadata tabulate and qiime taxa barplot I am actually generating an ASV table?

---

<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: [May 31, 2019, 11:57am UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/9 "2019-05-31T11:57:39Z")

</div>

> [@Clara](#):
>
> Thanks for the link. By looking at the vsearch plugin, I have 2 questions at the moment and answers for it, but need some validation if my analysis is right:

BIG differences! DO not equate the two, and the overview tutorial should make that clear. The first clusters sequences into OTUs to make an OTU table and list of OTUs (sequences). The second assigns taxonomy to those sequences.

> [@Clara](#):
>
> **an OTU table will be generated** (since this is an OTU clustering step) and I will then use this rep-seqs-or-85.qza file in the **qiime feature-classifier classify-sklearn** to generate a taxonomy.qzv file and subsequently in a qiime taxa barplot to obtain the final OTU table?

Yes, you have the sequence correct. But adding taxonomy does not make it a "final" OTU table. The OTU table is an observation matrix of OTUs in each sample, and in QIIME 2 that is final. Taxonomy is always kept separate — this is a type of "feature metadata", information about the features (OTUs) stored inside a feature table (OTU table).

> [@Clara](#):
>
> 1. If I do not use qiime vsearch cluster-features-open-reference (or de novo or closed reference), and straight away classify using qiime feature-classifier classify-sklearn followed by qiime metadata tabulate and qiime taxa barplot I am actually generating an ASV table?

It depends. Are you using dada2/deblur to denoise your data? If yes, then you have an ASV table, and should not need to cluster those ASVs (I discourage it, and there is lots of discussion on this forum for the reasons for/against and I suggest reading those discussions for more info).

---

<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: [July 1, 2019, 5:57pm UTC](https://forum.qiime2.org/t/how-to-generate-an-asv-table-from-dada2-output/9973/10 "2019-07-01T17:57:42Z")

</div>

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