I was able to figure out a solution. There are a few factors at play. First, when using set style data histogram it only expects one column for the plotting command. Second, the natural x tic values start at 0 and increment by 1.
So we simply have to simply add each column 1 by 1.
For the x labels, if using gnuplot script, one can say plot for [COL=2:4:2] 'file.dat' using COL:xticlabels(1) where the x-axis label is then calculated from the first column of the data file. But given that we are in Julia code, I am not sure how to write the equivalent. Instead, a solution is to replace the xtic values manually.