[maker-devel] QI codes insufficient - how to get frac exons with EST only?

[Matt Simenc]

Hey MAKER people,

I would like to make a Venn diagram showing the kinds of evidence
supporting gene models in my MAKER annotation where the left side shows
number of genes with EST support only, the right side shows number of genes
with protein support only, and the intersection shows number of genes with
EST and protein support.

QI summary has:

Fraction of exons that overlap an EST alignment
Fraction of exons that overlap EST or Protein alignments

Please correct me if I'm wrong, because I am interpreting the first to be
fraction of exons that overlap an EST alignment and possibly also a protein
alignment. If that is the case then we can't calculate the number of genes
that overlap only EST or (EST and protein) from the QI information.

Anyone have a way to do this or have a script to parse the MAKER GFF3 to
get this?

Thanks!!!
Matt Simenc
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://yandell-lab.org/pipermail/maker-devel_yandell-lab.org/attachments/20171011/c838a742/attachment-0002.html>