DNA double-strand breaks (DSBs), are a major threat to genomic stability and may lead to cancer. Several technologies to accurately detect DSBs genome-wide have been developed recently, but still lacking publicly available tools for analysis of the resulting data. Here, we present a step-by-step iSeq package (http://breakome.utmb.edu/software.html), custom designed for analysis and interpretation of DSB-sequencing data. iSeq performs barcode trimming and read counting, and identifies DSB-enriched regions by statistical test and annotate them to the desired genomic features. Applying this package, users can identify and annotate DSB-enriched regions from base pair (eg. Cas9 cleavage sites) up to megabase (eg. DNA replication stress-induced) resolution, and if possible quantify DSB frequencies per cell genome-wide by combining with qDSB-Seq. iSeq can be used for any sequencing-based DSB detection techniques. The analysis for Steps 1-19 can be performed within ~4 hours.Recently, we presented a general method, qDSB-Seq 13 , to quantify DSBs genome-wide, which provides both DSB frequencies per cell and their precise genomic coordinates. qDSB-Seq has been applied to BLESS 4 and i-BLESS 5 , it is also available for END-Seq 6 , Break-Seq 7 , DSBCapture 8 , BLISS 15 , and so on. In this work, we induced spike-in DSBs by a site-specific endonuclease and used these induced DSBs to quantify studied DSBs. qDSB-Seq method has