bigwig

Get test data

data url: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM4604189

!curl https://ftp.ncbi.nlm.nih.gov/geo/samples/GSM4604nnn/GSM4604189/suppl/GSM4604189_1360_CUT_TAG_H3K27ac.MACS2.nodup_x_ctl.pval.signal.bigwig -O GSM4604189_H3K27ac.bw

import pyBigWig

import trackc as tc

AML_1360_H3K27ac = pyBigWig.open(
    "../../trackc_data/tutorials/bigwig/GSM4604189_H3K27ac.bw"
)

ten = tc.tenon(figsize=(8, 1))
ten.add(pos="bottom", height=1, hspace=0.1)
ten.add(pos="bottom", height=1, hspace=0.2)
ten.add(pos="bottom", height=0.5, hspace=0.2)


regions = ["chr8:127000000-129200000", "chr14:96500000-99300000"]
tc.pl.bw_track(
    AML_1360_H3K27ac,
    ten.axs(0),
    regions=regions,
    vmax=20,
    label="H3K27ac",
    binsize=10000,
    primary_col="tab:blue",
)
tc.pl.bw_track(
    AML_1360_H3K27ac,
    ten.axs(1),
    regions=regions,
    vmax=5,
    label="H3K27ac",
    binsize=10000,
    invert_y=True,
    primary_col="tab:purple",
)
tc.pl.multi_scale_track(
    ten.axs(2), regions=regions, scale_adjust="Mb", intervals=1, tick_fontsize=10
)

# tc..savefig('trackc_bigwig.pdf')