Comments (2)
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.
Related Issues (20)
- convert_gff3_to_ncbi_tbl HOT 5
- Syntax error on gff.py HOT 4
- Exclude mRNA features in bacterial TBL exports
- Attribute error for update_selected_column9_values.py HOT 1
- Biocode.gff module error HOT 2
- [convert_genbank_to_gff3.py] key_error: locus_tag HOT 5
- AttributeError: type object 'str' has no attribute 'maketrans' HOT 2
- AttributeError: 'Gene' object has no attribute 'add_CDS' HOT 4
- Insert EC numbers into chado database issue HOT 5
- convert_augustus_to_gff3.py error HOT 6
- Conda based install HOT 2
- convert_gff3_to_ncbi_tbl.py HOT 2
- convert_gff_to_ncbi_tbl.py HOT 3
- Incorrect parent features from convert_tRNAScanSE_to_gff3.pl HOT 2
- biocode error HOT 16
- product info not printout in tbl HOT 2
- fasta/fasta_simple_stats.py fails on any file with only one sequence
- [convert_genbank_to_gff3.py] No Locus_tag present in my genbank file HOT 2
- Python upgrade/conversion? 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 biocode.