I would like to plot aligned reads from a
bam file over a specific region. Everything works fine with the following command:
alTrack <- AlignmentsTrack(my.bam.file, isPaired=TRUE) snp.pos <- 308632 plotTracks(alTrack, chromosome="chr1", from=snp.pos-200, to=snp.pos+100)
But shortening the region (by increasing
from) results in an error:
plotTracks(alTrack, chromosome="chr1", from=snp.pos-100, to=snp.pos+100) Error in data.frame(x1 = start(pairGaps) - 1, y1 = gy, x2 = end(pairGaps) + : arguments imply differing number of rows: 0, 1
How can I fix this?