expected standard distribution of l in block 310 and has mean u03bc and standard deviation u03c3 the normalization is performed using coverage c on a by filtering out the vectors v that are unlikely to represent the coverage c using an upper and lower cut off given by the user the library is recomputed using the output of the last step positions i not used to obtain d u2032 are reserved in block 310 for each position i in the assembly a a subset of reads s i u25a1 l that overlap the position i is collected in vector v i the mean u03bc i and standard deviation u221a c i u00b7 u03c3 i are calculated from s i in block 312 in block 314 of fig3 b the deviation of u03bc i from u03bc of the library is computed in block 316 the deviation of u221a c i u00b7 u03c3 i from u03c3 of the library is determined thresholds are used to determine unusual deviations i e deviations outside the thresholds in u03bc i and u221a c i u00b7 u03c3 i in block 318 the results may be output to a display device for user analysis in block 320 for each position i in the assembly when mean u03bc i deviates from the expected by more than a given threshold or standard deviation u221a c i u00b7 u03c3 i is above a given threshold the position i is flagged as.