[maker-devel] error: training genemodel with SNAP and GeneMark & run time to generate AUGUTUS species file

[Daniel Ence]

HI Parul, Training augustus does take a long time. Much longer than for the other two predictors that you mentioned. Have you tried using the webAugustus web portal? The team that made augustus run it and can probably help you with trouble-shooting their page for creating training sets: http://bioinf.uni-greifswald.de/webaugustus/training/create

The error that you got regarding genemark is saying that maker can’t find the genemark and probuild executable files. These are specified in the maker_exe.ctl file, not the “opts” file. You need to put valid paths to those executable files in for the given parameters. This is something that is usually specified during installation of MAKER.

Hope that helps,
Daniel


Daniel Ence
Graduate Student
Eccles Institute of Human Genetics
University of Utah
15 North 2030 East, Room 2100
Salt Lake City, UT 84112-5330

On Nov 29, 2016, at 10:13 AM, Kudtarkar, Parul V. <parulk at caltech.edu<mailto:parulk at caltech.edu>> wrote:



Dear Maker developers,

1. We use assembled RNAseq(from same species) and protein evidence(from evolutionary close species) to generate training gene structure(1st iteration, est2genome=1,protein2genome=1 ).

2. This is than used to train abinito gene predictors, SNAP and AUGUSTUS.

3. GeneMarkES( version: GeneMark-ES / ET v.4.32) is used to produce training data-set with the command

gmes_petap.pl --sequence pmin_jelly.fa

4. We would be predicting genes using results from SNAP, Genemark and AUGUSTUS(2nd iteration, est2genome=0, protein2genome=0)

I have couple of questions relating to Genemark and AUGUSTUS

1. AUGUSTUS

We do not have a species file for species file of our interest or evolutionary closer species

following command is used to generate species file

/autoAug.pl --genome=pmin_jelly.fa --species=pminiata --cdna=pmin_transcripts.fa --trainingset=genome.gff3 --singleCPU -v --useexisting
AUGUSTUS is taking too long to compute species file, is there a solution for this issue. Using species file from other organism might generate false positives. Is it advised in such situations to not used AUGUSTUS model?


2. Genemark

I used the gmhmm file generated in the genemark output directory, however I encounter following error


-------------------------

STATUS: Parsing control files...
ERROR: You have failed to provide a value for 'gmhmme3' in the control files.
ERROR: You have failed to provide a value for 'probuild' in the control files.
---------------------
FYI


-----

maker_opts.ctl

#-----Gene Prediction
snaphmm=/home/parul/Pmin_new/maker_snap/pmin1.hmm #SNAP HMM file
gmhmm=/home/parul/Pmin_new/maker_snap/gmhmm.mod #GeneMark HMM file


-----

Using SNAP for training gene model yields over 6000-7000 additional gene. The model has good cumulative AED value.

I was hoping in addition to SNAP, if I could use AUGUSTUS and GeneMark to train the gene model to fuse dispersed models so that the gene count is within the expected range.


Thanks and regards,

Parul


Sent from my iPhone
_______________________________________________
maker-devel mailing list
maker-devel at box290.bluehost.com<mailto:maker-devel at box290.bluehost.com>
http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org

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