perl script for BS-Seq analysis