Use of several gene predictors

257 views
Skip to first unread message

Juan O A

unread,
Jul 25, 2017, 3:38:00 AM7/25/17
to EVidenceModeler-users
Hi,

I'm trying to include as an input output files from SNAP and Augustus, but I cannot see any other way to include both but to duplicate the --gene_predictions parameter, or to run EVM twice, one for each predictor. As I'm having some errors and I haven't seen any examples of those two possibilities, I would like to know whether any of those are possible or not.

Thanks,

Brian Haas

unread,
Jul 25, 2017, 8:47:29 AM7/25/17
to Juan O A, EVidenceModeler-users
Hi,

You would include all the gene predictions in a single input file but set the 2nd column of the gff3 format to indicate the prediction type.  You then identify the gene prediction types in the weights file:


and that prediction type should match the 2nd column value in the gff3 file.

Be sure to look through the example data to examine formatting, etc.

good luck!

~b


--
You received this message because you are subscribed to the Google Groups "EVidenceModeler-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to evidencemodeler-users+unsub...@googlegroups.com.
To post to this group, send email to evidencemodeler-users@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/evidencemodeler-users/348a0299-77d4-4df6-95d4-00bc00b504bc%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.



--
--
Brian J. Haas
The Broad Institute
http://broadinstitute.org/~bhaas

 

Juan O A

unread,
Jul 26, 2017, 5:37:04 AM7/26/17
to EVidenceModeler-users
Hi,

Thank you very much for your help. I'll do that then. 

Regards,


El martes, 25 de julio de 2017, 14:47:29 (UTC+2), Brian Haas escribió:
Hi,

You would include all the gene predictions in a single input file but set the 2nd column of the gff3 format to indicate the prediction type.  You then identify the gene prediction types in the weights file:


and that prediction type should match the 2nd column value in the gff3 file.

Be sure to look through the example data to examine formatting, etc.

good luck!

~b

On Tue, Jul 25, 2017 at 3:38 AM, Juan O A <jocho...@gmail.com> wrote:
Hi,

I'm trying to include as an input output files from SNAP and Augustus, but I cannot see any other way to include both but to duplicate the --gene_predictions parameter, or to run EVM twice, one for each predictor. As I'm having some errors and I haven't seen any examples of those two possibilities, I would like to know whether any of those are possible or not.

Thanks,

--
You received this message because you are subscribed to the Google Groups "EVidenceModeler-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to evidencemodeler-users+unsub...@googlegroups.com.
To post to this group, send email to evidencemo...@googlegroups.com.

Juan O A

unread,
Jul 27, 2017, 9:36:41 AM7/27/17
to EVidenceModeler-users
(I am not sure if the question was sent, so I'll write it again)

Hi again,

There's an error I keep getting in all the partition folders:

error, no chainID in attributes intron_id 2 ; splice_site "TG" of line scaffold1 exonerate:protein2genome:local splice3 33098 33099 . + . intron_id 2 ; splice_site "TG"
 at /home/jotxoteko/soft/EVidenceModeler-1.1.1/EvmUtils/.././evidence_modeler.pl line 1714, <FILE> line 1.
main::parse_evidence_chains('+', 'protein_alignments_exon.gff3') called at /home/jotxoteko/soft/EVidenceModeler-1.1.1/EvmUtils/.././evidence_modeler.pl line 1217
main::load_evidence_data('+', 'protein_alignments_exon.gff3') called at /home/jotxoteko/soft/EVidenceModeler-1.1.1/EvmUtils/.././evidence_modeler.pl line 725
main::process_features('+') called at /home/jotxoteko/soft/EVidenceModeler-1.1.1/EvmUtils/.././evidence_modeler.pl line 385

Always happens with an intron, but I do not really know what chainID info is missing from the column.

Thanks again,

Brian Haas

unread,
Jul 27, 2017, 8:24:16 PM7/27/17
to Juan O A, EVidenceModeler-users
Hi,

best thing to do is to compare your alignment formatting to what's in the example data provided.

best,

~b

--
You received this message because you are subscribed to the Google Groups "EVidenceModeler-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to evidencemodeler-users+unsub...@googlegroups.com.

For more options, visit https://groups.google.com/d/optout.
Reply all
Reply to author
Forward
Message has been deleted
Message has been deleted
0 new messages