library(rbamtools)
bam_file <- system.file("extdata", "accepted_hits.bam", package="rbamtools")
reader<-bamReader(bam_file)
idx_file<-"index.bam.bai"
createIndex(reader,idx_file)
loadIndex(reader,idx_file)
index.initialized(reader)
bamClose(reader)
Run the code above in your browser using DataLab