saturated_rainbow_pct(..., na.value = "white", threshold = 0.4, scale = 1)
Adapted from qiimer::saturated_rainbow, this provides a ggplot2 ale_fill_gradientn that adopts the conventions from that scale. Namely: The lowest-abundance taxa are in dark blue, empty taxa are white, and highly-abundant taxa (> threshold) are in red.