rearranged interaction#
CLI mode#
conf-hicmap_3.yml
trackc:
- ax: t1
height: 2
track_type: hicmap
track_para:
mat:
method: extractContactRegions
clr: /path/GSM4417639_SK-N-DZ_C6BC81F2_b38d5.mcool::/resolutions/5120000
mapC:
cmap: PuBu
map_type: tri
maxrange: 10000
- ax: t2
height: 0.6
track_type: multi_scale_track
track_para:
scale_adjust: Mb
intervals: 2
trackc cli conf-hicmap_3.yml -s 4,1 -r 'chr6:0-25600000 chr8:122880000-143360000'
API mode#
[1]:
import trackc as tc
[6]:
ten = tc.tenon(figsize=(4, 1))
ten.add(pos="bottom", height=2)
ten.add(pos="bottom", height=0.6)
regions = ["chr6:0-25600000", "chr8:122880000-143360000"]
clr = "/Users/yuanzan/Documents/RD/GSM4417639_SK-N-DZ_C6BC81F2_b38d5.mcool::/resolutions/5120000"
mat = tc.tl.extractContactRegions(clr=clr, row_regions=regions)
tc.pl.mapC(
ax=ten.axs(0),
mat=mat.cmat,
map_type="tri",
ax_on=False,
cmap="PuBu",
maxrange=10000,
)
tc.pl.multi_scale_track(ax=ten.axs(1), regions=regions, scale_adjust="Mb", intervals=2)
# tc.savefig('trans_interactions.pdf')
maxrange: 10000 minrange: 77.0