How plot GO termns?

[berumeng]

Dear all

I can’t plot the GO termns due to the Trinotate.xls report includes separated GO term columns called "gene_ontology_BLASTX", "gene_ontology_BLASTP" and "gene_ontology_Pfam".

Therefore, the split_GO function doesn't work anymore.

Could the script be modify to plot the GO termns? Or could be fixed?? Any suggestion???

Thank you in advance!!!

[matevzl533]

You can modify the original function to eg. for BLASTP terms:
split_GOp <- function(x, hit = "gene_ontology_BLASTP"){

All the changes are in the first line. I have renamed the function to split_GOp for BLASTP terms and the name of the column that it takes from the Trinotate table (from gene_ontology_blast to gene_ontology_BLASTP).