Hi,

I have a set of samples and log10 total read counts of them in a data frame like below

    > head(a)
        log10_total_counts
    A1            6.468503
    A10           6.565213
    A11           6.752139
    A12           5.078598
    A2            6.277342
    A3            6.473411
    >

![enter image description here][1]


  [1]: https://i.ibb.co/P6mmKc5/Picture1.png


For instance red dots have less than 1000000 counts and black dots have more thsn 1000000 counts. 

How could I repreduce this plot please?
Something like the following would work with base R functions:


    # make some data
    a <- data.frame(sample=paste0("A", 1:50), log10_total_counts=rnorm(50, 6))
    
    # get x points
    x <- 1:length(a$log10_total_counts)
    
    # create plot, drop x axis
    plot(x, a$log10_total_counts, ylim=c(0,8), pch=19, xaxt="n", xlab="sample")
    # add dashed line and X-axis labels with rotation
    abline(h=6, lty="dashed")
    axis(1,1:50, labels=a$sample, las=3)
    
    # select points to color differently
    LessThan6 <- a$log10_total_counts < 6
    points(x[LessThan6], a$log10_total_counts[LessThan6], pch=19, col="red")

