# Importing and Classifying already quality filtered, de-noised .fastq sequences

**URL:** <https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954>\
**Category:** User Support\
**Tags:** import\
**Created:** [July 9, 2018, 6:35pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954 "2018-07-09T18:35:44Z")\
**Posts on this page:** 8\
**Page:** 1

<div class="post-metadata">

**Author:** ![Purrsia\_Felidae](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/p/85e7bf/32.png) [@Purrsia\_Felidae](https://forum.qiime2.org/u/Purrsia_Felidae)\
**Post date:** [July 9, 2018, 6:35pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/1 "2018-07-09T18:35:44Z")

</div>

Hi there!!

I have Illumina 2x150bp paired end sequences .fastq sequences. These are not joined. These sequences are already ready to go, they are denoised, quality filtered, etc. I need to by-pass the Dada2/Deblur steps bc this isn't necessary. What I am trying to do is eventually get taxonomic information using the SILVA\_132 database.

In Qiime1 I used to run the following command:  
pick\_open\_reference\_otus.py -i seqs.fna -p SILVA\_123\_db.parameters.txt -r SILVA\_128\_QIIME\_release/rep\_set/rep\_set\_all/97/97\_otus.fasta -o M40A\_silva123\_out

I am trying to run the equivalent in Qiime2

I followed this for import: [Importing data — QIIME 2 2017.12.0 documentation](https://docs.qiime2.org/2017.12/tutorials/importing/#fastq-manifest-formats)  
following the prompts under: “Fastq manifest” formats"

Here is my command:  
qiime tools import --type SampleData[PairedEndSequencesWithQuality] --input-path manifest\_test.txt --output-path demux.test.qza --source-format PairedEndFastqManifestPhred33

This generated a demux.test.qza.

From here, I found this suggestion to by-pass DADA2:

> [@Is there a method to skip denoising and make table.qza and rep-seqs.qza?](https://forum.qiime2.org/t/is-there-a-method-to-skip-denoising-and-make-table-qza-and-rep-seqs-qza/2994):
>
> I am using some mock data that has non-standard quality profiles and I want to just skip the denoising step that wouldn't make sense, but I can't work out how to do this without using dada2 or deblur. Lots of experience with qiime1, first time using qiime2. Thanks

This was the 'general' CLI suggestion:  
qiime vsearch dereplicate-sequences   
--i-sequences demux.qza   
--o-dereplicated-table table   
--o-dereplicated-sequences rep-seqs

But - as you can see, I need a 'table' to pass to the --o-dereplicated-table command as it gives me an 'Error: Missing option: --o-dereplicated-table' if I don't pass this flag. Only way I seem to find to create this table is to go through the DADA2 steps, which is unnecessary for my data set.

**May I ask if there are any other options on how to generate a rep-seqs.qza?**

Next, I wanted to then run the following command after I generate the rep-seqs.qza:

qiime feature-classifier classify-sklearn   
--i-classifier silva-132-99-nb-classifier.qza   
--i-reads rep-seqs.qza   
--o-classification taxonomy.qza

-I found the silva-132-99-nb-classifier.qza here: [Silva 132 classifiers](https://forum.qiime2.org/t/silva-132-classifiers/3698)

From here, I can attempt to move forward.

Many thanks for your advice!!

---

<div class="post-metadata">

**Author:** ![ebolyen](https://forum.qiime2.org/user_avatar/forum.qiime2.org/ebolyen/32/11_2.png) [@ebolyen](https://forum.qiime2.org/u/ebolyen)\
**Post date:** [July 10, 2018, 4:46pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/2 "2018-07-10T16:46:13Z")

</div>



---

<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:** [July 17, 2018, 1:46pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/3 "2018-07-17T13:46:03Z")

</div>



---

<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:** [July 17, 2018, 1:46pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/4 "2018-07-17T13:46:05Z")

</div>



---

<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:** [July 17, 2018, 1:46pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/5 "2018-07-17T13:46:27Z")

</div>

Hey there @Purrsia_Felidae!

> [@Purrsia\_Felidae](#):
>
> I need to by-pass the Dada2/Deblur steps bc this isn’t necessary

I would advise against skipping that level of QA/QC - perhaps this is a good opportunity for you to compare a few methods on one dataset!

> [@Purrsia\_Felidae](#):
>
> But - as you can see, I need a ‘table’ to pass to the --o-dereplicated-table command as it gives me an ‘Error: Missing option: --o-dereplicated-table’ if I don’t pass this flag. Only way I seem to find to create this table is to go through the DADA2 steps, which is unnecessary for my data set.

Looks like you have a typo or a missing value in that command --- the `--o-*` flags are for the _output_ of the command, meaning, the thing the command is generating or creating. In this case, sequences in, feature table and rep seqs out. Make sense? You won't be able to use your demuxed seqs as input to that command the, you could convert these reads to the [`post_split_libraries.py` QIIME 1 format](http://qiime.org/documentation/file_formats.html#post-split-libraries-fasta-file-overview), and then import and dereplicate them.

I would highly encourage you to take some time reviewing the extensive QIIME 2 documentation, perhaps starting with the [Overview tutorial!](https://docs.qiime2.org/2018.6/tutorials/overview/)

---

<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:** [July 17, 2018, 1:46pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/6 "2018-07-17T13:46:41Z")

</div>



---

<div class="post-metadata">

**Author:** ![Purrsia\_Felidae](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/p/85e7bf/32.png) [@Purrsia\_Felidae](https://forum.qiime2.org/u/Purrsia_Felidae)\
**Post date:** [July 20, 2018, 5:21pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/7 "2018-07-20T17:21:52Z")

</div>

Thank-you for your suggestions and insight.

---

<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:** [August 20, 2018, 11:21pm UTC](https://forum.qiime2.org/t/importing-and-classifying-already-quality-filtered-de-noised-fastq-sequences/4954/8 "2018-08-20T23:21:56Z")

</div>

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