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.