Script from Steve to reanalyse old wgbs aligned data with updates script - add # cytosine site info