Comments (12)
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.
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.
Related Issues (20)
- faToTwoBit: cannot execute binary file HOT 5
- gt gff3 error HOT 7
- liftover gtf file HOT 1
- Present as a set of scripts
- gt gff3: error: Parent ... was not defined (via "ID=") HOT 6
- mitochondrial genes HOT 3
- flo compilation error HOT 1
- Workflow manager? HOT 1
- Flo failed at the GenomTools section HOT 3
- flo failed on Large genome HOT 9
- file "-" does not contain 9 tab (\t) separated fields HOT 7
- Syntax issue with rake HOT 2
- Error when blat begins HOT 5
- gff_recover error HOT 4
- undefined method `last' for nil:NilClass (NoMethodError)
- Expecting number line ERROR HOT 2
- cc command not found HOT 1
- GFF format issue with CoGe input
- install.sh fails with a make error for libgenometoools.so HOT 2
- Issue a release & publish on Bioconda? HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from flo.