<div dir="auto">Actually, it's a fungal genome. Although not very typical, almost half of it are repeats. Worth mention that Genemark generates a lot of predictions that overlap LTRs and other complex repeats, something that neither SNAP or Augustus do. Have you seen this before? </div><div class="gmail_extra"><br><div class="gmail_quote">On 14 Oct. 2017 02:56, "Carson Holt" <<a href="mailto:carsonhh@gmail.com">carsonhh@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Both transcript and protein evidence will go into the AED calculation for overlap support. So in both cases the chosen model had better overlap (protein evidence will not count toward the eAED overlap calculation if it is out of frame with the model it is supposed to be supporting). The larger merged model generates a clutering affect on it’s evidence, so it’s evidence set for AED calculation is slightly different than the SNAP and Augustus model would generate. In both cases, I think GeneMark is hurting more than it is helping. You may want to just drop it from the analysis (unless it’s a fungi, I often find GeneMark can have that affect).<div><br></div><div>—Carson<br><div><br></div><div><br><div><blockquote type="cite"><div>On Oct 12, 2017, at 12:09 AM, Xabier Vázquez-Campos <<a href="mailto:xvazquezc@gmail.com" target="_blank">xvazquezc@gmail.com</a>> wrote:</div><br class="m_-1089358979162214852Apple-interchange-newline"><div><div dir="ltr"><div>Hi there,<br><br></div>I was visualising the annotations and I realised that in some cases, what it seems to be a gene is splitted according to one of the gene models, despite that the other 2, est2genome and prot2genome suggest that it isn't the case.<br><div><br><span id="m_-1089358979162214852cid:ii_j8o0g4b70_15f0f021abda3b14"><split-gene.png></span><br><br></div><div>Although the opposite also happens.<br></div><div><br><span id="m_-1089358979162214852cid:ii_j8o0g4bw1_15f0f021abda3b14"><merged-gene.png></span><br>​</div><div>For some reason, the "out of place" model is always (or almost) the one from Genemark.</div><div><br></div><div>How much weight does carry the RNAseq and protein data on this decision (if any)?</div><div>How exactly is the final gene selected?</div><div><br></div><div>Cheers,</div><div>Xabi<br></div><div><div><br>-- <br><div class="m_-1089358979162214852gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Xabier Vázquez-Campos, <i>PhD</i><br><i>Research Associate</i><br>NSW Systems Biology Initiative<br>School of Biotechnology and Biomolecular Sciences<br>
The University of New South Wales<br>Sydney NSW 2052 AUSTRALIA<br></div></div></div></div></div></div></div></div></div></div></div>
</div></div></div>
______________________________<wbr>_________________<br>maker-devel mailing list<br><a href="mailto:maker-devel@box290.bluehost.com" target="_blank">maker-devel@box290.bluehost.<wbr>com</a><br><a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" target="_blank">http://box290.bluehost.com/<wbr>mailman/listinfo/maker-devel_<wbr>yandell-lab.org</a><br></div></blockquote></div><br></div></div></div></blockquote></div></div>