Entering edit mode
2.2 years ago
noodle
▴
590
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