bam_file <- system.file("extdata", "accepted_hits.bam", package="rbamtools")
reader<-bamReader(bam_file)
isOpen(reader)
header<-getHeader(reader)
headLine<-getHeaderLine(header)
getVal(headLine,"VN")
getVal(headLine,"SO")
setVal(headLine,"SO","unsorted")
l<-as.list(headLine)
txt<-getHeaderText(headLine)
bamClose(reader)
Run the code above in your browser using DataLab