<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif; "><div>I've released 2.28 on the website. This is one of the bugs that was fixed. It happens under a very specific set of circumstances. You need to run maker with the -a command line flag to get it to recalculate upstream variables after upgrading. Alternatively you can also just give maker your old GFF3 file (make all other options blank exempt for the *_pass= options), and maker will just rebuild it.</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; color:black; BORDER-BOTTOM: medium none; BORDER-LEFT: medium none; PADDING-BOTTOM: 0in; PADDING-LEFT: 0in; PADDING-RIGHT: 0in; BORDER-TOP: #b5c4df 1pt solid; BORDER-RIGHT: medium none; PADDING-TOP: 3pt"><span style="font-weight:bold">From: </span> Hung-Wei Hsu <<a href="mailto:ares711122@gmail.com">ares711122@gmail.com</a>><br><span style="font-weight:bold">Date: </span> Monday, 20 May, 2013 9:19 PM<br><span style="font-weight:bold">To: </span> Carson Holt <<a href="mailto:carsonhh@gmail.com">carsonhh@gmail.com</a>><br><span style="font-weight:bold">Cc: </span> <<a href="mailto:maker-devel@yandell-lab.org">maker-devel@yandell-lab.org</a>><br><span style="font-weight:bold">Subject: </span> Re: [maker-devel] Why are some complete gene predictions not present in the final results?<br></div><div><br></div><div dir="ltr">Thanks a lot for your helps.<div style="">Your suggestions will be greatly helpful for our analysis.</div><div style=""><br></div><div style=""><div><font face="arial,sans-serif"><span style="font-size:14px">I've tried to add EST sequences to improve gene predictions.</span></font></div><div style=""><font face="arial,sans-serif"><span style="font-size:14px">The EST sequences I used were CDS sequences of the same organism.</span></font></div><div><font face="arial,sans-serif"><span style="font-size:14px">But I got an error as below.</span></font></div><div><font face="arial,sans-serif"><span style="font-size:14px"><br></span></font></div><div><font face="arial,sans-serif"><span style="font-size:14px">substr outside of string at .../TranslationMachine.pm line 162</span><br></font></div><div><font face="arial,sans-serif"><span style="font-size:14px">ERROR: Failed while polishig ESTs</span></font></div><div><font face="arial,sans-serif"><span style="font-size:14px">ERROR: Chunk failed at level:2, tier_type:3</span></font></div><div><font face="arial,sans-serif"><span style="font-size:14px"><br></span></font></div><div style=""><font face="arial,sans-serif"><span style="font-size:14px">What's wrong with my analysis?</span></font></div><div style=""><span style="font-size: 14px; font-family: arial, sans-serif; ">The EST sequences I used are wrong?</span><br></div><div style=""><font face="arial,sans-serif"><span style="font-size:14px">Thank you.</span></font></div><div style=""><font face="arial,sans-serif"><span style="font-size:14px"><br></span></font></div><div style=""><font face="arial,sans-serif"><span style="font-size:14px">Hung-Wei</span></font></div><div><font face="arial,sans-serif"><span style="font-size:14px"><br></span></font></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/5/21 Carson Holt <span dir="ltr"><<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="font-size:14px;font-family:Calibri,sans-serif;word-wrap:break-word"><div>On default settings MAKER will only put ab initio predictions that have some sort of evidence support (EST or protein) in the final gene set. The rejected predictions are still in the GFF3 for reference purposes as match/match_part features, but not as gene/mRNA/exon/CDS features. So a lack of evidence might be why it is not there. You can add all rejected models that don't overlap an accepted model by setting keep_preds=1 (this usually brings a lot more into the final gene set than you really want though (lots of false positives). But for some organisms like fungi, which have high gene densities, this approach is relatively safe.</div><div><br></div><div>Alternatively the gene is missing because it overlaps another gene model that was accepted. MAKER won't allow overlapping models on the same strand in eukaryotes. The only way to force that kind of overlap is to give MAKER the reference models in model_gff and not let it call it's own models (then maker is really just aligning evidence and scoring the reference models).</div><div><br></div><div>One final note. If there is no evidence supporting the model, and that is why it is rejected, you can also try adding more evidence to the maker run or you can consider the possibility that the gene model in the reference is not real to being with (i.e. a false positive gene model called during the initial annotation process and not supported by protein or expression data from any source).</div><div><br></div><div>Thanks,</div><div>Carson</div><div><br></div><div><br></div><div><br></div><span><div style="border-right:medium none;padding-right:0in;padding-left:0in;padding-top:3pt;text-align:left;font-size:11pt;border-bottom:medium none;font-family:Calibri;border-top:#b5c4df 1pt solid;padding-bottom:0in;border-left:medium none"><span style="font-weight:bold">From: </span> Hung-Wei Hsu <<a href="mailto:ares711122@gmail.com" target="_blank">ares711122@gmail.com</a>><br><span style="font-weight:bold">Date: </span> Monday, 20 May, 2013 12:16 AM<br><span style="font-weight:bold">To: </span> <<a href="mailto:maker-devel@yandell-lab.org" target="_blank">maker-devel@yandell-lab.org</a>><br><span style="font-weight:bold">Subject: </span> [maker-devel] Why are some complete gene predictions not present in the final results?<br></div><div><div class="h5"><div><br></div><div dir="ltr">Hi MAKER developers,<div><br></div><div>I was exploiting MAKER to perform gene prediction and annotation on my contigs.</div><div>I used Artemis to examine gff and found some CDS with complete structure were absent in the final results.</div><div>They are really predicted and annotated on the ref genome.</div><div>I'm wondering if they were discarded due to overlapping with another CDS.</div><div>How can I preserve these CDS?</div><div>
Thanks a lot in advance.</div><div><br></div><div>Hung-Wei</div></div></div></div>
_______________________________________________
maker-devel mailing list
<a href="mailto:maker-devel@box290.bluehost.com" target="_blank">maker-devel@box290.bluehost.com</a><a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" target="_blank">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a></span></div></blockquote></div><br></div></span></body></html>