8000 Fix mate sequence when reading sam/bam by jbedo · Pull Request #31 · dcjones/quip · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

Fix mate sequence when reading sam/bam #31

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

Open
wants to merge 1 commit into
base: master
Choose a base branch
from

Conversation

jbedo
Copy link
@jbedo jbedo commented May 31, 2022

Previously upon reading the case of tid == mtid was detected and the
sequence name mapped to "=". This causes missing sequence name errors
upon decompression. As the case of tid == mtid is handled during writing
of sam/bam, this patch simply records the full mate sequence name,
resolving the matching issues.

Example read after decompression pre patch:

SL1344_1_530_0:0:0_0:0:0_6c9    163     SL1344  1       60      70M     *       461     530     AGAGATTACGTCTGGTTGCAAGAGATCATGACAGGGGGAATTGGTTGAAAATAAATATATCGCCAGCAGC  IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII       MQ:i:60 AS:i:70 RG:Z:mysample1  NM:i:0  MC:Z:70M        MD:Z:70 ms:i:2800       XS:i:0

and post patch:

SL1344_1_530_0:0:0_0:0:0_6c9    163     SL1344  1       60      70M     =       461     530     AGAGATTACGTCTGGTTGCAAGAGATCATGACAGGGGGAATTGGTTGAAAATAAATATATCGCCAGCAGC  IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII       MQ:i:60 AS:i:70 RG:Z:mysample1  NM:i:0  MC:Z:70M        MD:Z:70 ms:i:2800       XS:i:0

Previously upon reading the case of tid == mtid was detected and the
sequence name mapped to "=". This causes missing sequence name errors
upon decompression. As the case of tid == mtid is handled during writing
of sam/bam, this patch simply records the full mate sequence name,
resolving the matching issues.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant
0