I am trying to recreate this plot without a terrible 3d bar chart and an obscure x axis (these are different time points, and it's hard to say when they are).

(from Science 291, No. 5513 (2001): 2606-8, otherwise good paper.)
My first instinct is to do something similar to what they did, with a 2d-line graph and various x-axis labels, using biased stripes for the genotype and then folded stripes to get a black and white split on the front panels, but a few other good questions here say you can't do this.
, ( ), , , . ? ?
: , , ( m n, -), , .

library(tidyverse)
library(cowplot)
data = tribble(
~Timepoint, ~`Ancestral genotype`, ~Mutator, ~`Mean % of auxotrophs`,
100, 'mutS-', 'o', 10.5,
150, 'mutS-', 'o', 16,
220, 'mutS-', 'o', NA,
300, 'mutS-', 'o', 24.5,
100, 'mutS+', 'n', 1,
150, 'mutS+', 'n', NA,
220, 'mutS+', 'n', 1,
300, 'mutS+', 'n', 1,
100, 'mutS+', 'm', 0,
150, 'mutS+', 'm', NA,
220, 'mutS+', 'm', 2,
300, 'mutS+', 'm', 5
)
data <- data %>% mutate(Timepoint = as.character(Timepoint))
data %>% ggplot(aes(x = Timepoint, y = `Mean % of auxotrophs`)) +
geom_col(aes(fill = Mutator), position = 'stack') + facet_grid(~`Ancestral genotype` ) +
guides(fill=FALSE)