How can I selectively colour the entire area of a plot based on specified points on the x-axis?

Question

I would like to colour the background of a line plot based on manually inserted coordinates on the x-axis.

As a simple example:

x <- 1:10
y <- c(