cegma : output report
4
0
Entering edit mode
6.7 years ago
du.ma • 0

Hello,

I am trying to understand the cegma's report call : .completeness_report

I am focussing on the column #Prots (Prots = number of 248 ultra-conserved CEGs present in genome). For exemple I obtained 237 (in partial) which mean that in 248 ultra-conserved CEGs 237 are predicted in my genome.

In the other output files, the number of protein is much more because it contain all of the KOG (not 237 but 458 protein). But I am just interest in the ultra-conserved CEGs so I had filtre all of my files (the reference 248 ids come from the file completeness_cutoff.tbl in cegma/data) and I was expected to generate an output with 237 ids, but I obtained 234 proteins ! How is it possible ?

 

My second question is : What is the .number (e.g. KOG0002.2) after KOG ids in cegma output ?

 

Thanks

cegma alignment kog • 3.1k views
ADD COMMENT
0
Entering edit mode

I had a lot of trouble understanding CEGMA's output, and like in your case, reading their documentation and various pages has not helped. From their completeness report, I was interested in finding out exactly what KOGs were found in my dataset, and where were they found. However, the outputted files did not match in number of KOGs with the completeness report. I ended up giving up... but if someone has a little more knowledge on this I am still interested.

ADD REPLY
0
Entering edit mode
ADD COMMENT
0
Entering edit mode

Thank you for your answer but I already read this page and this is not solving my problem

ADD REPLY
0
Entering edit mode
6.7 years ago
Reema Singh ▴ 160

Hello dubarry,

Question :- What is the .number (e.g. KOG0002.2) after KOG ids in cegma output ?

Answer :- The number after KOG id represent the second region considered by BLAST for KOG0002. You can find the detail description in CEGMA FAQ [http://korflab.ucdavis.edu/Datasets/cegma/faq.html#link6] question 5 [ What do the numerical suffix on KOG IDs represent?]

However, for the first question - Could you please describe how did you filter?

 

 

ADD COMMENT
0
Entering edit mode

In my completness report it's written that I have 237 partial proteins predicted in my genome. But in all of my output I have 438 proteins ( wc -l myFile.cegma.id). So I took the 248 references ids in cegma/data/completeness_cutoff.tbl and I just keep them in my output. I was surprise because I didn't find 237 proteins in common between my output and completeness_cutoff.tbl as predicted, but 234.

ADD REPLY
0
Entering edit mode
6.7 years ago
arnstrm ★ 1.8k

Hi Dubarry,

Please check my answer on the other CEGMA related thread on Biostar.

A: CEGMA report meaning

 

ADD COMMENT
0
Entering edit mode
6.6 years ago
keith ▴ 130

There are two sets of core eukaryotic genes (CEGs), a larger set (458 CEGs) that are designed to be used to help train a gene finder in novel genomes. All of the CEGMA output except the completeness report file refer to this larger set of core genes. 

A subset of the 458 CEGs can be used to assess the completeness of the gene-space of your target genome. These 248 CEGs are taken from the larger set but CEGMA uses slightly different filtering criteria to determine whether these are present. So it is possible for CEGMA to report a CEG being present in the set of 458 CEGs but NOT in the subset of 248 CEGs.

Your original question refers to partial CEGs, these are candidate core genes that exceed a score threshold but which do not exceed a length threshold to be considered 'complete' (this is a somewhat arbitrary threshold…is 95% of a gene complete…how about 85%?). Your genome may contain many partial core genes none of which are complete and so none of which will be present in the other CEGMA output files.

 

ADD COMMENT
0
Entering edit mode

is there any way to get the sequences of the 248 CEGs reported in the completeness report, both partial and full length ones?

ADD REPLY
0
Entering edit mode

Not directly from CEGMA. You can turn on verbose mode to get extra output from the last step which would give you information to use in a custom script to then go back and extract the sequences.

ADD REPLY
0
Entering edit mode

Hi Keith,

Here's my question about CEGMA perform on hg19, could you please take a look at it (I wish I could @ you)? Thanks!

C: CEGMA on human assembly

ADD REPLY

Login before adding your answer.

Traffic: 2089 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6