Cross-links in GenomeDiagram
I’ve just finished writing up an example for the Biopython Tutorial of the new GenomeDiagram functionality added in Biopython 1.59. You can now control the start and end points of individual tracks, and you can add cross-links between regions of different tracks, as shown here:
This example attempts a simplified reproduction of Figure 6 in Proux et al. (2002), and shows three related phage genomes one above the other. Different classes of genes have been given different colors, while the strength of the red shaded cross-links indicates the percentage identity of the linked genes. Note there are some minor differences in the GenBank annotation we’ve used and the genes shown in the original figure.
Another motivating use case for this functionality was to produce vector images of whole genome alignments in the style of the Artemis Comparison Tool (ACT) or Mauve. We’ve got a poster printer in the building just crying out to be used for showing whole genome comparison of a dozen bacteria strains!