bed#

[1]:
import cooler

import trackc as tc
[2]:
region = "chr12:31000000-34000000"
clr = cooler.Cooler("../../trackc_data/tutorials/bed/chr12_20_35Mb.cool")
tadfile = "../../trackc_data/tutorials/bed/chr12_20-35Mb.tad.bed3"
isfile = "../../trackc_data/tutorials/bed/chr12_20-35Mb.insulation.bed5"

mat = tc.tl.extractCisContact(clr=clr, region=region)

bed styles#

[3]:
ten = tc.tenon(figsize=(6, 1))
ten.add(pos="bottom", height=1.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)
ten.add(pos="bottom", height=0.8, hspace=0.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)
ten.add(pos="bottom", height=0.5, hspace=0.1)

tc.pl.mapC(ax=ten.axs(0), mat=mat, label="cmap", map_type="tri", height=40)

tc.pl.bed_track(
    ax=ten.axs(1),
    bed=tadfile,
    regions=region,
    style="rec",
    intervals=2,
    label="TAD\nstyle:rec",
)
tc.pl.bed_track(
    ax=ten.axs(2),
    bed=isfile,
    regions=region,
    style="line",
    label="insulation score\nstyle:line",
)

tc.pl.bed_track(
    ax=ten.axs(3), bed=tadfile, regions=region, style="tri", label="style:tri"
)
tc.pl.bed_track(
    ax=ten.axs(4),
    bed=tadfile,
    regions=region,
    style="link",
    invert_y=True,
    label="link:invert_y",
)

tc.pl.bed_track(
    ax=ten.axs(5), bed=tadfile, regions=region, style="link", label="style:link"
)
tc.pl.bed_track(
    ax=ten.axs(6),
    bed=tadfile,
    regions=region,
    style="tri",
    invert_y=True,
    label="tri:invert_y",
)

tc.pl.bed_track(
    ax=ten.axs(7), bed=isfile, regions=region, style="bar", label="style:bar"
)
tc.pl.bed_track(
    ax=ten.axs(8),
    bed=isfile,
    regions=region,
    style="rec",
    cmap="Blues",
    invert_y=True,
    label="invert_y",
)
tc.pl.bed_track(
    ax=ten.axs(9),
    bed=isfile,
    regions=region,
    vmin=0.015,
    vmax=0.025,
    tick_fl="%0.3f",
    style="bar",
    invert_y=True,
    label="style:bar",
)

tc.pl.scale_track(
    ax=ten.axs(9),
    region=region,
    scale_adjust="Mb",
    tick_pos="bottom",
    ratio2ax=1.5,
    tick_rotation=45,
)
# tc.savefig('bed_track.pdf')
no max min range
maxrange: 0.040769440000000004 minrange: 0.0003320906
../_images/track_types_bed_4_1.png