Question: Rsamtools: error with quickBamFlagSummary
0
gravatar for bsmith030465
16 months ago by
bsmith030465140
United States
bsmith030465140 wrote:

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!

rna-seq rsamtools • 326 views
ADD COMMENTlink written 16 months ago by bsmith030465140
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1842 users visited in the last hour