Comments (4)
Would you be able to upload a couple of files you are having issues with? I've tested mostly with GFFs downloaded from NCBI, but files generated by Augustus should also work.
from clinker.
Here is the output. First I use the gff that comes straight out of Augustus. It complains about
ben@iMac-van-iMac locus_5 % /Users/ben/Library/Python/3.8/bin/clinker -p gff.html *.gff [13:35:11] INFO - Starting clinker [13:35:11] INFO - Parsing GenBank files: ['locus5_A1163.gff', 'locus5_Af293.gff', 'locus5_p20.gff', 'locus5_p21.gff'] Traceback (most recent call last): File "/Users/ben/Library/Python/3.8/bin/clinker", line 10, in <module> sys.exit(main()) File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/main.py", line 208, in main clinker( File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/main.py", line 71, in clinker clusters = parse_files(paths) File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/classes.py", line 200, in parse_files cluster = parse_gff(path) File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/classes.py", line 92, in parse_gff seqid = feature.attributes["ID"][0] File "/Users/ben/Library/Python/3.8/lib/python/site-packages/gffutils/attributes.py", line 63, in __getitem__ v = self._d[k] KeyError: 'ID'
Then I try with using the --gff3=on option, and get a similar error message about uid
ben@iMac-van-iMac locus_5 % /Users/ben/Library/Python/3.8/bin/clinker -p gff.html *.gff3 [13:36:14] INFO - Starting clinker [13:36:14] INFO - Parsing GenBank files: ['locus5_A1163.gff3', 'locus5_Af293.gff3', 'locus5_p20.gff3', 'locus5_p21.gff3'] /Users/ben/Library/Python/3.8/lib/python/site-packages/Bio/Seq.py:2334: BiopythonWarning: Partial codon, len(sequence) not a multiple of three. Explicitly trim the sequence or add trailing N before translation. This may become an error in future. warnings.warn( [13:36:14] INFO - Starting cluster alignments Traceback (most recent call last): File "/Users/ben/Library/Python/3.8/bin/clinker", line 10, in <module> sys.exit(main()) File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/main.py", line 208, in main clinker( File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/main.py", line 81, in clinker globaligner = align.align_clusters(*clusters, cutoff=identity, jobs=jobs) File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/align.py", line 51, in align_clusters aligner.add_clusters(*args) File "/Users/ben/Library/Python/3.8/lib/python/site-packages/clinker/align.py", line 325, in add_clusters self._genes[gene.uid] = gene AttributeError: 'NoneType' object has no attribute 'uid'
I can get clinker to work, using the gff2gbSmall.pl script from the Augustus distribution, but the output does not recognize that they are coming from the same DNA strand, and things they are all separate contigs.
I have attached the .gff files, gff3, gb and fast files if that helps.
from clinker.
Thanks for the files, think I've narrowed down to the issue.
By default, Augustus produces GTF files, which is the source of the first error (see daler/gffutils#65). In each row, the parser used by clinker expects an 'ID' key/value pair, which that file does not have. Perhaps in the future I'll try to add a fallback to try and catch this.
Using the -gff3 option like you did fixes that, but then you hit a clinker bug, where gene coordinates were read in incorrectly from GFF files. I've fixed that in #44 which should land in the next release.
from clinker.
Thanks, very excited to use this!
from clinker.
Related Issues (20)
- Errors in using gff and faa file
- The output result has no information HOT 4
- some issues for -mo --matrix_out parameter HOT 2
- Locus inversion showing incorrect locations
- Possible bug in cluster similarity output (-mo option)
- -gf/--gene_functions not working?
- Problem with -gf and -cm parameteres HOT 3
- try example encouter error
- How to only generate html files? HOT 2
- Clinker issue: Blocks instead of arrows HOT 2
- protein groups and similarity HOT 2
- Conda package errors HOT 3
- Similarity score and plotting HOT 1
- Can't save changes to session HOT 2
- update 0.0.28 HOT 3
- the legend HOT 4
- How to visualize a gene on a large chromosome?
- Use clinker to visualize individual exons rather than whole gene?
- ValueError: not enough values to unpack (expected 2, got 1)
- Error in clinker HTML plotting function
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 clinker.