Problem with temporary GFF file about flo HOT 12 CLOSED

I can take a look if you can send me "lifted.gff3".

from flo.

That would be great!
Since the file has a size of around 30MB, I sent you a dropbox link to the email adress stated in your profile.

from flo.

I am seeing the same problem as well... need anymore test data?

from flo.

I think the specific issue with these parents not being defined happens due to them being in the unlifted file

For example I had child features with Parent=SP_0.1_T008586-R3 in lifted.gff3 but then unlifted.gff3 had the actual parent where ID=PKINGS_0.1_T008586-R3

from flo.

from flo.

Ah...I think I remember at one point writing a script to synthesize a parent features for features without parents for something like this...is that what process_gff does?

from flo.

from flo.

Gotcha...I was considering maybe using crossmap, but it looks like it has the same issue

Maybe need to convert from gff to something else, bed12 or similar

from flo.

@cmdcolin:

I am seeing the same problem as well... need anymore test data?

There was a bug. I have made some changes. Can you give it a spin?

@blackFirefly - please see my email

from flo.

@yeban I believe it is working better now, it now gets to the genometools stage, but the genometools ends up crashing

Could maybe ask their team about it, error message isn't easy to interpret

$ rake
mkdir annotations.gff-liftover-target
liftOver -gff annotations.gff run/liftover.chn annotations.gff-liftover-target/lifted.gff3 annotations.gff-liftover-target/unlifted.gff3
Reading liftover chains
Mapping coordinates
WARNING: -gff is not recommended.
Use 'ldHgGene -out=<file.gp>' and then 'liftOver -genePred <file.gp>'
/home/me/flo/gff_recover.rb annotations.gff-liftover-target/lifted.gff3 | gt gff3 -tidy -sort -addids -retainids - > annotations.gff-liftover-target/annotations.gff-liftover-target.gff3
warning: line 1 in file "-" does not begin with "##gff-version" or "##gvf-version", create "##gff-version 3" line automatically
Assertion failed: (elemidx >= q->front), function gt_queue_remove, file src/core/queue.c, line 135.
This is a bug, please report it at
https://github.com/genometools/genometools/issues
Please make sure you are running the latest release which can be found at
http://genometools.org/pub/
You can check your version number with `gt -version`.
Aborted (core dumped)
/home/me/flo/gff_recover.rb:60:in `write': Broken pipe @ io_write - <STDOUT> (Errno::EPIPE)
        from /home/me/flo/gff_recover.rb:60:in `puts'
        from /home/me/flo/gff_recover.rb:60:in `puts'
        from /home/me/flo/gff_recover.rb:60:in `<main>'
rake aborted!

from flo.

At least one thing that could be suspicious is that there are still lines that exist without parents. If I save the file from

gff_recover.rb annotations.gff-liftover-target/lifted.gff3 > out.gff then out.gff (first feature in file) has an mRNA that references a parent gene that is not in out.gff

from flo.

@blackFirefly's problem was partly flo and partly the gff. The former is now fixed.

@cmdcolin I can't be sure what the problem is without looking at the input / lifted gff. Please could you open a new issue with test data?

from flo.