I have results of differential expressed genes from a bacterial transcriptome generated through Deseq2.
The reference genome deposited in ENA (European Nucleotide Archive) is functionally annotated with GO terms.
Now I want to know if there is an specific functional category that is being activated in the test condition.
For this purpose, firstly I want to relate the locus_tag of the upregulated genes (I have a list of locus-tags) to the GO terms of the genome in embl file format.
Do you know a script that could perform this task? Or, at least a script that transform the embl file in a table that relates locus-tag and GO annotation.
1 answer
I found a great tool for this purpose: Go retriever http://agbase.msstate.edu/cgi-bin/tools/goretriever_select.pl
It allows retrieving GO terms using as input Genbank acession numbers.
With these data I can utilize Bingo of Cytoscape for GO enrichment analysis.
Log in to answer this question.
-
Gene ontology with locus tags
written by sea.joson •
-
Number of genes involved in a metabolic pathway as a function of the attributes of the pathway
written by new_bioinfo •
-
How to combine EMBL annotation with sequence file
written by clairebendix •
-
Incompatibility between bbmap and htseq?
written by fhsantanna
-
adding fields to a genbank (translation)
written by Joe
-
How to download GI number and locus tag number for organism?
written by fernan23 •
-
Are ENA and SRA archives in sync?
written by dariober
-
News: Managing and Exploring Next Generation Sequencing Data course at EMBL-EBI, UK
written by rajeshr •
-
Biopython Embl Parser Only Reads One Entry
written by sinanugur •
-
Gff Converter For Use With Mummer/Promer (Microbial Genome)
written by Mat •
