plot improved error rate confidence interval

[tdelhomme]

Plot a better estimated confidence interval (0.01, 0.99) according to the data.
At this time, we compute only the confidence interval for the maximum coverage with qnbinom(), and plot a line from 0 to this value (see big blue dotted lines on the following plot).
The issue is that the shape of this interval is not necessary linear.
Now we would like to compute this confidence interval for a sufficient number of points to be sufficiently accurate. Let's take 1000 equidistant points in {0,max(coverage)} if max(coverage)>1000, otherwise all coverage between 0 and the maximum coverage.

Use this opportunity to change colors for error rate and confidence interval from blue to black to erase conflict with the contour for qvalue=10.

[tdelhomme]

why minor bug?

[tdelhomme]

I guessed "bug" was only for pipeline crashing...

[mfoll]

The current line is wrong, so it's a bug. It's still minor as it's only visual.