[maker-devel] keep_preds option?

Christoph Hahn chrisi.hahni at gmail.com
Fri Aug 31 06:43:21 MDT 2012


Hello maker users and developers,

I am new to gene prediction and I am trying to use maker 2.25 on a newly 
assembled non-model organisms draft genome. Within maker I use genemark, 
SNAP and Augustus. I have a few questions:

1. I was wondering what the keep_preds option means exactly.

I found two slightly different explanations on the option
#Add unsupported gene prediction to final annotation set (maker2.25)
#Add non-overlapping ab-inito gene prediction to final annotation set 
(found on the net - probably older maker version)

As far as I understood keep_preds=0 only retains gene models for which 
the ab initio predictions agree. But how many, all three? two of three?
keep_preds=1 instead keeps all gene models regardless if the different 
programs agree, right?

In my case I get substantial differences in the number of gene models 
found between the two settings, while with =1 I get a number that is 
close to what we would expect. How would you interpret that? What would 
you recommend me to do? Obiously =0 is the saver option.

2. I tried to use EST data of an alternative organism in altest= 
(#EST/cDNA sequence file in fasta format from an alternate organism). 
The organism is quite distantly related, but its the closest I have so I 
thought I d give it a shot. I ran maker twice with identical settigs 
expect in altest and est2genome=0/1. The number of genes predicted is 
identical with both approaches, so I am not sure whether or not the EST 
data was actually used or its just to distantly related. Any easy way to 
assess this?

3. I am running maker in several passes and after each pass I am 
training SNAP using the result of the previous pass. Then for every pass 
I run maker from scratch. Would you recommend to supply the gff of the 
previous pass in "#-----Re-annotation Using MAKER Derived GFF3
maker_gff= #re-annotate genome based on this gff3 file", instead?

Thanks in advance for any thoughts/advice on these things! I appreciate 
your help!

much obliged,
Christoph




More information about the maker-devel mailing list