I was interested in how the LCBs called by SibeliaZ change as I include more genomes, so to start, I ran it with two ~4.5Mb genomes that differed only by a single ~10kb inversion that I introduced. I expected to see something like 3 LCBs since the sequences are otherwise identical, but I was getting around 2100. I can't find anything in the algorithm description in your paper to explain this. There should only be bubbles in the vicinity of the inversion that I introduced, but nowhere else. Is the large number of LCBs for these near-identical genomes expected behavior, and how?
I was interested in how the LCBs called by SibeliaZ change as I include more genomes, so to start, I ran it with two ~4.5Mb genomes that differed only by a single ~10kb inversion that I introduced. I expected to see something like 3 LCBs since the sequences are otherwise identical, but I was getting around 2100. I can't find anything in the algorithm description in your paper to explain this. There should only be bubbles in the vicinity of the inversion that I introduced, but nowhere else. Is the large number of LCBs for these near-identical genomes expected behavior, and how?