We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

SplitNCigarReads Unsupported Cigar Operator: =

Hi,

I got this error when trying to run SplitNCigarReads on a BAM file. I checked the source that generates this exception, and it appears that (as may be expected from the name), there is no handling for some Cigar String Operators.

The BAM specification states that '=' is perfectly fine. So... what gives?

Issue · Github
by Sheila

Issue Number
1191
State
closed
Last Updated
Assignee
Array
Milestone
Array
Closed By
vdauwera

Best Answer

Answers

  • SheilaSheila Broad InstituteMember, Broadie ✭✭✭✭✭

    @cacampbell
    Hi,

    Can you please tell us the exact command you ran and the version you are using?

    -Sheila

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    @cacampbell To specifically address your original question: you're correct that this tool does not handle the '=' operator, even though the spec allows for it. This is a limitation of the current incarnation of the tool due to an oversight at the time it was written -- we rarely encounter those operators in the data we see so we did not think to test for them . I believe the updated version that is going into GATK4 is not subject to that limitation, so feel free to give that a spin when it goes in if you'd like.

Sign In or Register to comment.