Comments (8)
Came to post an issue but I think I'm having the same problem noted above, so will just add a concrete example of why this is a problem. I am trying to extract 16S sequences that are annotated in a GenBank file (example). The fact that a gene is the 16S sequence is identified by the product name in the GenBank file,
gene 517900..517988
/locus_tag="SAMN05444282_102329"
rRNA 517900..517988
/locus_tag="SAMN05444282_102329"
/product="16S ribosomal RNA . Bacterial SSU"
However, the product name doesn't make it into the GFF3 file and so it is impossible to select the 16S sequences downstream separately from other rRNA's,
FNQD01000002 GenBank gene 517900 517988 . + . ID=SAMN05444282_102329;locus_tag=SAMN05444282_102329
FNQD01000002 GenBank rRNA 517900 517988 . + . ID=SAMN05444282_102329.rRNA.1;Parent=SAMN05444282_102329
from biocode.
I'll see if I can get this added tonight.
from biocode.
Last night has shifted into today.
from biocode.
@mikemc Is it possible to attach your GBK file so I can test with it, or is it private?
from biocode.
@jorvis The example I gave is from this GenBank file
from biocode.
@mikemc - The current version of the code should fix your issue. The tRNAs now export with anticodon reported and rRNAs with product. I'm not closing this ticket yet, as what @jonathancrabtree reported is actually the reverse conversion, going from GFF3 -> GBK.
from biocode.
Closing. I've now confirmed retention of annotation of tRNAs and rRNAs from source Genbank Flat file, converted to GFF3, then converted back to Genbank.
from biocode.
Great, thanks @jorvis! I haven't had a chance to test yet but sounds like this covers my issue.
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
- Formatting Issue? 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.