Added PASSeq support to call_peaks. Also cleaned up the gene regions discovery code.