Possible inconsistencies in pathway completeness percentages.

[avandieren]

for bug reports and errors please report the output of: ./gapseq test
Hi,

This is an amazing tool!

I am doing a comparative metabolic analysis among several species of symbiotic bacteria, and when I look the individual all-Pathways.tbl files, for several of the pathways, the percent completeness based on reactions found is above the 66% threshold, yet the "Prediction" is still labeled as "FALSE." Is this a glitch or could you explain why this may be?

Secondly, for some of the pathways, there will be a completeness percentage given, but no corresponding reactions found. I thought it may have to do with the "vague reactions found", but this result seems to be inconsistent from the rest. Is this a glitch or am I missing something?

Thanks

[Calvin2077]

From my experiences a potential answer to your first question is due to your pathways within your species missing key enzymes or having vague reactions (see below)

A potential answer to your second is 'vague reactions' are "Number of reactions without available sequences" (obtained from GapSeq tutorial pages) More specifically, "In cases in which no sequence data is available for specific reactions, the status of the reactions is set to ‘vague’ and these reactions do not count as missing if they account for less than vagueCutoff of the total reactions of a pathway or subsystem. We used a value of 1/3 for this threshold." (GapSeq publication article)

These are just my (a big user of GapSeq) understanding of the results I am not a developer of this application but I totally wish I was