# ITS analysis for Pacbio generated data sets

**URL:** https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317
**Category:** User Support
**Created:** [January 26, 2023, 4:57pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317 "2023-01-26T16:57:34Z")
**Posts on this page:** 13
**Page:** 1

<div class="post-metadata">

### Author: ![divyaprince321](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/d/13edae/32.png) [@divyaprince321](https://forum.qiime2.org/u/divyaprince321)
#### Post date: [January 26, 2023, 4:57pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/1 "2023-01-26T16:57:34Z")

</div>

Hi all?  
I am for the first time analyzing the fungal metagenome. I have obtained my dataset from Pacbio.  
Now my question is that I need to know what Pipe line I have to follow.  
Secondly, I am not able to understand the data file as both forward and reverse reads are in the same file. So, how can I process my files. Do I need to split the data into forward and reverse reads or no. If so how to split.  
Thanks and regards in Advance

---

<div class="post-metadata">

### Author: ![Keegan-Evans](https://forum.qiime2.org/user_avatar/forum.qiime2.org/keegan-evans/32/10422_2.png) [@Keegan-Evans](https://forum.qiime2.org/u/Keegan-Evans)
#### Post date: [February 1, 2023, 4:58pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/2 "2023-02-01T16:58:14Z")

</div>

@divyaprince321,

Any there is a question of workflows and pipelines come up, there are a couple of great resources on our website that often answer basic questions and help get you oriented to the problem: [QIIME 2 for Experienced Microbiome Researchers](https://docs.qiime2.org/2022.11/tutorials/qiime2-for-experienced-microbiome-researchers/) and [Overview of QIIME 2 plugin workflows](https://docs.qiime2.org/2022.11/tutorials/overview/).

That being said, your case is somewhat complicated by a few things that are a little outside of typical QIIME 2 workflows, namely your PacBio sequencing data and your targeting of ITS.

For quality control and deneoising with PacBio sequencing data, the best tool that we have right now is `denoise-ccs` found in q2-DADA2([DOCS](https://docs.qiime2.org/2022.11/plugins/available/dada2/denoise-ccs/)). As far as I am aware, there are not really "forward" or "reverse" reads, rather the genetic material is sequenced continuously in a loop until the quality scores rise to an acceptable level, so what you are seeing with your raw sequencing data makes a lot of sense, and the DADA2 method does dereplicate them as needed as well. You will want to import your raw sequencing data with the type `SampleData[SequencesWithQuality]` to feed it into q2-DADA2.

After running through q2-DADA2, you should be left with a `FeatureTable[Frequency]` and a `FeatureData[Sequence]` that you can use for any of the typical downstream QIIME 2 workflows. To look at ITS specifically, checkout the [ITS Fungal Analysis tutorial](https://forum.qiime2.org/t/fungal-its-analysis-tutorial/7351) and the [q2-ITSexpress tutorial](https://forum.qiime2.org/t/q2-itsxpress-a-tutorial-on-a-qiime-2-plugin-to-trim-its-sequences/5780).

---

<div class="post-metadata">

### Author: ![divyaprince321](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/d/13edae/32.png) [@divyaprince321](https://forum.qiime2.org/u/divyaprince321)
#### Post date: [February 17, 2023, 8:27am UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/3 "2023-02-17T08:27:08Z")

</div>

Thank You Keegan  
For the timely reply.  
Actually, I have both the forward and the reverse reads in a single file, So I am confusing how to proceed. Do I need to sperate the reads, if So how to do it.  
For your convenience, I tried to attach a file but it exceeded the limits.

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [February 17, 2023, 8:56pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/4 "2023-02-17T20:56:02Z")

</div>



---

<div class="post-metadata">

### Author: ![Keegan-Evans](https://forum.qiime2.org/user_avatar/forum.qiime2.org/keegan-evans/32/10422_2.png) [@Keegan-Evans](https://forum.qiime2.org/u/Keegan-Evans)
#### Post date: [February 27, 2023, 5:01pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/5 "2023-02-27T17:01:15Z")

</div>

@divyaprince321,

Were your sequences read in CLR instead of CCS mode?

---

<div class="post-metadata">

### Author: ![divyaprince321](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/d/13edae/32.png) [@divyaprince321](https://forum.qiime2.org/u/divyaprince321)
#### Post date: [February 27, 2023, 7:40pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/6 "2023-02-27T19:40:15Z")

</div>

Thank you Keegan  
So how to process with these files and which pipeline to follow.  
Thanks in Advance

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [March 6, 2023, 11:18pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/7 "2023-03-06T23:18:48Z")

</div>

Hi @divyaprince321,

Please take a look at @Keegan-Evans response - he is asking whether your sequences were read in CLR or CCS mode. It is important that you read and respond to your posts with any information requested from you in order for us to further assist you.

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [March 6, 2023, 11:18pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/8 "2023-03-06T23:18:51Z")

</div>



---

<div class="post-metadata">

### Author: ![divyaprince321](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/d/13edae/32.png) [@divyaprince321](https://forum.qiime2.org/u/divyaprince321)
#### Post date: [March 7, 2023, 7:37am UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/9 "2023-03-07T07:37:39Z")

</div>

Hi @ lizgehret and Keegan-Evans  
Sorry for inconvenience,  
I am not able to understand in what mode my sequences are.  
For your reference, I have shared one file.  
Thanks and Regards

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [March 10, 2023, 8:31pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/11 "2023-03-10T20:31:09Z")

</div>



---

<div class="post-metadata">

### Author: ![Keegan-Evans](https://forum.qiime2.org/user_avatar/forum.qiime2.org/keegan-evans/32/10422_2.png) [@Keegan-Evans](https://forum.qiime2.org/u/Keegan-Evans)
#### Post date: [March 10, 2023, 9:20pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/14 "2023-03-10T21:20:13Z")

</div>

@divyaprince321,  
It sounds like the next step is to contact your sequencing center and confirm the mode of your sequences. If they were read as CLR, you will not be able to use Q2-DADA2 to perform your denoising 😢 If this is the case, you will need to denoise(if necessary) using a PacBio recommended tool(I do not know what their recommended processing for CLR reads is), after which, you can import your data into QIIME2 following the guidance in the [importing tutorial](https://docs.qiime2.org/2023.2/tutorials/importing/) and proceed with your analysis.

---

<div class="post-metadata">

### Author: ![Keegan-Evans](https://forum.qiime2.org/user_avatar/forum.qiime2.org/keegan-evans/32/10422_2.png) [@Keegan-Evans](https://forum.qiime2.org/u/Keegan-Evans)
#### Post date: [March 10, 2023, 9:20pm UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/15 "2023-03-10T21:20:17Z")

</div>



---

<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: [April 11, 2023, 3:20am UTC](https://forum.qiime2.org/t/its-analysis-for-pacbio-generated-data-sets/25317/16 "2023-04-11T03:20:23Z")

</div>

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