I get the following error from snp_store (on output from splazers):
WARNING! Read TESTING: cigar alignment length does not match genome coordinates. Discarding read..
I'm fairly certain this error stems from splazers as the alignment length computed from the cigar string does not match the start and end match position (off by 1 bp). Interestingly, when I truncate 300 bps same sequence from 478 bps to 178 bps, the cigar string appears to be correct. Note that 32M is now 31M which is the correct length.
cigar of the longer sequence (478 bps):
cigar of the truncated sequence (178 bps):
Has anyone else encountered this issue and is there a fix?
Tina Hu, Ph.D.
1 Little West 12th St.
New York, NY 10014
This email message and any attachments being sent by Celmatix, Inc., are confidential, and may be privileged. If you are not the intended recipient, please notify us immediately—by replying to this message—and destroy all copies of this message and any attachments. Thank you.