panaroo icon indicating copy to clipboard operation
panaroo copied to clipboard

Error with panaroo-generate-gffs

Open sarapita opened this issue 3 years ago • 15 comments

I'm trying to update my bakta gffs after panaroo 1.3.0 proccessing with the following command:

panaroo-generate-gffs -o panaroo_out -i *gff3 -t 1 -f 'gff3'

However, I keep getting the following error:

Traceback (most recent call last): File "/services/tools/panaroo/1.3.0/bin/panaroo-generate-gffs", line 10, in sys.exit(main()) File "/services/tools/panaroo/1.3.0/lib/python3.10/site-packages/panaroo/post_run_gff_output.py", line 277, in main parsed_gffs = parse_all_gffs(isolate_names, args.input_files, args.verbose) File "/services/tools/panaroo/1.3.0/lib/python3.10/site-packages/panaroo/post_run_gff_output.py", line 84, in parse_all_gffs input_file = next((s[1] for s in zip(base_input, input_list) if isolate_id == s[0])) StopIteration

Any idea on how I should address it or what might be causing it?

Thanks in advance!

Best, Sara

sarapita avatar Nov 11 '22 07:11 sarapita

Hi,

I am afraid I have only tested the updating of GFFs on Prokka output. I hope to extend this to Bakta GFFs soon so I will leave this issue open as a feature request and update it when the functionality is added.

gtonkinhill avatar Nov 16 '22 02:11 gtonkinhill

Ah, I see thanks.

sarapita avatar Nov 16 '22 06:11 sarapita

Can you run panaroo on bakta annotated genomes?

vonMentzer-Group avatar May 12 '23 17:05 vonMentzer-Group

Yes, it should be possible using the --remove-invalid-genes option. We are hoping to enable panaroo to take advantage of some of the extra features of Bakta soon.

gtonkinhill avatar May 14 '23 17:05 gtonkinhill