# Possible to get kingdom for each read?

**URL:** <https://forum.biobakery.org/t/possible-to-get-kingdom-for-each-read/1571>\
**Category:** MetaPhlAn\
**Created:** [January 19, 2021, 6:19pm UTC](https://forum.biobakery.org/t/possible-to-get-kingdom-for-each-read/1571 "2021-01-19T18:19:35Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![Sidduppal](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/sidduppal/32/585_2.png) [@Sidduppal](https://forum.biobakery.org/u/Sidduppal)\
**Post date:** [January 19, 2021, 6:19pm UTC](https://forum.biobakery.org/t/possible-to-get-kingdom-for-each-read/1571/1 "2021-01-19T18:19:35Z")

</div>

Hey I was thinking of using MetaPhlA3 to find out which kingdom each pair of reads belong to. The purpose of doing this is to subset my dataset by removing eukaryotic, viral, and archaeal reads and hopefully get a better assembly by the remaining reads.  
MetaPhlA3 is able to generate an intermediate file telling which unique marker each read is hitting to. Is there a way to use that file, reverse engineer it and tell the kingdom that marker belongs to?

Thanks

---

<div class="post-metadata">

**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\
**Post date:** [February 3, 2021, 10:31am UTC](https://forum.biobakery.org/t/possible-to-get-kingdom-for-each-read/1571/2 "2021-02-03T10:31:31Z")

</div>

MetaPhlAn is not designed to assign taxonomy to all the reads in the metagenome, but it reports presence and relative abundance of the detectables species by mapping the metagenome to a set of marker genes, so the resulting alignment reports only the reads that are matching to a marker gene.