<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="">Any hint based predictions rejected for one of the non-AED filters will still be rejected regardless of keep_preds=1.<div class=""><br class=""></div><div class="">However you are right in that any raw predictions that were rejected by other filters may still be kept by using this workflow, since keep_preds=1 essentially rescues all non-overlapping raw predictions that would have been rejected.</div><div class=""><br class=""></div><div class="">—Carson<br class=""><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Feb 17, 2017, at 2:57 PM, Quanwei Zhang <<a href="mailto:qwzhang0601@gmail.com" class="">qwzhang0601@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class=""><div class="">Hello:<br class=""><br class=""></div>I am following the Maker protocol (2014) to rescue rejected gene models. I set keep_preds=1 when I ran Maker, and then used "<a href="http://quality_filter.pl/" class="">quality_filter.pl</a> -d" to get the "default" Maker report. According to the annotation it will "Prints transcripts with an AED <1". <br class=""><br class=""></div>But in a previous post (the link below), it was said "there are models with AED less than 1 that get rejected for other reasons. ... For example if the only evidence is a protein alignment that has deep
overlapping HSPs (extremely low complexity alignment) it will be
filtered out even though AED is not technically equal to 1. Also if the
overlapping protein evidence is in a different reading frame than the
model it is supposed to support then the AED will be less than 1 but
eAED will be 1 (extended AED), and the model will be rejected."<br class=""><br class=""></div><b class="">So I wonder, whether the "default" result obtained by "<a href="http://quality_filter.pl/" class="">quality_filter.pl</a>" are really the same as those achieved by setting keep_preds=0? Will the script also consider other reasons besides "</b><b class="">AED <1"?<br class=""></b><div class=""><br class=""><a href="https://groups.google.com/forum/#!searchin/maker-devel/The$20reason$20there$20are$20differences$20between$20the$20runs$20is$20that$20there$20are$20models$20with$20AED$20less$20%7Csort:relevance/maker-devel/97aNJkT3bgk/W0MEyqEZAAAJ" class="">https://groups.google.com/forum/#!searchin/maker-devel/The$20reason$20there$20are$20differences$20between$20the$20runs$20is$20that$20there$20are$20models$20with$20AED$20less$20%7Csort:relevance/maker-devel/97aNJkT3bgk/W0MEyqEZAAAJ</a> <br class=""><br class=""></div><div class="">Thanks<br class=""><br class=""></div><div class="">Best<br class=""></div><div class="">Quanwei<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=""></div></div></body></html>