Rsamtools: error with quickBamFlagSummary

Question

[bsmith030465]

Hi,

I was trying to determine if my bam file was single or pair ended. However, when I try to use Rsamtools, I get an error:

> quickBamFlagSummary(file="/Users/bsmith/Documents/run_s1_Aligned_Reads.bam")
Error: identical(length(N_1seg_rec_per_uqname) + length(N_mseg_rec_per_uqname),  .... is not TRUE

> traceback()

9: stop(msg, call. = FALSE, domain = NA)
8: stopifnot(identical(length(N_1seg_rec_per_uqname) + length(N_mseg_rec_per_uqname), 
   length(N_rec_per_uqname)))
7: quickBamFlagSummary(res0, param = param, main.groups.only = main.groups.only)
6: quickBamFlagSummary(res0, param = param, main.groups.only = main.groups.only)
5: quickBamFlagSummary(file, param = param, main.groups.only = main.groups.only)
4: quickBamFlagSummary(file, param = param, main.groups.only = main.groups.only)
3: .local(file, ..., param = param, main.groups.only = main.groups.only)
2: quickBamFlagSummary(file = "/Users/bsmith/Documents/run_s1_Aligned_Reads.bam")
1: quickBamFlagSummary(file = "/Users/bsmith/Documents/run_s1_Aligned_Reads.bam")

Is there something I need to do differently? Its a RNA seq data, if that helps any...

thanks!