<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>That's good to know :-)</div><div><br></div><div>--Carson</div><div><br></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> Jeremy Semeiks <<a href="mailto:jeremy.semeiks@utsw.edu">jeremy.semeiks@utsw.edu</a>><br><span style="font-weight:bold">Date: </span> Friday, 21 September, 2012 11:39 AM<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] How to preserve human-friendly IDs when reannotating<br></div><div><br></div>For the record: After analyzing these runs, I have confirmed that neither error I described (ie, "DBD::SQLite::db do failed" and "<span><span>Warning: unable to close filehandle DF properly"</span></span>) affects maker's protein output in any way I can detect.<br><br>Thanks,<br>J<br><br><div class="gmail_quote">On Thu, Sep 13, 2012 at 8:11 AM, Carson Holt <span dir="ltr"><<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>></span> wrote:<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>The error --> DBD::SQLite::db do failed: database is locked at /home/jrs/maker-2.26-beta/bin/../lib/GFFDB.pm line 186.</div><div><br></div><div>The location of the specific error you are getting is probably benign. It is a failure to alter the default cache size for the database. The database should already be populated. I'm planning removing the SQLlite database entirely in the future. Perhaps in favor of something like tabix based indexing of the GFF3 file.</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> Jeremy Semeiks <<a href="mailto:jeremy.semeiks@utsw.edu" target="_blank">jeremy.semeiks@utsw.edu</a>><br><span style="font-weight:bold">Date: </span> Tuesday, 11 September, 2012 7:56 PM<div><div class="h5"><br><span style="font-weight:bold">To: </span> Carson Holt <<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>><br><span style="font-weight:bold">Cc: </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> Re: [maker-devel] How to preserve human-friendly IDs when reannotating<br></div></div></div><div><div class="h5"><div><br></div>maker 2.26.<br><br>And I have verified for myself that the three options I mentioned below suffice to preserve human-friendly IDs when reannotating.<br><br>Thanks,<br>J<br><br><div class="gmail_quote">On Tue, Sep 11, 2012 at 3:33 PM, Carson Holt <span dir="ltr"><<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>></span> wrote:<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>Which version of MAKER are you running (maker --version)?</div><div><br></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> Jeremy Semeiks <<a href="mailto:jeremy.semeiks@utsw.edu" target="_blank">jeremy.semeiks@utsw.edu</a>><br><span style="font-weight:bold">Date: </span> Monday, 10 September, 2012 11:02 AM<br><span style="font-weight:bold">To: </span> Carson Holt <<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>><br><span style="font-weight:bold">Cc: </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> Re: [maker-devel] How to preserve human-friendly IDs when reannotating<br></div><div><div><div><br></div>OK, thanks. So if I understand correctly, to preserve human-friendly IDs requires setting just three options: map_forward=1, maker_gff=<my_human_friendly.gff>, and model_pass=1. (Or instead of the last two I could equivalently just set model_gff to a GFF containing only models.)<br><br>A couple new issues came up when I tried to run with these options. I started maker like this:<br><br>/usr/bin/time mpiexec -n 10 maker -q < /dev/null > maker.oe 2>&1<br><br>1. I get a bunch of messages as follows, but with variable line number:<br><br>DBD::SQLite::db do failed: database is locked at /home/jrs/maker-2.26-beta/bin/../lib/GFFDB.pm line 186.<br><br>I saw that this came up in another thread <<a href="https://groups.google.com/forum/?fromgroups=#!topic/maker-devel/TscBgbQfBX4" target="_blank">https://groups.google.com/forum/?fromgroups=#!topic/maker-devel/TscBgbQfBX4</a>>, but I'm not sure it was ever resolved, nor whether it will affect my reannotation results (as I'm not sure what "your GFF3 results will not be integrated" means). This error did not come up the last time I ran maker for reannotation with similar options in a different directory. And both my current directory and my tmp directory are locally mounted, ie not NFS.<br><br>2. Both in this run and in previous runs, I get a lot of lines like this, seemingly at random:<br><br>Warning: unable to close filehandle DF properly.<br><br><br><div class="gmail_quote">On Mon, Sep 10, 2012 at 6:01 AM, Carson Holt <span dir="ltr"><<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>></span> wrote:<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>The map_forward option requires that the pass option for the gene models be turned on. Otherwise you will have to do some spacial overlap test outside of MAKER.</div><div><br></div><div>If you have a new assembly, you can try mapping the old models onto the new assembly using the old transcripts as input to the est= and setting est2genome=1 (nothing else set, i.e no repeat masking etc.). Then there is an undocumented option that is still a little buggy (hence why it is still undocumented). Add the line est_forward=1 to your control files. This tells MAKER to copy names from the ESTs, build the models directly from their alignment, and to do other things to try and make a 1 to 1 match across the genome. You will have to manually check that it is 1 to 1 in the end (as I said still a little buggy and hence undocumented). Use the resulting file as input to the model_gff option on a separate run with map_forward=1 for additional reannotation wil more evidence, etc. where you want to still be able to map names forward. </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> Jeremy Semeiks <<a href="mailto:jeremy.semeiks@utsw.edu" target="_blank">jeremy.semeiks@utsw.edu</a>><br><span style="font-weight:bold">Date: </span> Sunday, 9 September, 2012 3:49 PM<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] How to preserve human-friendly IDs when reannotating<br></div><div><div><div><br></div>Hi all,<br><br>I have sequenced some novel fungal genomes, and I am annotating them with maker-2.26-beta. The entire project is pretty iterative, in the sense that I first get some seemingly-sane annotation sets, then analyze and compare the proteomes biologically, then reannotate when new data comes in or as I learn more about how maker works. Because I have already attached biological meaning to some of my proteins, I would like to retain the same human-friendly IDs across annotations. Eg, if maker suddenly finds 1,000 new proteins on a reannotation run because I turned on keep_preds, then I don't want the transcript formerly known as mymold_09652T0 to become mymold_10698T0 when I run maker_map_ids; I want to keep it named mymold_09652T0.<br><br>So, is there any built-in way to preserve human-friendly IDs, or do I need to write my own script for this? I have tried setting map_forward=1 and maker_gff=<the GFF file output by the previous run of maker_map_ids>, but setting these seems to preserve neither the human-friendly IDs nor even the original IDs. (Eg, protein "genemark-scaffold353-processed-gene-0.9-mRNA-1" changed its name to "genemark-scaffold353-processed-gene-0.6-mRNA-1" when reannotated.) I haven't turned on any of the *_pass options, eg protein_pass; would this be relevant?<br><br>Extra credit question: I am making some mate-pair libraries for these fungi; when I re-assemble, that will completely change my scaffold names. Is there any easy way to preserve human-friendly transcript names in this case? As with the above simpler case, I think it would be pretty easy to transfer 90% of the names just by doing an all-vs-all blastp between two annotation sets and fishing out the best hits, but the remaining 10% might be a headache.<br><br>Thanks,<br>Jeremy<br>Grad student, Grishin lab<br>UT Southwestern, Dallas TX<br><a href="tel:510.385.8959" value="+15103858959" target="_blank">510.385.8959</a><br></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></div></span></div></blockquote></div><br></div></div></span></div></blockquote></div><br></span></body></html>