<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class="">Most single exon alignments from ESTs or mRNA-seq will actually be spurious in nature. It’s the reason why single exon alignments are ignored by default. Also as Daniel mentioned, gene predictors do not like to call single exon genes. Maker gets around this by supplying hints to the predictor based on alignment evidence which increases the probabilities of the HMM used by the predictor within a given region, but given the spurious nature of single exon nucleotide alignments, you will also need protein alignment support for single exon genes for this to work.</div><div class=""><br class=""></div><div class="">Also the region you gave as an example does not appear to be an open reading frame. In fact almost the entire region that overlaps the Augustus call is labeled as UTR. It looks like a classic spurious alignment. Likely an unmasked repetitive element.</div><div class=""><br class=""></div><div class="">—Carson</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><br class=""><div><blockquote type="cite" class=""><div class="">On Aug 4, 2015, at 10:03 AM, Daniel Ence <<a href="mailto:dence@genetics.utah.edu" class="">dence@genetics.utah.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class="">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Nicola, I don’t think that August will predict single-exon genes in most cases. That’s probably responsible for your results.
<div class=""><br class="">
</div>
<div class="">~Daniel</div>
<div class=""><br class="">
</div>
<div class="">Daniel Ence</div>
<div class="">
<div apple-content-edited="true" class="">Graduate Student<br class="">
Eccles Institute of Human Genetics<br class="">
University of Utah<br class="">
15 North 2030 East, Room 2100<br class="">
Salt Lake City, UT 84112-5330 </div>
<br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Aug 4, 2015, at 8:46 AM, Nicola Palmieri <<a href="mailto:palmierinico@gmail.com" class="">palmierinico@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class=""><span style="font-size:12.8000001907349px" class="">Dear MAKER developers,</span>
<div style="font-size:12.8000001907349px" class=""><br class="">
</div>
<div style="font-size:12.8000001907349px" class="">I am using Maker to annotate a new protozoan species (Cystoisospora suis), I am trying to incorporate TopHat/Cufflinks evidence from RNA-Seq but some genes still do not get annotated.</div>
<div style="font-size:12.8000001907349px" class=""><br class="">
</div>
<div style="font-size:12.8000001907349px" class="">I attached a screenshot from IGV, in which I use separately Augustus, Cufflinks and then various run of Maker incorporating TopHat and /or Cufflinks as evidence. The single exon genes are constantly not predicted
(ie. gene under the label g.38.t1). I have also played around with the single_exon option. I didn't find any similar post in the wiki.</div>
<div style="font-size:12.8000001907349px" class=""><br class="">
</div>
<div style="font-size:12.8000001907349px" class="">I would really appreciate some hints on how to proceed.</div>
<div style="font-size:12.8000001907349px" class=""><br class="">
</div>
<div style="font-size:12.8000001907349px" class="">Kind regards,<br class="">
Nicola</div>
<div class=""><br class="">
</div>
-- <br class="">
<div class="gmail_signature">
<div dir="ltr" class="">Nicola Palmieri
<div class="">Postdoctoral fellow<br class="">
Institut für Parasitologie
<div class="">Vetmeduni Vienna</div>
</div>
</div>
</div>
</div>
<span id="cid:8889F48F-DA7C-4588-917B-0D00CD2C5D53@genetics.utah.edu" class=""><Screen Shot 2015-08-04 at 16.39.26.png></span>_______________________________________________<br class="">
maker-devel mailing list<br class="">
<a href="mailto:maker-devel@box290.bluehost.com" class="">maker-devel@box290.bluehost.com</a><br class="">
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" class="">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a><br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
_______________________________________________<br class="">maker-devel mailing list<br class=""><a href="mailto:maker-devel@box290.bluehost.com" class="">maker-devel@box290.bluehost.com</a><br class="">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org<br class=""></div></blockquote></div><br class=""></body></html>