import gosling as gos
data = gos.multivec(
url="https://resgen.io/api/v1/tileset_info/?d=UvVPeLHuRDiYA3qwFlm7xQ",
row="sample",
column="position",
value="peak",
categories=["sample 1"],
)
base = gos.Track(data).encode(y="peak:Q").properties(width=800, height=200)
top = gos.overlay(
base.mark_line().encode(
x=gos.X("position:G", domain=gos.GenomicDomain(chromosome="chr1"), axis="top")
),
base.mark_brush().encode(
x=gos.X(
"position:G",
domain=gos.GenomicDomain(chromosome="chr1"),
axis="top",
linkingId="linking-with-brush"
),
color=gos.value("steelBlue")
),
)
bottom = base.mark_line(background="steelBlue", backgroundOpacity=0.1).encode(
x=gos.X(
"position:G",
domain=gos.GenomicDomain(chromosome="chr1", interval=[200000000, 220000000]),
axis="top",
linkingId="linking-with-brush",
)
)
gos.stack(top, bottom)