Entering edit mode
9.7 years ago
tonja.r
▴
600
I am using gviz from a bioconductor to plot SNPs data using chromosome and genome axis:
My problem is that it plots dense located SNPs under each other and I could not find a parameter to avoid. I also could not find if it is possible to plot the name of the SNPs to the corresponding bar or mark one specific SNP with another color.
gr99 <- GRanges(seqnames = c("chr7","chr7","chr7","chr7","chr7","chr7"),
IRanges(start=c(128580042,128580052,128575552,128576086,128579202,128579666),
end=c(128580042,128580052,128575552,128576086,128579202,128579666)))
atrack2=AnnotationTrack(gr99, name = "SNP")
gen = c()
gen[chr] =genome(atrack2)="hg19"
gtrack <- GenomeAxisTrack()
itrack <- IdeogramTrack(genome = gen, chromosome = "chr7")
plotTracks(list(itrack,gtrack, atrack2),panel.only == TRUE)
Not the best solution. Make two or more tracks and plot them together