Problem processing non-NCBI generated GBK files

[HilbigA]

Problem processing non-NCBI generated GBK files

[HilbigA]

My apologies, I misunderstood that files need to be provided as a list for all filetypes that arent GFF3. While providing a list with the gbk files worked well on a test-subset, I am now running into this error:

.gbk
Traceback (most recent call last):
File "/home/ubuntu/miniconda3/bin/panaroo", line 10, in
sys.exit(main())
File "/home/ubuntu/miniconda3/lib/python3.9/site-packages/panaroo/main.py", line 298, in main
files.append(create_temp_gff3(line[0], None, temp_dir))
File "/home/ubuntu/miniconda3/lib/python3.9/site-packages/panaroo/prokka.py", line 84, in create_temp_gff3
convert_gbk_gff3(gff_file, temp_dir + "temp_gffs/" + prefix + '.gff', True)
File "/home/ubuntu/miniconda3/lib/python3.9/site-packages/panaroo/biocode_convert.py", line 150, in convert_gbk_gff3
gene.add_mRNA(mRNA)
UnboundLocalError: local variable 'gene' referenced before assignment

Any ideas whats causing that? thanks a lot.

[HilbigA]

I solved my issue by chucking gff and gbk files together and using a list of those as input, which worked.

However I am leaving the issue open as it is not clear to me why one subset of the .gbk files was rejected (see 2nd post). Is this connected to the --compliance tag in Prokka?

[gtonkinhill]

Hi,

Sorry for the slow response and I'm glad you got it running.

If it's possible, it would be great if you could share one of the gbk files that caused this error along with the exact command you ran. This would help me a lot in working out what went wrong.

[HilbigA]

Hi, now my apologies for the late reply, I must have missed the notification that you responded. Here is the command I used:
ls .g > list.txt
panaroo -i list.txt -o panaroo_out --clean-mode strict --remove-invalid-genes ->

Github is saying that I mustnt upload gbk format - how do I provide this?

The gbk files that worked fine were produced with the fq2dna pipeline (https://gitlab.pasteur.fr/GIPhy/fq2dna), here however the gff files are not routinely provided with the output but gbk only. The gbk files that did not work were generated by me when running Prokka on contigs that were pre-assembled. From this prokka run I did myself I had the gff files that worked with panaroo, command (where var is my contig.fasta file and pref the basename):
prokka $var --prefix $pref --usegenus --outdir ../ProkkaOut/$pref --kingdom Bacteria --genus Streptococcus --species agalactiae --strain $pref

Thanks a lot for looking into this!

[gtonkinhill]

Hi,

Apologies, I lost track of this issue. I have now added extra error messages to v1.3.4 which indicate that Panaroo requires GBK files that are compliant with Genbank/ENA/DDJB.