<div dir="ltr">Fantastic. I thought that might have been the problem.<div><br></div><div>Looking forward to 2.28. Thanks!</div><div><br></div><div style>Rob</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Wed, May 22, 2013 at 10:32 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>It looks like the phase was calculated from the wrong strand orientation. I believe I have corrected this now. I'm checking a few more things, but I'll have 2.28 as the latest release likely tomorrow with the cumulative bug fixes since the last release.</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> Rob Syme <<a href="mailto:rob.syme@gmail.com" target="_blank">rob.syme@gmail.com</a>><br><span style="font-weight:bold">Date: </span> Tuesday, 21 May, 2013 1:57 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] Maker-derived CDS GFF3 phase column<br></div><div><div class="h5"><div><br></div><div dir="ltr"><div>Hi all</div><div><br></div><div>By my reading of the GFF3 spec (<a href="http://sequenceontology.org/resources/gff3.html" target="_blank">http://sequenceontology.org/resources/gff3.html</a>), I'm getting gff3 from Maker that has odd data in the phase column.</div>
<div><br></div><div>For example, see some example Maker output at <a href="https://gist.github.com/robsyme/5617399" target="_blank">https://gist.github.com/robsyme/5617399</a></div><div><br></div><div>There are two exons, 5617 <- 5737 and 5793 <- 5953 with phases 0 and 2, respectively. Both exons are in the reverse strand.</div>
<div><br></div><div>From the spec, phase indicates "the number of bases that should be removed from the beginning of this feature to reach the first base of the next codon", and for "reverse strand features, phase is counted from the end field". </div>
<div><br></div>In the case of the 3' exon (5793 <- 5953), the end field (the 5th column) is 5953. <div>The base at the end field is the first base of the translated CDS, so there should be no bases removed "to reach the first base of the next codon". I suggest that this phase should be 0, not 2.</div>
<div><br></div><div>There is an illustration of the feature at <a href="http://i.imgur.com/DKLxnSf.png" target="_blank">http://i.imgur.com/DKLxnSf.png</a>.</div><div><br></div><div>The output gff3 is correct if "the number of bases that should be removed from the beginning of this feature to reach the first base of the next codon" is measured from the 'left-hand' end of this feature (the start field) rather than the end field.</div>
<div><br></div><div>Has anybody else ran into this problem or am I misreading the gff3 spec?</div><div><br></div><div>Rob Syme</div><div>PhD Student</div><div>Curtin University</div><div><br></div><div><br><div><br></div>
<div><br></div></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>