I've assembled the genome from the NGS reads and did the annotation using the PGAP. Now, I want to know if there are any unique genes or proteins present in my sample by comparing the annotations with the reference genome. Is there any automated way to do the deep analysis? I've tried to do it manually by reannotating unique sequences found in Mauve.