<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="font-family: Calibri, sans-serif; font-size: 14px; "><blockquote style="margin: 0px 0px 0px 40px; border: none; padding: 0px; "><div><font color="#00007f">So my question is how Maker handles pass-throughs, exactly. Can it</font></div><div><font color="#00007f">reject pass-throughs, or should it always use such models over any other</font></div><div><font color="#00007f">data source? Is there any scenario were it wouldn't?</font></div></blockquote><div style="color: rgb(0, 0, 0); "><br></div></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; ">pred_gff is treated the same as any other ab initio prediction. It is just one among several candidate gene models. The model that is kept is the one with the lowest AED score (lower means better evidence match/support). Any model with no evidence support or AED=1 will be rejected (no evidence support) unless keep_preds=1 is set. There is also another score eAED which takes into account protein reading frame (protein evidence must be in same reading frame as the gene model). If eAED =1 it will also cause models to be rejected.</div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="font-family: Calibri, sans-serif; font-size: 14px; "><blockquote style="margin: 0px 0px 0px 40px; border: none; padding: 0px; "><div><font color="#00007f">I understand that Maker uses some internal scoring system to estimate</font></div><div><font color="#00007f">the accuracy of an annotation - could that be a reason?</font></div></blockquote></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; ">Possibly. Look at the AED score of the pass-through model in the final MAKER GFF3 to see what the AED score was. If you want to send me GFF3 to look at with a list of regions you are concerned about I can tell you more. Also consider giving PASA results to est_gff as well to bias the scoring algorithm to maintain those models (I.e. Model supports itself, which is reasonable since these are EST derived anyways and not just ab initio predictions). Also the model_gff option will always keep an input model (with or without evidence support) and will only replace it with something else if that something else has a better AED score.</div><blockquote style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; margin: 0px 0px 0px 40px; border: none; padding: 0px; "><div><br></div><div><br></div><div><br></div></blockquote><div style="font-family: Calibri, sans-serif; font-size: 14px; "><blockquote style="margin: 0px 0px 0px 40px; border: none; padding: 0px; "><div><font color="#00007f">And a related question: Is there a comprehensive documentation where I</font></div><div><font color="#00007f">can get more information on the internal decision making process of</font></div><div><font color="#00007f">Maker? Or do I have to dig into the code for that?</font></div></blockquote></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; ">Look at these two papers --></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><div style="font-family: Calibri; font-size: medium; "><font face="Calibri,sans-serif">Holt, C., and Yandell, M. (2011). MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinformatics 12, 491.</font></div></div><div><font face="Calibri,sans-serif">Eilbeck, K., Moore, B., Holt, C., and Yandell, M. (2009). Quantitative measures for the management and comparison of annotated genomes. BMC Bioinformatics 10, 67.</font></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; ">Thanks,</div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; ">Carson</div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; ">On 9/13/13 4:15 AM, "Marc P. Hoeppner" <<a href="mailto:marc.hoeppner@imbim.uu.se">marc.hoeppner@imbim.uu.se</a>> wrote:</div><div style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; "><br></div><blockquote id="MAC_OUTLOOK_ATTRIBUTION_BLOCKQUOTE" style="color: rgb(0, 0, 0); font-family: Calibri, sans-serif; font-size: 14px; border-left-color: rgb(181, 196, 223); border-left-width: 5px; border-left-style: solid; padding: 0px 0px 0px 5px; margin: 0px 0px 0px 5px; "><div>Dear list,</div><div><br></div><div>I have started using Maker to explore its use for a number of genome </div><div>projects we are planning on running. One of the tools we intend on </div><div>incorporating into our pipeline is PASA (Since we will be using Trinity </div><div>etc). The (cleaned) output with predicted gene structures I would like </div><div>to pass to Maker as pass-through annotation (I am optimistic that </div><div>way...) - but I noticed that doing so does not always result in the </div><div>incorporation of the PASA gene model into the final maker annotation </div><div>track. Sometimes it seems to be superseded by an Augustus/Maker model, </div><div>sometimes the region stays empty (even tho a protein alignment is present).</div><div><br></div><div>So my question is how Maker handles pass-throughs, exactly. Can it </div><div>reject pass-throughs, or should it always use such models over any other </div><div>data source? Is there any scenario were it wouldn't?</div><div><br></div><div>I understand that Maker uses some internal scoring system to estimate </div><div>the accuracy of an annotation - could that be a reason? It would be a </div><div>bit odd tho, since a lift-over from chicken (to our bird genome) seems </div><div>to support gene models produced by PASA, yet they are nowhere to be </div><div>found in the final models.</div><div><br></div><div>And a related question: Is there a comprehensive documentation where I </div><div>can get more information on the internal decision making process of </div><div>Maker? Or do I have to dig into the code for that?</div><div><br></div><div>Cheers,</div><div><br></div><div>Marc</div><div><br></div><div>PS I have attached a screenshot of such an example - the green track is </div><div>Maker with proteins + augustus (chicken models) + PASA pass-through of a </div><div>cleaned-up gene structure file. (Orange: Cleaned ORFs directly from PASA </div><div>output, Grey: PASA ORFs without cleaning, Dark red: Maker with proteins </div><div>and trinity transcripts as EST evidence, Black: chicken lift-overs from </div><div>EnsEMBL)</div><div>_______________________________________________</div><div>maker-devel mailing list</div><div><a href="mailto:maker-devel@box290.bluehost.com">maker-devel@box290.bluehost.com</a></div><div><a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a></div><div><br></div></blockquote></body></html>