# Qiime Analysis pipeline

**URL:** https://forum.qiime2.org/t/qiime-analysis-pipeline/19334
**Category:** User Support
**Tags:** qiime1, qiime2-r
**Created:** [April 22, 2021, 9:16am UTC](https://forum.qiime2.org/t/qiime-analysis-pipeline/19334 "2021-04-22T09:16:57Z")
**Posts on this page:** 3
**Page:** 1

<div class="post-metadata">

### Author: ![mallika](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/m/e5b9ba/32.png) [@mallika](https://forum.qiime2.org/u/mallika)
#### Post date: [April 22, 2021, 9:16am UTC](https://forum.qiime2.org/t/qiime-analysis-pipeline/19334/1 "2021-04-22T09:16:57Z")

</div>

hello  
I'm new to metagenomic pipelines. Currently, I have 500 sample of crohns patients(feces) and 467 samples of healthy patients (feces). All these sample files are in fasta format as follows:

> > GVJ07HB01AWJS1\_cs\_nbp\_rc cs\_nbp=28-337 sample=C\_4022\_01\_S1 rbarcode=GACTCTGA primer=V1-V2 subject=4022 body\_site=stool center=UPENNBL barcode\_mismatch=0 primer\_mismatch=0  
> > ACTAGGCGTTAACACATAGCAAGCGAGGGGACGAGCATCATCAAAGCTTGCTTTGATGGATGGCGACCGGCGGCACGGTGAGTAACACGTATCCAACCTGCCGACAACACTGGGATAGCCTTTCGAAAGAAAGATTAATACCGGATGGCATAATTTTCCCGCATGGGATAATTATTAAAGAATTTCGGTTGTCGATGGGGGATGCGTTCCATTAGGCAGTTGGCGGGGTAACGGCCCACCAAGACAACGATGGATAGGGGTTCTGAGAGGAAGGTCCCCCACATTGGAACTGAGACACGGTCCAAACTCC  
> > GVJ07HB01C9SIH\_cs\_nbp\_rc cs\_nbp=28-381 sample=C\_4022\_01\_S1 rbarcode=GACTCTGA primer=V1-V2 subject=4022 body\_site=stool center=UPENNBL barcode\_mismatch=0 primer\_mismatch=0  
> > TGGTAAGAAGTTTGTAGTCCTGGCGTCAGGATGAACGCTGCGGCGTGCCTAACACATGCAAGTCGAGCGTAAGCGGTTTTAGGAAGTTTTCGGATGGATTAAACTGACTGAGCGGCGGACGGGTGAGTAACGCGTGGGTACCTGCCTCATACAGGGGGATAACAGTTAGAAATGGCTGCTAATACCGCATAAGCACACAGCTTCGCATGGAGCAGTGTGAAAAACTCCGGTGGTATGAGATGGACCCGCGTCTGATTAGCTAGTTGGTAAGGTAACGGCTTACCAAGGCGACGATCAGTAGCCGACCTGAGAGGGTGACCGGCCACATTGGGACTGAGACACGGCCCAAACTCC

After reading up on the forum and the qiime2 website, I came to conclusion that I can combine the 967 fasta files into a single fasta file using qiime1 via virtual machine - add\_qiime\_labels.py, which looks as follows:

> ![image](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/2X/0/02e23afe5cba864ae0d6ee37dc99435a8990ab70.png)

Reading up further on the forum, I'm to analyze the combined file following through steps from [Clustering sequences into OTUs using q2-vsearch — QIIME 2 2021.2.0 documentation](https://docs.qiime2.org/2021.2/tutorials/otu-clustering/)  
And this is where I have questions and please let me know if the above logic behind the pipeline is correct.  
In the OTU clustering tutorial, it essentially clusters based on similarity in denovo step but I believe this tutorial from UCLA ([https://qcb.ucla.edu/wp-content/uploads/sites/14/2017/12/QCB\_W11-Metagenomics-Analysis\_BS\_day2.pdf](https://qcb.ucla.edu/wp-content/uploads/sites/14/2017/12/QCB_W11-Metagenomics-Analysis_BS_day2.pdf) ) asks me to know a few things like which region it is based on before I choose between different OTU methods. I am unaware as to how to find the region of my data (data is from HMP)

Is there a filtering step that I'm missing?

The end product from OTU step via either of three methods would be a rep-seq file and table artifact, which would be piped into the moving tutorial picture at featuretable and featuresummaries step ? Would this pipeline be accurate for a foolproof analysis?

Thanks!

---

<div class="post-metadata">

### Author: ![thermokarst](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/t/8e7dd6/32.png) [@thermokarst](https://forum.qiime2.org/u/thermokarst)
#### Post date: [May 10, 2021, 2:38pm UTC](https://forum.qiime2.org/t/qiime-analysis-pipeline/19334/2 "2021-05-10T14:38:59Z")

</div>

Hi @mallika, I am so sorry for the slow reply, I'm not sure how this post got missed by our moderation team (myself included). The UCLA docs you linked to are for QIIME 1, but much of the reasoning outlined there still holds true in QIIME 2.

> [@mallika](#):
>
> which region it is based on

This is really good information to have, and you will almost certainly want/need this in other aspects of your analysis. As well, for any publications that might come out of this analysis, you'll most likely need to report that information, as well. If you don't know it, please contact your sequencing center. In the case of the HMP, I think you'll need to reference the source publication(s) to learn about the specific samples you might have downloaded.

> [@mallika](#):
>
> The end product from OTU step via either of three methods would be a rep-seq file and table artifact, which would be piped into the moving tutorial picture at featuretable and featuresummaries step ?

Yes!

> [@mallika](#):
>
> Would this pipeline be accurate for a foolproof analysis?

It's not really much of a pipeline, just one initial step along the path of a full analysis. Unfortunately there are many more places that things can go wrong - if you run into problems though, you know where to find us!

:qiime2:

---

<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: [June 10, 2021, 8:39pm UTC](https://forum.qiime2.org/t/qiime-analysis-pipeline/19334/3 "2021-06-10T20:39:14Z")

</div>

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