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'

heatmap-3.png

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
../_images/cli_heatmap-3_4_1.png