Figure 1

Flowchart of MUSIC. H3K36me3 ChIP and control data in the region chr1:55,170,679-55,240,996 of the K562 cell line is used for illustrating the signal processing steps. ChIP and control reads (represented by short horizontal lines) are filtered for duplicates (red colored) and control signal is normalized with respect to ChIP signal (step 1). RD profiles are generated (step 2). The ChIP-Seq profile is corrected for mappability (labeled 'Mappability Corrected Signal') using the multi-mappability profile. Note the region indicated between the dashed lines that has low signal because of low mappability is filled with correction (step 3). Seven scale decomposition of the ChIP-Seq signal is computed. Under each decomposition, the ERs with the corresponding local minima are shown (step 4). The connected window shows the processing performed for generating the SSERs at each scale. The mappability-corrected signal is smoothed, the local minima are identified and the candidate ERs are formed (shown in grey), then candidate ERs are trimmed and filtered (shown in red) with respect to significance to identify the SSERS (shown in green.) SSERs for each scale are shown under the corresponding decomposition (step 5). The final set of ERs is formed by merging the SSERs (step 6).