eggnog-mapper注释结果COG functional

首先下载COG functional categories代码对应的功能描述信息 fun2003-2014.tab，文件位于COGS数据库中。

这里eggnog-mapper注释结果eggnog-mapper软件emapper.py运行的原始结果，如下图：

eggnog-mapper注释结果

COG functional categories频数统计代码：

```

library(tidyverse)

cogcat <- read_tsv("/home/database/cog/cog2014/fun2003-2014.tab") %>% rename(Code = `# Code`)

emap_cogcat_summary <- read_tsv("rice.emapper.annotations") %>%

  select(Code = `COG cat`) %>% na.omit()%>%

  separate(Code, paste0("X", 1:(max(str_count(.$Code,","))+1), seq = "")) %>%

  gather(key = "X", value = "Code") %>% select(Code) %>% na.omit() %>%

  group_by(Code) %>% summarise(Count  = n()) %>%

  full_join(cogcat) %>% replace_na(replace = list(Count = 0))

```

cogcat 即COG functional categories代码对应的功能描述信息 ，如下

COG functional categories

统计结果如下：

COG functional categories频数统计