Open Jungal10 opened 3 years ago
what is the issue? not clear
the small gff fragment that is pasted can't be loaded because it refers to parent features that don't exist and it does seem that flatfile-to-json fails on the contigs containing a pipe symbol in them...not sure if there is anything we can do about that though
This is the error message when trying to run flatfile-to-json:
at /Users/josedias/Documents/PhD_Cologne/projects/PopAmaranth/popamabrowser/bin/../src/perl5/Bio/JBrowse/ConfigurationManager.pm line 7.
Bio::JBrowse::ConfigurationManager::__ANON__("\x{a}GFF3 parse error: some features reference other features tha"...) called at /Users/josedias/Documents/PhD_Cologne/projects/PopAmaranth/popamabrowser/bin/../src/perl5/../../extlib/lib/perl5/Bio/GFF3/LowLevel/Parser.pm line 195
Bio::GFF3::LowLevel::Parser::_buffer_all_under_construction_features(Bio::GFF3::LowLevel::Parser=HASH(0x7f879e42b7c0)) called at /Users/josedias/Documents/PhD_Cologne/projects/PopAmaranth/popamabrowser/bin/../src/perl5/../../extlib/lib/perl5/Bio/GFF3/LowLevel/Parser.pm line 168
Bio::GFF3::LowLevel::Parser::_buffer_items(Bio::GFF3::LowLevel::Parser=HASH(0x7f879e42b7c0)) called at /Users/josedias/Documents/PhD_Cologne/projects/PopAmaranth/popamabrowser/bin/../src/perl5/../../extlib/lib/perl5/Bio/GFF3/LowLevel/Parser.pm line 73
Bio::GFF3::LowLevel::Parser::next_item(Bio::GFF3::LowLevel::Parser=HASH(0x7f879e42b7c0)) called at /Users/josedias/Documents/PhD_Cologne/projects/PopAmaranth/popamabrowser/bin/../src/perl5/Bio/JBrowse/FeatureStream/GFF3_LowLevel.pm line 15
Bio::JBrowse::FeatureStream::GFF3_LowLevel::next_items(Bio::JBrowse::FeatureStream::GFF3_LowLevel=HASH(0x7f879d83c348)) called at /Users/josedias/Documents/PhD_Cologne/projects/PopAmaranth/popamabrowser/bin/../src/perl5/Bio/JBrowse/Cmd/NCFormatter.pm line 52
Bio::JBrowse::Cmd::NCFormatter::_format(Bio::JBrowse::Cmd::FlatFileToJson=HASH(0x7f879d01cea8), "trackConfig", HASH(0x7f879d1e70f8), "featureStream", Bio::JBrowse::FeatureStream::GFF3_LowLevel=HASH(0x7f879d83c348), "featureFilter", CODE(0x7f879d1f3318), "trackLabel", ...) called at /Users/josedias/Documents/PhD_Cologne/projects/PopAmaranth/popamabrowser/bin/../src/perl5/Bio/JBrowse/Cmd/FlatFileToJson.pm line 128
Bio::JBrowse::Cmd::FlatFileToJson::run(Bio::JBrowse::Cmd::FlatFileToJson=HASH(0x7f879d01cea8)) called at bin/flatfile-to-json.pl line
what is the issue? not clear The issue is that I want to have all the attributes I manually added to the gff file (Panther, Pfam, KO (..), but when I load the indexed file, instead of having these attributes, it only adds an attribute called unique with an offset and a number (like the example on the picture with uniqueID offset-1968619510)
the small gff fragment that is pasted can't be loaded because it refers to parent features that don't exist and it does seem that flat-file-to-json fails on the contigs containing a pipe symbol in them...not sure if there is anything we can do about that though
With the original gff file, it works even with the pipes. And I only pasted new items to the already existing column with attributes. I did not add any parents or changed ID's. Not sure how can then parents point to non-existent IDs.
In the alternative, do you any other method to add more attributes to a gff file without eventually breaking it?
Thanks for your help, super fast as always
Can you share the whole GFF? The error I think is that the Parent= attribute in column 9 refers to another feature's ID that is then not in the file
Ahypochondriacus_459_v2.1.gene.gff3.gz
another_attempt.gff.sorted.gff.gz Ahypochondriacus_459_v2.1.annotation_info.txt
The ‘annotation_info’ is the file having extra attributes that I want to add. The Ahypochondriacus_459_v2.1 is the original fully functional gff. And the another_attempt is the file I created that is not working. Thank you
@Jungal10
In what way is the another_attempt file not working?
When I load the another_attempt.gff.sorted.gz it appears to look ok in the browser
If the issue is strictly about the uniqueID appearing in the feature details then that is probably a small issue, but otherwise I am not sure what the issue is?
My track config is just this
[tracks.test]
urlTemplate=another_attempt.gff.sorted.gff.gz
I tried now with the simple configuration like you shown and it worked.
I had it this way before which resulted in not showing any of the new attributes:
[tracks.Gene annotation A hypochondriacus v2 complete]
urlTemplate=another_attempt.gff.sorted.gff.gz
storeClass=JBrowse/Store/SeqFeature/GFF3Tabix
type=CanvasFeatures
category=Annotation
I just had a bad definition of the track, then.
Any idea why also the uniqueID insists on appearing?
Thank you one more time for your help!
Good to know. This branch might fix the uniqueID appearing.... https://github.com/GMOD/jbrowse/tree/uniqueID_remove_gff3tabix
If you have any interest feel free to test...can probably get merged and released later on too
Also not sure, that track definition looks ok so surprising that mine helped at all! Maybe just some weird luck or something? Hopefully didnt just gloss over a bug but the another_attempt file looked like it worked ok
I have a gff file that contains some info on the annotation of my genome. E.g:
I have another file with extra information on the annotation that I parsed and pasted to the 9th column of the gff file. I followed the same pattern, having the feature identifier, equal sign, feature and semi-colon. e.g
"Panther=PTHR10668,PTHR10668:SF80;"
The new gff is looking like this:
When I tried to use
bin/flatfile-to-json.pl
I had several error messages, so I opted to follow the indexed file format to add the annotation.My problem now is that instead of having all these features that I added no the gff file column, I only have a feature stating `"uniqueID offset-1968619510"
`
Can you help me with this issue, please? Thank you