CIGAR String error

[ashokrags]

Hello, I am trying to run bayesembler on a bam generated by GSNAP and have run into a CIGAR string error as below:

bayesembler -s first -b ../GsnapAlignments/Gsnap.T1990_ctx.33080.cmu059.19.Aug.2014/T1990_ctx.all.mm.rg.dup.srtd.bam

You are using the Bayesembler v1.2.0. For more information go to bayesembler.binf.ku.dk

bam_nd_pe_plus_file_nameT1990_ctxallmmrgdupsrtd_nd_plus.bam [24/06/2015 16:27:56] Removing duplicate reads ERROR: Unhandled cigar string symbol 'S'!

I will be much obliged for any insight into this issue Cheers Ashok

[lassemaretty]

Hi Ashok,

Thank you for posting. The Byesembler currently only supports alignments generated by TopHat, however the next release (expected to appear soon) will also provide support for the STAR and HISAT aligners. If you need it, we will try to also include support for GSNAP.

Best regards,

Lasse

[ashokrags]

hi Lasse, Thanks so much for your quick reply. It will be great to include support for GSNAP. However, i just wanted to check with you as to why soft clipped reads are not recognized by the program. I can probably just write a wrapper to modify the bam so that soft-clipping is converted to hard clips if you think that will work.

[ashokrags]

Hi Lasse, FYI. I removed soft-clipped reads from the bam file I was working on. Now it appears the even hard clipped reads pose a problem. Just wanted to let you know [ar474@cn022 bayesembler]$ bayesembler -b T1990_ctx.all.mm.rg.dup.srtd.nosoft.bam -s first -f 4

You are using the Bayesembler v1.2.0. For more information go to bayesembler.binf.ku.dk

bam_nd_pe_plus_file_nameT1990_ctxallmmrgdupsrtdnosoft_nd_plus.bam [01/07/2015 13:25:09] Removing duplicate reads ERROR: Unhandled cigar string symbol 'H'!

[lassemaretty]

Hi,

Sorry for not getting back to you before now. Neither soft nor hard clips are currently supported. My best advice is to wait for the next release or use TopHat2 for mapping. I'm sorry for the inconvenience.

Best,

Lasse

[ashokrags]

hi Lasse, Not an issue. I already had a script to remove soft clip from the cigar and correct the read and quality in the bam accordingly. I didn't realize the hard clips were an issue either. I am currently removing hard clipped reads and will try to run bayesembler

[lassemaretty]

sounds great, let us know if you run into problems again.

[ashokrags]

hi Lasse, sorry to be a bother, but does the program not support multiply mapped reads? I know I don't have an issue with Cufflinks when I use the same bam, after fixing the soft clip issue. Also is there an option to not remove duplicate reads?