mapc_markline

The trackc.mapc_markline function used for mark contact regions on the heatmap.

It should be used in conjunction with trackc.mapc.

import trackc as tc

BxPC3_chr18_25k = "../../trackc_data/examples/BxPC3.chr18.mcool::/resolutions/25000"

regions = ["18:47400000-48280000", "18:75280000-74030000"]
neo_domain_regions = ["18:47950000-48280000", "18:75280000-74850000"]
tumor_zoom = tc.tl.extractContactRegions(clr=BxPC3_chr18_25k, row_regions=regions)

ten = tc.tenon(figsize=(6, 1))
ten.add(pos="bottom", height=1.4, hspace=0.05)
ten.add(pos="bottom", height=1.4, hspace=0.05)

ten.add(pos="bottom", height=1, hspace=0.05)


tc.pl.mapC(
    ax=ten.axs(0),
    mat=tumor_zoom.cmat,
    map_type="triangle",
    maxrange=200,
    minrange=10,
    ax_on=False,
    height=40,
)
tc.pl.mapc_markline(
    ax=ten.axs(0),
    row_regions=tumor_zoom.row_regions,
    map_type="triangle",
    mark_regions=neo_domain_regions,
    binsize=25000,
    map_order=0,
    symmetric=False,
    show_regions_edge=False,
)


tc.pl.mapC(
    ax=ten.axs(1),
    mat2=tumor_zoom.cmat,
    map_type="triangle",
    cmap=tc.pa.washu,
    maxrange=200,
    minrange=10,
    ax_on=False,
    height=40,
)
tc.pl.mapc_markline(
    ax=ten.axs(1),
    row_regions=tumor_zoom.row_regions,
    map_type="triangle",
    mark_regions=neo_domain_regions,
    binsize=25000,
    map_order=1,
    symmetric=False,
    show_regions_edge=True,
)

tc.pl.multi_scale_track(
    ax=ten.axs(2), regions=regions, scale_adjust="Mb", intervals=2, tick_rotation=0
)

# tc.savefig('trackc_markline.pdf')

maxrange: 200 minrange: 10
maxrange: 200 minrange: 10