Dear All, I'm new in the field. I want to identify microbial sequences in metatranscriptomic data deprived of rRNA. For this purpose I selected some available databases with microbial genes and aligned my reads to the databases using DIAMOND. These databases comprise of files with sequences in fasta format and annotations in .txt or .tsv format. How can I add the information about biological functions of the identified genes from annotation files to my aligned seqs?