script to plot GC skew