bam_file <- system.file("extdata", "accepted_hits.bam", package="rbamtools")
reader<-bamReader(bam_file)
isOpen(reader)
header<-getHeader(reader)
refSeqs<-getRefSeqDict(header)
dim(refSeqs)
ht<-getHeaderText(refSeqs)
df<-as.data.frame(refSeqs)
refSeqs[1,1]
refSeqs[1,1]<-"newSeqName"
addSeq(refSeqs,SN="nextSeqName",LN=1000)
bamClose(reader)
Run the code above in your browser using DataLab