Box plot for rna seq data

Hi friends
I plotted this box-wisker for TCGA HTSeq data in R.
I want to have harf of them as red and half as blue (control vs treatment groups). or is there any better way for boxplot?
How can I do that?
I just used simple code:

data <- read.table("mydata.csv", header=TRUE, row.names = 1)
library(DESeq2)
mat <- as.matrix(data)
log <- rlog(mat)
boxplot(log)

my plot:
ibb.co/xXsHtC0