def addPeaks(bedDir,timepoint,direction,expression):
records=BedReader.readAll(bedDir+"/"+timepoint+"_"+direction+".bed")
for rec in records:
array=expression.get(rec.chr,None)
if(array is None): array=expression[rec.chr]=[]
interval=rec.interval
interval.begin=interval.intCenter()
interval.end=interval.begin+1
interval.type="peak"
interval.dir=direction
array.append(interval)
