<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;"><br>
<div><br>
<div class="BodyFragment"><font size="2">
<div class="PlainText">Daniel Ence<br>
Graduate Student<br>
Eccles Institute of Human Genetics<br>
University of Utah<br>
15 North 2030 East, Room 2100<br>
Salt Lake City, UT 84112-5330</div>
</font></div>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF42143" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> berenice.benayoun@gmail.com [berenice.benayoun@gmail.com] on behalf of Bérénice Benayoun [benayoun@stanford.edu]<br>
<b>Sent:</b> Friday, June 07, 2013 3:50 PM<br>
<b>To:</b> Daniel Ence<br>
<b>Subject:</b> Re: [maker-devel] Maker and mono-exonic genes ?<br>
</font><br>
</div>
<div></div>
<div>Hi Daniel,
<div><br>
</div>
<div>Thanks for the quick answer !</div>
<div><br>
</div>
<div>I used SNAP, and trained from a hmm model made with the CEGMA output on my genome (240 gene models) plus a first run of maker of 1/3 of the genome. I tried GenemarkES and Augustus, but for some reason they don't run, so I stopped indicating their existence
 to maker.</div>
<div><br>
</div>
<div>Should I do something in particular to train it "better" ? Is there any other predictor that would be worth running ? </div>
<div><br>
</div>
<div>Thanks so much for your help !</div>
<div><br>
</div>
<div>Berenice<br>
<br>
<div class="gmail_quote">2013/6/7 Daniel Ence <span dir="ltr"><<a href="mailto:dence@genetics.utah.edu" target="_blank">dence@genetics.utah.edu</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div style="direction:ltr; font-size:10pt; font-family:Tahoma">Hi Berenice, Thank you for sending that screenshot and the maker_opts.log file. Those are exactly what we need to understand how to expect MAKER to perform. 
<div><br>
In looking at the screenshot, it doesn't look like any of the gene predictors gave a prediction in this region. Uses the predictions from ab-initio tools as a basis for models and considers models that are supported by evidence. It won't by default create a
 model when there isn't a prediction in the region. </div>
<div><br>
Can I ask which gene predictors you used and how they were trained? You might consider training one or more of them on the specific evidence that you expect to support these genes and then rerunning maker with the retrained predictors. </div>
<div><br>
Thanks,</div>
<div>Daniel<br>
<div><br>
<div><font>
<div>Daniel Ence<br>
Graduate Student<br>
Eccles Institute of Human Genetics<br>
University of Utah<br>
15 North 2030 East, Room 2100<br>
Salt Lake City, UT 84112-5330</div>
</font></div>
</div>
<div style="font-size:16px; font-family:Times New Roman">
<hr>
<div style="direction:ltr"><font face="Tahoma" color="#000000"><b>From:</b> maker-devel [<a href="mailto:maker-devel-bounces@yandell-lab.org" target="_blank">maker-devel-bounces@yandell-lab.org</a>] on behalf of Bérénice Benayoun [<a href="mailto:benayoun@stanford.edu" target="_blank">benayoun@stanford.edu</a>]<br>
<b>Sent:</b> Friday, June 07, 2013 11:17 AM<br>
<b>To:</b> <a href="mailto:maker-devel@yandell-lab.org" target="_blank">maker-devel@yandell-lab.org</a><br>
<b>Subject:</b> [maker-devel] Maker and mono-exonic genes ?<br>
</font><br>
</div>
<div>
<div class="h5">
<div></div>
<div>Dear maker developers,
<div><br>
</div>
<div>I am currently annotating a de novo fish genome, and have started looking for genes of interest in particular in Maker's output to verify that it's outputting proper gene sets.</div>
<div><br>
</div>
<div>While many of the genes I look for seem to be correctly annotated by the pipeline, I have noticed that important genes that do have strong evidentiary support but are monoexonic are NOT reported by maker. </div>
<div><br>
</div>
<div>I am attaching a screenshot for the contig that I know should contain the <i>
Foxl2</i> gene (notoriously monoexonic across evolution), and highlighted the corresponding evidence for it.</div>
<div><br>
</div>
<div>Is there any setting I can give to maker to force it to output monoexonic genes ? I already set "single_exon=1" with no success. I attached my config file FYI.</div>
<div><br>
</div>
<div>Thank you so much in advance for your answer !!!</div>
<div><br>
</div>
<div>Best,</div>
<div><br>
</div>
<div>Berenice.</div>
<div>-- <br>
<font>Bérénice A. BENAYOUN, Ph.D.</font><br>
Stanford University/Genetics Department<br>
<i>BRUNET Laboratory</i>, 'Molecular Basis of Longevity and Age Related Diseases'<br>
M312 Alway Building<br>
300, Pasteur Drive<br>
<span>MC 5120 <br>
Stanford, CA 94305-5120</span><br>
USA<br>
Email: <a href="mailto:benayoun@stanford.edu" target="_blank">benayoun@stanford.edu</a><br>
Web: <a href="http://www.stanford.edu/group/brunet/" target="_blank">www.stanford.edu/group/brunet/</a><a href="http://www.stanford.edu/group/brunet/" target="_blank"></a>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<font>Bérénice A. BENAYOUN, Ph.D.</font><br>
Stanford University/Genetics Department<br>
<i>BRUNET Laboratory</i>, 'Molecular Basis of Longevity and Age Related Diseases'<br>
M312 Alway Building<br>
300, Pasteur Drive<br>
<span>MC 5120 <br>
Stanford, CA 94305-5120</span><br>
USA<br>
Email: <a href="mailto:benayoun@stanford.edu" target="_blank">benayoun@stanford.edu</a><br>
Web: <a href="http://www.stanford.edu/group/brunet/" target="_blank">www.stanford.edu/group/brunet/</a><a href="http://www.stanford.edu/group/brunet/" target="_blank"></a>
</div>
</div>
</div>
</div>
</body>
</html>