Hi, I'm just started to learn how to analyze RNAseq data. I have done everything up until the differential expression analysis. I used R package "edgeR" which is quite simple to follow as they have examples of case study. Now I have a list of genes that I sorted following their log fold change . My result looks like this .
twd logFC logCPM PValue PValue_fdr
Tb927.4.4670 0.024126445 4.950238554 6.298519468 3.23E-53 9.04E-50
Tb927.8.3630 0.018775214 4.442670862 9.299881422 5.19E-79 2.17E-75
Tb09.160.4310 0.054947744 4.314437112 6.441043897 1.52E-24 8.52E-22
Tb927.8.3650 0.016756885 4.192394821 8.989082116 2.47E-79 2.07E-75
Tb927.8.3560 0.110638812 4.006356514 5.12998029 2.41E-11 2.10E-09
Tb927.8.3620 0.03023286 3.745937149 8.791866075 5.55E-39 6.66E-36
Tb927.6.4970 0.029706571 3.402814705 6.098852272 4.56E-28 3.19E-25
I need to get the possible function for all genes in my list so that I can have a result that looks like this
function
Tb927.4.4670 hypothetical protein, conserved
Tb927.8.3630 folate transporter, putative,expression site-associated gene 10 (ESAG10) protein, putative
Tb09.160.4310 glutamate dehydrogenase (GDH)
Tb927.8.3650 folate transporter, putative,expression site-associated gene 10 (ESAG10) protein, putative
Tb927.8.3560 DNA repair and recombination helicase protein PIF5 (PIF5)
I tried to use DAVID but not all genes were returned and the result list was sorted differently so I cannot relate to my original DE result i.e. which I sorted according to logFC. Is there a simple web tool that allow me to upload my list of gene and return the function according my list of genes (without being sorted differently ) ? The organism that I'm looking at is T.Brucei.
