To give context to my methods so far, I have run a differential expression analysis of nonmodel organisms using DESeq2. I was able to use the coding sequences (cds) from a genome published on Ensembl as my reference for DESeq2, though the organism is not functionally annotated.
Is there a way I can create a function annotation or discover gene ontology and connect that to my DESeq2 results? I was considering BLASTing the cds fasta, but I am not sure if this is a viable method.
To give an update, my reference was blasted using uniprot's database and the additional information was gathered using biomaRt. It provided enough information to have a predicted annotation with functional support