counts = data.frame(hic=dim(hic)[1]/303984633*4*10^8,hic_actd=dim(hic_actd)[1]/171680362*4*10^8,
                    nhic=dim(nhic)[1]/332710486*4*10^8,nhic_actd=dim(nhic_actd)[1]/333186514*4*10^8)
barplot(as.matrix(counts), main="",
        xlab="", col=c("#2BA809","#F97164"),
         beside=T,ylim=c(0,3000000),ylab = 'Total interactions with rDNA')