), Track('H3K9me3',1e5,g.sizes).from_bigWig( '/mnt/d/data/sniper/histone_marks/bigWigs_HCT116/Rao-2017-CHIP019-H3K9me3-untreated.bw' ), Track('H3K79me2',1e5,g.sizes).from_bigWig( '/mnt/d/data/sniper/histone_marks/bigWigs_HCT116/Rao-2017-CHIP023-H3K79me2-untreated.bw' ), Track('H4K20me1',1e5,g.sizes).from_bigWig( '/mnt/d/data/sniper/histone_marks/bigWigs_HCT116/Rao-2017-CHIP025-H4K20me3-untreated.bw' ), )) )), Sequential('subcompartments', tracks= ( Track('cohesin_depleted',1e5,g.sizes).from_bed('data/cd_cd_HCT116.bed'), Track('wild_type',1e5,g.sizes).from_bed('data/cd_wt_HCT116.bed'), )) ) """""" for track in tracks: g.add_track(track) g.save('genomes/HCT116.mg') pprint.pprint(g.get('chr2',1e5,5e5,to_numpy=True,tracks=[ 'histone_marks/cohesin_depleted/H3K27ac', 'histone_marks/wild_type/H3K27ac', 'histone_marks/cohesin_depleted/H3K36me3', 'histone_marks/wild_type/H3K36me3', ]))