<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">I looked at two cases in which the model_gff disappeared and they occurred in regions where there are multiple overlapping cufflinks features. One model that I'm looking at right now has overlapping protein2genome and a SNAP feature overlapping it but it was still not included in the output. it could be a problem in MAKER or it could be a problem with my RNA Seq data. I aligned the RNA Seq data using tophat/cufflinks and converted the transcripts.gtf file to gff using cufflinks2gff3 script.<div><br></div><div>Is it better to use RNA Seq feature from tophat or cufflinks?</div><div><br><div><br><div><div>On Oct 1, 2012, at 4:01 PM, Carson Holt <<a href="mailto:carsonhh@gmail.com">carsonhh@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; font-size: 14px; font-family: Calibri, sans-serif; "><div>They can be replaced under two circumstances. </div><div>1. If you provide two model_gff files (comma separated list), in which case MAKER thinks it is merging legacy annotations and will only keep one or the other if models overlap. </div><div>2. If you turn snap, augusutus, genemark, or est2genome on. MAKER sees this as a cue that if these other programs produce a better model, it can replace the current model. If you set map_forward=1, MAKER will conserve the name of the previous model (so models change structure but names are conserved); otherwise, it gets a new name. Sometimes groups like to rename models every time their is a structural change. I think you are supposed to get the Alias attribute set when you don't get names mapped forward though (I can't remember if I added this or just planned on adding the Alias mapping though).</div><div><br></div><div>MAKER should never drop a model_gff model. It can only replace it if something better comes along, but it should not disappear.</div><div><br></div><div>Thanks,</div><div>Carson</div><div><br></div><div><br></div><span id="OLK_SRC_BODY_SECTION"><div style="font-family: Calibri; font-size: 11pt; text-align: left; border-width: 1pt medium medium; border-style: solid none none; padding: 3pt 0in 0in; border-top-color: rgb(181, 196, 223); "><span style="font-weight:bold">From: </span> Michael Thon <<a href="mailto:mike.thon@gmail.com">mike.thon@gmail.com</a>><br><span style="font-weight:bold">Date: </span> Monday, 1 October, 2012 1:53 AM<br><span style="font-weight:bold">To: </span> <<a href="mailto:maker-devel@yandell-lab.org">maker-devel@yandell-lab.org</a>><br><span style="font-weight:bold">Subject: </span> [maker-devel] model_gff question<br></div><div><br></div><div><meta http-equiv="Content-Type" content="text/html charset=us-ascii"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Under what circumstances will maker not include a gene model from the model_gff file in its final output? It was my understanding from this post:<div> <a href="https://groups.google.com/d/topic/maker-devel/Y5jSdZ1Olcc/discussion">https://groups.google.com/d/topic/maker-devel/Y5jSdZ1Olcc/discussion</a></div><div><br></div><div>That maker will keep or replace models in model_gff and never remove them. I'm reannotating a fungal genome and in model_gff I'm providing the gene models originally made by the sequencing center. I have 12006 models in the file I specify in model_gff but maker's final annotation has only 10727 models in it. </div><div>-Mike</div><div><br></div></div></div>_______________________________________________
maker-devel mailing list
<a href="mailto:maker-devel@box290.bluehost.com">maker-devel@box290.bluehost.com</a>
<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>
</span></div>
</blockquote></div><br></div></div></body></html>