circlize icon indicating copy to clipboard operation
circlize copied to clipboard

Can we plot SNPs on a circos plot ?

Open Krithika-Bhuvan opened this issue 2 years ago • 3 comments

Is it possible to plot single nulcletotide polymorphisms (SNPs) on a circos plot ? A SNP has the same start and end position. For example, 1:207620289 (rs12081383) 1:207624395 (rs72644189)

Krithika-Bhuvan avatar Oct 14 '22 19:10 Krithika-Bhuvan

Then you can draw them as points.

jokergoo avatar Oct 15 '22 15:10 jokergoo

Hello, I tried that but it doesn't work.

I put the same location in the "start" and "end" points. I also tried with "start" and "start + 1 " . Both times the track was plotted but no points. It seems like it will only work with a region where "start" and "end" points are far apart. Is there any workaround for this that you can suggest ?

Thanks, Krithika

On Sat, Oct 15, 2022 at 11:34 AM Zuguang Gu @.***> wrote:

Then you can draw them as points.

— Reply to this email directly, view it on GitHub https://github.com/jokergoo/circlize/issues/342#issuecomment-1279768110, or unsubscribe https://github.com/notifications/unsubscribe-auth/AANXTHC4A4Z6RQKESB3OYJ3WDLFHDANCNFSM6AAAAAARFQWN7E . You are receiving this because you authored the thread.Message ID: @.***>

Krithika-Bhuvan avatar Oct 15 '22 23:10 Krithika-Bhuvan

No, it works for any type of genomic intervals:

circos.initializeWithIdeogram(plotType = "axis")
bed = data.frame(chr = "chr1", start = 1e8, end = 1e8, value = 1)
circos.genomicTrack(bed, ylim = c(0, 2), panel.fun = function(gr, value, ...) circos.genomicPoints(gr, value))
image

Maybe you can send me some example dataset?

jokergoo avatar Oct 27 '22 10:10 jokergoo