Create histogram of number of transposon per gene

%%%%histogramme #TN per #gene

figure(2)

hist(tnpergene, [0:1:250],'b')

h.FaceColor = [1 0 0];

h.EdgeColor = [1 0 0];

hold on

%%%%histogramme #TN per #essential

hist(tnperessential,[0:1:250])

h.FaceColor = [0 1 0];

h.EdgeColor = [0 1 0];

hold on

%%%%histogramme #TN per #gene - histogramme #TN per #geneessential = #TN

%%%%per non-essential

histnonessential=hist(tnpergene,[0:1:250])-hist(tnperessential,[0:1:250]);

bar([0:1:250],histnonessential,'g');