Hi all, was hoping for a bit of advice with something I'm struggling with on my current project.

I have 21 Pectobacterium atrosepticum genomes on which I have ran: PROKKA, ROARY, and METABOLIC, I'm looking to analyse my genomes further using bash and python to identify all bacteriocin genes present on the genomes and pull they sequences onto a new file and determining if any are novel.

Can anyone please outline some potential steps I could take to being able to do this / advise on a tool they'd recommend?

Any help is appreciated!