Entering edit mode
Vinicius Henrique da Silva
▴
40
@vinicius-henrique-da-silva-6713
Last seen 17 months ago
Brazil
Is it possible to change the fill color of the reads in a track obtained from the AlignmentsTrack()
function, library(Gviz)
package based on its mapping quality (i.e. mapq)?
I searched for something similar to groupAnnotation from AnnotationTrack()
but it seems that such a parameter is not present in for AlignmentsTrack()
.
I tried the below code but it does not work:
afrom <- 44945200
ato <- 44947200
alTrack <- AlignmentsTrack(
system.file(package = "Gviz", "extdata", "snps.bam"), isPaired = TRUE,
groupAnnotation="mapq")
plotTracks(alTrack, chromosome = "chr21", from = afrom, to = ato)
Any idea to get around this?