Hi all,

I'm trying to retrieve the duplicate tags (DT:Z:SQ and DT:Z:LB) picard adds with RSamtools via the below in R, but I'm not successful, someone might know how to change this to retrieve this info?

DT = ScanBamParam(tag=c("DT"), what="flag")
bam.DT = scanBam(bamFile, param=DT)
str(bam.DT[[1]][["tag"]])
List of 1

$ DT: NULL