frame = \
BigWig(bw4, alpha=0.6) + Color("#ff9c9c") + \
BigWigCoverage(bw1, alpha=0.6, color="#fcbe07") + \
Title("RNA / H3K27ac") + \
BigWig(bw4, alpha=0.6) + Color("#ff9c9c") + \
BigWigCoverage(bw3, alpha=0.6, color="#07fcbe") + \
Title("RNA / H3K4me3") + \
BigWig(bw4, alpha=0.6) + Color("#ff9c9c") + \
BigWigCoverage(bw2, alpha=0.6, color="#66ccff") + \
Title("RNA / H3K27me3")
frame *= MaxValue(20)
frame *= MinValue(0)
frame * TrackHeight(4)
frame = XAxis() + frame
# frame.plot("chr1:1-5000000")
frame.plot("chr9:4000000-6000000")