[maker-devel] Can maker select a gene model based on #algoritham predicted it

[Barry Moore]

Carson,

How hard would it be to have maker take an option something like 'require_abinits=2' that would instruct maker to promote predictions that overlap with (2, 3 or more) other predictions?  Seems like the maker might have all that info in one place at some point already?

Gowthaman, your contributions to the maker tutorial would be most welcome.  I've got an offline copy of a newer tutorial wiki that is more up to date than the GMOD version.  It's on a server right now that we've got locked behind a firewall, but I'm hoping to move that to a public facing server in the next week and I'd be happy to give you an account on the wiki.

B

On May 30, 2012, at 6:54 AM, Carson Holt wrote:

> It's not an option in exactly the way you are specifying, but there is
> something I usually do for annotation that works well.  I run interproscan
> or rpsblast on the non_overlapping.proteins.fasta file and select just
> those non-overlapping models that have a recognizable protein domain (just
> searching the pfam doamin space is more than sufficient).  Then I provide
> the selected results to model_gff, and provide the previous maker results
> to the maker_gff option with (all reannotation pass options set to 1 and
> all analysis options turned off).  This adds models with at least
> recognizable domains (as even multiple gene predictors can overpredict in
> a similar way).
> 
> Attached is a script to help select predictions and upgrade them to models
> in GFF3 format.  If you have question let me know.
> 
> Thanks,
> Carson
> 
> 
> 
> On 12-05-29 5:54 PM, "Gowthaman Ramasamy"
> <gowthaman.ramasamy at seattlebiomed.org> wrote:
> 
>> Hi Carson,
>> Thanks for all the help during the long weekend, in spite of that long
>> drive. I am still trying to imagine that.
>> 
>> I now have maker to consider our own prediction via pred_gff, and use
>> augustus and gene mark (with our training model). And i was able to use
>> altest and protein evidences. Maker happily picks one gene model when
>> there is a overlap between three different predictions. But, when I look
>> at the gff, it seems like it picks a gene model only when there is an
>> est/protein evidence. It leaves out some genes even though, they are
>> predicted by all three algorithms. Of course, keep_pred=1 helps to keep
>> all the models. This kind of leads to over prediction.
>> 
>> But, I am looking for something in between. And would like to know if
>> that is possible?
>> 1) Pick a gene model if it has an evidence from (est/prot etc...)
>> irrespective of how many algorithms predicted it
>> 2) In the absence of extrinsic evidence (est/prot etc), pick a gene model
>> if that is predicted by at least two algorithms.
>> 
>> Or even simpler:
>> I have ab-initio predictions from three algorithms, Can I output, those
>> genes that is supported by at least two of them. I care less about
>> exactness of gene boundaries.
>> 
>> Thanks,
>> Gowthaman
>> 
>> PS: With my recent attempts, i learned couple things about maker/other
>> associated tools that is not documented in gmod-maker wiki. Is it
>> possible/ok if I add contents to it. I am okay with running it by you
>> before making it public.
> 
> <gff3_preds2models>_______________________________________________
> maker-devel mailing list
> maker-devel at box290.bluehost.com
> http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org

Barry Moore
Research Scientist
Dept. of Human Genetics
University of Utah
Salt Lake City, UT 84112
--------------------------------------------
(801) 585-3543




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://yandell-lab.org/pipermail/maker-devel_yandell-lab.org/attachments/20120601/1a78f472/attachment-0002.html>