Hi Gviz developers, thanks for creating this powerful genome visualisation package.
I used Gviz to visualise indels, however, I found that the annotation is not consistent with IGV (see figure); Gviz extended one more base to the right.
I wonder whether this is a bug or not.
Here is my code
read_track <- AlignmentsTrack(range = bam, stacking = "full", showIndels = TRUE, isPaired=TRUE)seq_track <- SequenceTrack(genome)plotTracks(list(seq_track, read_track), chromosome = chr, from = start, to = end)
Hi Gviz developers, thanks for creating this powerful genome visualisation package. I used Gviz to visualise indels, however, I found that the annotation is not consistent with IGV (see figure); Gviz extended one more base to the right.
I wonder whether this is a bug or not. Here is my code
read_track <- AlignmentsTrack(range = bam, stacking = "full", showIndels = TRUE, isPaired=TRUE)
seq_track <- SequenceTrack(genome)
plotTracks(list(seq_track, read_track), chromosome = chr, from = start, to = end)
Thanks!