script to use bsseq with output