Script for analysis of SEAseq data