[maker-devel] Model training with AED=0.7 made all contigs FAILED

Lahcen Campbell lcampbell at ebi.ac.uk
Thu Nov 9 04:13:35 MST 2017


Hi folks,

I would just like some insight into a recent round of MAKER annotation I 
performed and returned back 0 Finished contigs.

The genome is a white fly, which I successfully ran MAKER on initally 
with the first round of "Evidence in", so passing in EST evidence as 
aligned transcript gffs, protein homology evidence etc. The run was 
successful and produced a lot of good quality gene models

Statistics:
            24,613 genes with 49,547 transcripts containing 141130 cds.

Now, I know this count is very high for our species, so in the 2nd round 
(completed running over 1 night due to all contigs failing) I attempted 
to increase the threshold for support, by reducing AED to 0.7 from an 
initial 1. Prior to starting the second round I had trained SNAP on the 
first round results and also ran Augustus separately and  passed this 
via the snaphmm, pred_gff option. Finally I set min protein to be no 
less than 100Aa and set est2genome and prot2genome off to allow for gene 
model refinement.

I checked the run today and all ~8,000 contigs/scaffolds returned as 
FAILED with all having tried to be retried once each.

My initial feeling was, I feared I have just lost my initial set of 
24,613 gene models. I know believe that this won't be the case but Im 
not sure... Can anyone explain what might have happened here and what 
consequences will follow given they all returned as failed ? Have they 
been deleted from the MAKER data store ?

I had capturdD all 1st round MAKER output files (GFF, Fasta files etc) 
before attempting this 2nd round (i.e. 1st round of model training) of 
MAKER .

If I have irrevocably changed the datastore for MAKER and lost those 
genes, might I be able to restore to an earlier point (say back to the 
first round of evidence in gene models) by passing the first MAKER gff 
in as "maker_gff=" / "pred_pass=1" / "model_pass=1" ?

Any advice on this would be much appreciated
Lahcen








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


More information about the maker-devel mailing list