Formatting Issue? about biocode HOT 2 OPEN

Hmm, in the documentation of the convert_augustus_to_gff3.py script it has an example gene block from Augustus. Can you see how yours compares with that structure? Maybe it has been updated. Regarding the 2nd error, after you've used the -gff3 flag, it seems the features are out of order with some being referenced before they are actually defined. This may be another augustus issue but possibly corrected with this script:

biocode/gff/correct_gff_feature_order.pl

Could you possibly copy/paste a block representing an entire gene's rows here into the ticket?

from biocode.

The first GTF that I get an error on (I put the - in front of the # just to copy paste below):

-# overlap start --------------------------------------------------------------------------------
-# this overlap has 1 different transcripts
-# This transcript jg33035.t1 is derived from g24141.t1 from the input file /scratch/rketchu1/Dovetail_Genome_EM/ANNOTATION/BRAKER3/augustus.Ppri5.gtf
-# It is supported by 0 other predicted genes
-# the core of this joined transcript has priority 2
spez_2 AUGUSTUS gene 44188381 44210443 . + . jg33035
spez_2 AUGUSTUS transcript 44188381 44210443 . + . transcript_id "jg33035.t1"; gene_id "jg33035"
spez_2 AUGUSTUS start_codon 44188381 44188383 . + 0 transcript_id "jg33035.t1"; gene_id "jg33035";
spez_2 AUGUSTUS CDS 44188381 44188722 1 + 0 transcript_id "jg33035.t1"; gene_id "jg33035";
spez_2 AUGUSTUS exon 44188381 44188722 . + . transcript_id "jg33035.t1"; gene_id "jg33035";

The gff3 file looks like this:

spez_2 AUGUSTUS gene 44188381 44210443 . + . ID=jg33035;
spez_2 AUGUSTUS mRNA 44188381 44210443 . + . ID=jg33035.t1;Parent=jg33035;
spez_2 AUGUSTUS start_codon 44188381 44188383 . + 0 ID=jg33035.t1.start1;Parent=jg33035.t1;
spez_2 AUGUSTUS CDS 44188381 44188722 1 + 0 ID=jg33035.t1.CDS1;Parent=jg33035.t1;
spez_2 AUGUSTUS exon 44188381 44188722 . + . ID=jg33035.t1.exon1;Parent=jg33035.t1;
spez_2 AUGUSTUS intron 44188723 44189267 . + . ID=jg33035.t1.intron1;Parent=jg33035.t1;
spez_2 AUGUSTUS CDS 44189268 44189333 1 + 0 ID=jg33035.t1.CDS2;Parent=jg33035.t1;
spez_2 AUGUSTUS exon 44189268 44189333 . + . ID=jg33035.t1.exon2;Parent=jg33035.t1;
spez_2 AUGUSTUS intron 44189334 44189619 . + . ID=jg33035.t1.intron2;Parent=jg33035.t1;
spez_2 AUGUSTUS CDS 44189620 44189671 1 + 0 ID=jg33035.t1.CDS3;Parent=jg33035.t1;

I definitely see some differences but not sure how to resolve the issue.

I am trying to locate correct_gff_feature_order.pl, it looks like my install of biocode is missing this script.

from biocode.