This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Convert exonerate output to EVM gff3 format error

Hello,

Exonerate command:

exonerate --model est2genome --showvulgar no --showalignment no --showquerygff no --showtargetgff yes --percent 80 --softmasktarget yes -q transcripts.fasta -t genome.fasta > exonerate.gff

I used command provided by EVM to convert to gff3:

exonerate_gff_to_alignment_gff3.pl exonerate.gff3

Error:

Error, line has unexpected format: Chr_01       exonerate:est2genome    gene    769406  770204  2591    -       .       gene_id 1 ; sequence CBOT101-C20 ; gene_orientation + at /home/user/bin/EVM_r2012-06-25/EvmUtils/misc/exonerate_gff_to_alignment_gff3.pl line 57, <$fh> line 31.

There is same error when used with exonerate output file provided as example in EVM.

Thanks
SS

evidencemodeler evm gff3 exonerate

Welcome to BioStar! Hmm, yes, the error message is a bit cryptic. Two comments: first, I would check to make sure that exonerate ran correctly and had no error messages; second, if indeed there were no errors when you ran exonerate, I would recommend contacting the EVM developer Brian Haas. In my experience he is very responsive, and since this error message is so vague, he may be in a better position to help you.

Hope this helps!

0 answers

No answers yet.

Log in to answer this question.