-
Notifications
You must be signed in to change notification settings - Fork 12
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Pileup reads with adjacent indels. #28
Comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
CIGAR code with adjacent indels (e.g. 1M2D3I, 1M2N3I) is allowed by SAM file format specification.
Insertions and deletions are piled up at the base position just before them according to samtools.
I couldn't figure out the specification for mpileup format, but mpileup output of samtools 1.3.1 seems to be incorrect.
input SAM
samtools 1.3.1 output
cljam output
I think inserted sequence
+5GCATC
should be+5TGCAT
.The last base
C
of the samtools output is affected by upper nibble of the first quality score.There might be an off-by-one error & buffer overrun or something in htslib/samtools implementation.
What should cljam output? Should it be the same as samtools?
The text was updated successfully, but these errors were encountered: