# Weird Taxonomic classification 28S rRNA **URL:** https://forum.qiime2.org/t/weird-taxonomic-classification-28s-rrna/18469 **Category:** User Support **Created:** [February 12, 2021, 12:10pm UTC](https://forum.qiime2.org/t/weird-taxonomic-classification-28s-rrna/18469 "2021-02-12T12:10:18Z") **Posts on this page:** 1 **Showing post:** 5
### Author: ![SoilRotifer](https://forum.qiime2.org/user_avatar/forum.qiime2.org/soilrotifer/32/21071_2.png) [@SoilRotifer](https://forum.qiime2.org/u/SoilRotifer) #### Post date: [February 12, 2021, 3:10pm UTC](https://forum.qiime2.org/t/weird-taxonomic-classification-28s-rrna/18469/5 "2021-02-12T15:10:54Z")
Hi @HugoEira, welcome to :qiime2:! I think you are one of the first to inform us about using RESCRIPt for SILVA LSU data! πŸŽ† > [@HugoEira](#): > > Results of taxonomic classification seem to make sense, apart from one weird classification: > d\_\_Eukaryota; p\_\_Eukaryota; c\_\_Eukaryota; o\_\_Eukaryota; f\_\_Eukaryota; g\_\_Eukaryota > > I was expecting in this case that the classification would be just: > d\_Eukaryota If you use your LSU classifier to classify your reads, and observe as the final classification that the upper level taxonomy is propagated downward through all of the ranks, then this means that you scored a hit to a specific sequence in the database that had no information other than 'Eukaryota'. That is, during the initial steps of constructing your database, either through using the pipeline `get-silva-data` or the action `parse-silva-taxonomy`, the default setting is to propagate taxonomy with `--p-rank-propagation`. _This ensures there are no empty ranks in the output for some use cases and tools that do not like empty ranks._ This can be disabled by using the flag `--p-no-rank-propagation` to obtain the `d__Eukaryota` only label. That is, if the resulting classification returns a taxonomy that is not truncated as in the second case, then the classifier is simply returning a specific hit (or hits) to a sequence (or sequences) that all contain the full `d__Eukaryota; ... g__Eukaryota` string. That is, you may want to consider removing any sequences from your reference database, prior to making the classifier, that do not have at least a phylum or other taxonomy information, as they are not particularly helpful. > [@HugoEira](#): > > In the results I get both: > d\_\_Eukaryota; p\_\_Eukaryota; c\_\_Eukaryota; o\_\_Eukaryota; f\_\_Eukaryota; g\_\_Eukaryota > and > d\_Eukaryota So, in these cases you are observing two different classification results. The first is the case that we discussed above. The second, where you only observe `d__Eukaryota`, is due to the fact that the classifier could not identify, beyond the domain level, what your query sequence was. That is, there could have been several equivalent hits to different Eukaryotes within the database (i.e with different taxonomies), and only the lowest common ancestor taxonomy was returned. That is all the lower ranks are cropped except for the upper most rank that it had reasonable confidence. Hopefully this makes sense. πŸ‘¨β€πŸ« Anyway you are on the right track! πŸ’― _Here are some additional tips worth considering:_ > [@HugoEira](#): > > Filtering sequences by length and taxonomy > qiime rescript filter-seqs-length-by-taxon --i-sequences silva-132-LSU-nr99-seqs-cleaned.qza --i-taxonomy silva-132-LSU-nr99-tax.qza --p-labels Archaea Bacteria Eukaryota --p-min-lens 900 1000 1200 --o-filtered-seqs silva-138-ssu-nr99-seqs-filt.qza --o-discarded-seqs silva-138-ssu-nr99-seqs-discard.qza I would recommend playing around with the sequence lengths for this command. The numeric values were tailored to the 16S / 18S rRNA gene data. Since LSU is larger you may want to consider increasing these? I've not benchmarked these, but just something work considering. πŸ€” You can also pick and choose which ranks you'd like in your classifier too. For more details see the tutorial, and this thread: > [@creating a 12-rank SILVA taxonomy with RESCRIPt](https://forum.qiime2.org/t/creating-a-12-rank-silva-taxonomy-with-rescript/17361/2): > > Hi @aeriel.belk, We had recently updated our SILVA taxonomy parser to allow you to select the exact taxonomy ranks you'd like to use via the --p-ranks flag. You can select these ranks for either the parse-silva-taxonomy or get-silva-data actions: That is, you can use this flag like so: qiime rescript get-silva-data \ --p-rank-propagation \ --p-include-species-labels \ --p-ranks domain kingdom phylum class subclass order family genus \ --output-dir ... or qiime re… Also, do not forget to read our warning about `--p-include-species-labels` within the **Species-labels: caveat emptor!** drop menu in the tutorial. -Good luck and keep us posted! -Mike --- _[View the full topic](https://forum.qiime2.org/t/weird-taxonomic-classification-28s-rrna/18469)._