Map reads to KEGG/GO/funct.profile

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

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.

Unfortunately its genome is not available.