def __init__(self, tcc, config): #get the absolute profile for this stretch self.tcc = tcc #this uses the merged wig files... self.profile = svs.svCoord([tcc], config)
for block in stretch.blocks: print block ''' ''' #histogram test a = [5,5,5,5,5,5,4,4,4,3,3,2,1] cgPlot.plotHistogram(a, name = 'testHistogram') ''' #svCoord test tcc = 'chr1:1:1103298:1103298' eLevels = sv.svCoord([tcc], 'siPeaks.conf') for coord in eLevels: print coord, eLevels[coord] ''' tcc = 'chr1:1:1102539:1103319' eLevels = sv.profileAroundPoint(tcc, 50, 'agoProfile.conf') sorted = eLevels.keys() sorted.sort() for coord in sorted: print coord, eLevels[coord] ''' #Plotting '''