Extract subset (based on GO or KEGG categories) from Trinotate.xls for goseq analysis

Hi, I have assembled unigenes using Trinity pipeline and annotated them using Trinotate pipeline. How can I extract a subset of genes (based on GO or KEGG categories, for example oxidoreductases activity) from Trinotate.xls for goseq analysis?

r rna-seq gene

• 3,294 views • link

updated 9.3 years ago by h.mon

• written 9.3 years ago by sbchua.1990 •

1 answer

Using grep:

grep "someGOterm" trinotate.xls

In case you want to filter for several GO terms, put them into a file, one per line:

grep -f GOterm.file trinotate.xls

• 0 views • link

written 9.3 years ago by h.mon

More posts like this

KEGG Pathway analysis for non-model organism using goseq
written 2.8 years ago by Bushra Bibi

Dear all, I have to perform KEGG Pathway analysis of wheat DEGs (differentially expressed genes). I found goseq in R but wheat is not native …
RNA de novo assembly - blasts - KEGG - GO
written 6.4 years ago by dimitrischat

Hello, I am a phd candidate to bioninformatics and with (almost) 0 guidance. Seeking help here.. I was asked to do a de novo RNA …
GO enrichment analysis graphs after goseq: any easy solutions?
written 6.9 years ago by Raito92

Hello everyone! I've performed a differential gene expression analysis through Trinity, as suggested [here][1]. GO enrichment analysis was also performed through the analyze_diff_expr.pl script. The …
Is there a way to get the name or IDs of genes with the goseq table?
written 7.6 years ago by unawaz •

Hi all, I've performed a GO enrichment analysis using goseq. The output of ``goseq()`` tells you the enriched categories with how many significant genes are …
KEGG enrichment for non-model organism
written 9.1 years ago by sbchua.1990 •

Hi, just want to know if I have obtained KEGG annotation and mapped to pathway using KEGG mapper (http://www.genome.jp/kegg/mapper.html). Due to my lack of skill …
Functional annotation pipeline
written 9.2 years ago by sbchua.1990 •

I have assembled a transcriptome using Trinity and clustering it using CD-hit with 0.95 threshold to reduce redundancy. For the clustered transcripts, I performed blastx …
How to create meaningful statistics and chart using annotation file of eggnog mapper and trinotate.
written 9.3 years ago by sbchua.1990 •

Hi, I have finish run annotation using eggnogmapper. Outputs are .annotations, .hmm_hits and .seed_orthologs files. How can I create meaningful charts showing percentages of each …
How to cluster transcripts assembled using TRINITY into unigene after blast?
written 9.3 years ago by sbchua.1990 •

I have assembled a transcriptome using TRINITY with 100,000 transcripts. I have perform blastx against NR database (-outfmt 5). I want to remove redundancy in …
The way to get the all p-value of functionally annotated categories from DAVID?
written 12.2 years ago by kim •

I have an issue about get the info from DAVID. For clustering analysis, I have about DE genes to put them into DAVID as input. …
De Genes And Isoforms?
written 13.3 years ago by lzsph •

Hey guys, We de novo assembly our transcriptome using Trinity, then aligned and quantified using RSEM's Perl script <code>run_RSEM_align_n_estimate.pl</code> supplied with Trinity. This resulted two …