This is a test version of Biostars. For the public version, visit https://www.biostars.org.
E.C. numbers to GO terms or annotations

I have a list of E.C. numbers (~200) that I would like to 'translate' into either GO terms or, even better, annotation descriptions. Is there a platform/ software on which I can submit that list as input and get their description as output? Thanks!

genome assembly sequence go term e.c. numbers
• 3,329 views • link
updated by Biostar Biostar • written by shibl_a • shibl_a

Please post some example EC numebers.

• 1 views • link
written by cpad0112 cpad0112

here is part of the list:

EC:1.1.1.100 EC:1.1.1.27 EC:1.1.1.29 EC:1.1.1.35 EC:1.1.1.38 EC:1.1.1.39 EC:1.1.1.79 EC:1.1.1.8 EC:1.1.1.81 EC:1.1.1.86 EC:1.1.1.94 EC:1.1.1.95 EC:1.10.2.2

• 1 views • link
written by shibl_a • shibl_a

1 answer

There is a 'conversion' list provided by the EBI, you can find it by linking though from here

• 1 views • link
written by lieven.sterck lieven.sterck

indeed it is, but I pointed to the 'intro' page (which also gives a bit of background info). Will edit my answer

• 1 views • link
written by lieven.sterck lieven.sterck

you are correct that you provided correct link (with background and every thing). Link I posted was for impatient users :). @Lieven

• 1 views • link
written by cpad0112 cpad0112

Thanks! I came across this earlier but I would have to manually search for each EC number on it's own. That would work if I had a few to deal with, but given that my list is quite long I was looking for more of a program that would give me the output we see on the page from EBI.

• 1 views • link
written by shibl_a • shibl_a

Simply use grep or awk or such to select the lines (ECnumbers) you need based on your input list

• 1 views • link
written by lieven.sterck lieven.sterck

This command converts GO2EC, but you can easily modify it to convert EC2GO.

Provide the path to ec2go mapping file:

ec2go = home + "/data/ec2go"

awk 'FNR==NR{go2ec[$NF]=$1; next} { if ($2 in go2ec) { ec=go2ec[$2]; n=split(ec, t, "."); ecshow=substr(ec, 4); for(i=n+1;i<=4;i++) ecshow=ecshow".-"; print $1 " " ecshow; }}'  ec2go  gofile   >  ecfile

https://sourceforge.net/p/fun4me/code/ci/master/tree/fun4me.py

Also, you can use Minpath tool to convert EC numbers to metabolic pathways. It's included in the above package.

• 0 views • link
written by Fatima Fatima

Log in to answer this question.

More posts like this