Adición de leyenda clave a un gráfico de múltiples histogramas en R

¿Cómo agrego una leyenda clave a la siguiente gráfica?

Desearía tener una leyenda clave en algún lugar en la esquina superior derecha con dos barras de color horizontales cortas, donde la roja debería decir "La cirugía plástica salió mal" y la azul debería decir "Alemania".

Usé el siguiente código para producir la trama:

bar2 <- read.table("div/ana-mut[...]/barriers-set-2.dat", sep=" ")
bar2val <- c(bar2$V1, bar2$V2)
bar3 <- read.table("div/ana-mut[...]/barriers-set-3.dat", sep=" ")
bar3val <- c(bar3$V1, bar3$V2)
p1 <- hist(subset(bar2val, bar2val < 30), breaks=30)
p2 <- hist(subset(bar3val, bar3val < 30), breaks=30)
plot(p1, col=rgb(1,0,0,8/9), main="Barrier distribution", xlab="Barrier [kcal/mol]", ylab="Mutant count")
plot(p2, col=rgb(0,0,1,8/9), add=T)

Cualquier consejo sería muy apreciado.

Respuestas a la pregunta(2)

Su respuesta a la pregunta