<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Also, you could probably build these overlap sets on the command line by subsetting the MAKER GFF3 file and then using BedTools intersect for overlap queries.
<div class=""><br class="">
</div>
<div class="">Barry</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Oct 11, 2017, at 10:19 PM, Matt Simenc <<a href="mailto:mcsimenc@gmail.com" class="">mcsimenc@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">Very good, thank you!
<div class=""><br class="">
</div>
<div class="">Matt</div>
</div>
<div class="gmail_extra"><br class="">
<div class="gmail_quote">On Wed, Oct 11, 2017 at 8:22 AM, Carson Holt <span dir="ltr" class="">
<<a href="mailto:carsonhh@gmail.com" target="_blank" class="">carsonhh@gmail.com</a>></span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word" class="">Also look at GAL for building GFF3 feature queries —>
<a href="https://github.com/The-Sequence-Ontology/GAL" target="_blank" class="">https://github.com/The-<wbr class="">Sequence-Ontology/GAL</a><span class="HOEnZb"><font color="#888888" class="">
<div class=""><br class="">
</div>
<div class="">—Carson</div>
</font></span>
<div class="">
<div class="h5">
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""> <br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Oct 11, 2017, at 9:18 AM, Michael Campbell <<a href="mailto:michael.s.campbell1@gmail.com" target="_blank" class="">michael.s.campbell1@gmail.com</a><wbr class="">> wrote:</div>
<br class="m_1718221875160280577Apple-interchange-newline">
<div class="">
<div class="">Hi Matt,<br class="">
<br class="">
I have a hacky way that I’ve done it. It requires running MAKER two more times but they are quicker runs.<br class="">
<br class="">
To identify the genes that have protein support I pass all of the annotation back to MAKER using the model_gff option in the maker_opts.ctl file. Then I pull out all of the protein2genome features from the big MAKER GFF3 file and pass them in using the protein_gff
option. I turn off all repeat masking and run MAKER. It runs fast because it doesn’t have to run any gene finders, align evidence, or repeatmask. In the output any gene with an AED less than 1 has protein support. Then I do the same thing with est2genome lines
from the big GFF3 file and put them in as est_gff. The output of that one gives you genes with EST support. Then the genes with an AED of less than one in both sets have support from protein and EST.<br class="">
<br class="">
Hope this helps,<br class="">
Mike<br class="">
<br class="">
<blockquote type="cite" class="">On Oct 11, 2017, at 10:53 AM, Matt Simenc <<a href="mailto:mcsimenc@gmail.com" target="_blank" class="">mcsimenc@gmail.com</a>> wrote:<br class="">
<br class="">
Hey MAKER people,<br class="">
<br class="">
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.<br class="">
<br class="">
QI summary has:<br class="">
<br class="">
Fraction of exons that overlap an EST alignment<br class="">
Fraction of exons that overlap EST or Protein alignments<br class="">
<br class="">
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.<br class="">
<br class="">
Anyone have a way to do this or have a script to parse the MAKER GFF3 to get this?<br class="">
<br class="">
Thanks!!!<br class="">
Matt Simenc<br class="">
______________________________<wbr class="">_________________<br class="">
maker-devel mailing list<br class="">
<a href="mailto:maker-devel@box290.bluehost.com" target="_blank" class="">maker-devel@box290.bluehost.<wbr class="">com</a><br class="">
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" target="_blank" class="">http://box290.bluehost.com/<wbr class="">mailman/listinfo/maker-devel_<wbr class="">yandell-lab.org</a><br class="">
</blockquote>
<br class="">
<br class="">
______________________________<wbr class="">_________________<br class="">
maker-devel mailing list<br class="">
<a href="mailto:maker-devel@box290.bluehost.com" target="_blank" class="">maker-devel@box290.bluehost.<wbr class="">com</a><br class="">
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" target="_blank" class="">http://box290.bluehost.com/<wbr class="">mailman/listinfo/maker-devel_<wbr class="">yandell-lab.org</a><br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
_______________________________________________<br class="">
maker-devel mailing list<br class="">
<a href="mailto:maker-devel@box290.bluehost.com" class="">maker-devel@box290.bluehost.com</a><br class="">
http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org<br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>