# Map reads to KEGG/GO/funct.profile

**URL:** <https://help.galaxyproject.org/t/map-reads-to-kegg-go-funct-profile/11427>\
**Category:** Uncategorized\
**Tags:** reference-annotation\
**Created:** [January 7, 2024, 11:07am UTC](https://help.galaxyproject.org/t/map-reads-to-kegg-go-funct-profile/11427 "2024-01-07T11:07:58Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![vebaev](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/vebaev/32/843_2.png) [@vebaev](https://help.galaxyproject.org/u/vebaev)\
**Post date:** [January 7, 2024, 11:07am UTC](https://help.galaxyproject.org/t/map-reads-to-kegg-go-funct-profile/11427/1 "2024-01-07T11:07:58Z")

</div>

Is there a tool I can use to map reads to make functional profile? (Like kraken is doing for taxonomics)

---

<div class="post-metadata">

**Author:** ![jennaj](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/jennaj/32/27_2.png) [@jennaj](https://help.galaxyproject.org/u/jennaj)\
**Post date:** [January 9, 2024, 7:18pm UTC](https://help.galaxyproject.org/t/map-reads-to-kegg-go-funct-profile/11427/2 "2024-01-09T19:18:00Z")

</div>

Hi @vebaev

The data processing logic goes something like this:

1. Map reads to a genome
2. Compare read genomic coordinates to annotation coordinates to associate known annotation “gene” identifiers with the reads
3. Then: read → gene identifier → functional information

**AnnotateMyIDs** is one tool to do the last step.

---

<div class="post-metadata">

**Author:** ![vebaev](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/vebaev/32/843_2.png) [@vebaev](https://help.galaxyproject.org/u/vebaev)\
**Post date:** [January 9, 2024, 7:30pm UTC](https://help.galaxyproject.org/t/map-reads-to-kegg-go-funct-profile/11427/3 "2024-01-09T19:30:10Z")

</div>

Unfortunately its genome is not available.
