Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Renaming coding sequence back to gene number #110

Open
timeresistance1996 opened this issue May 15, 2024 · 3 comments
Open

Renaming coding sequence back to gene number #110

timeresistance1996 opened this issue May 15, 2024 · 3 comments

Comments

@timeresistance1996
Copy link

Hi @gamcil ,

Clinker helps excellent visualization of gene comparisons. However, I encounter a slight issue when using the Identity/Similarity output. How can I rename the coding sequences back to gene numbers?

code
clinker 1_2.gff 1_1.gff -o clinker.out

error:
[14:50:20] WARNING - Could not find parent gene of 1_1. Using coding sequence coordinates instead.
[14:50:20] WARNING - Could not find parent gene of 1_2. Using coding sequence coordinates instead.

output
$ head clinker.out
1_2 vs 1_1

Query Target Identity Similarity
cf905c13-59a1-48fa-93da-e2074e048ddb a9c8c031-1b68-41bb-9303-35cb41674d1f 0.63 0.63
7f2eb39f-2337-40eb-87ad-8f650ca39158 ab28d9d1-b904-4a5f-b67d-24cf95c58a92 1.00 1.00
0e1a809f-e5d7-4d15-bb77-fdc927289453 00075d31-ceb7-4e04-ab93-4cc9b5f00211 1.00 1.00

@gamcil
Copy link
Owner

gamcil commented May 22, 2024

Could you upload one of the GFF files you are using here?

@timeresistance1996
Copy link
Author

Yes, this is header of the 1_2.gff i used which is generated by Prodigal.

1_2_gff.txt

@gamcil
Copy link
Owner

gamcil commented Jun 4, 2024

Looks like the error was clinker not looking for uppercase ID fields. This should be fixed with clinker v0.0.29.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants