Skip to content
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

gencore seems to have wrongly modified my SAM alignment record #47

Open
jammywwh opened this issue Feb 21, 2023 · 0 comments
Open

gencore seems to have wrongly modified my SAM alignment record #47

jammywwh opened this issue Feb 21, 2023 · 0 comments

Comments

@jammywwh
Copy link

jammywwh commented Feb 21, 2023

Hi,
Gencore seems to have wrongly modified my SAM alignment record.
I have SAM record as below:
NS500830:1150:HWNN2AFX3:2:11202:24875:15467:GGT_GGC 161 chr7 55268915 60 71M chrUn_gl000220 132375 0 ACTCCAACTTCTACCGTGCCCTGATGGATGAAGAAGACATGGACGACGTGGTGGATGCCGACGAGTACCTC /AEEEAEEEEAEEEEEEEEEEEAEEEAAEEEEEEEEEEEEEEEEE/AEEEAEEEEEEEEEEAEEEE/EEE< NM:i:0 MD:Z:71 MC:Z:71M AS:i:71 XS:i:19 RG:Z:BT2302040049LNCBC NS500830:1150:HWNN2AFX3:2:11202:24875:15467:GGT_GGC 81 chrUn_gl000220 132375 57 71M chr7 55268915 0 TCTCTCTCTCTCTCTCTCTCTGCCTCTCTCACTGTGTCTGTCTTCTGTCTTATTCTCTTTCTCTCTCTGTC EEEEEEAEEEEEEEEEEEEEEEAEEEEEEEEAEAEE/EEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEE NM:i:0 MD:Z:71 MC:Z:71M AS:i:71 XS:i:52 RG:Z:BT2302040049LNCBC
but, after gencore remove duplicate reads using s=2, I got
NS500830:1150:HWNN2AFX3:2:11202:24875:15467:GGT_GGC 161 chr7 55268915 60 71M chrUn_gl000220 132375 0 ACTCCAACTTCTACCGTGCCCTGATGGATGAAGAAGACATGGACGACGTGGTGGATGCCGACGAGTACCTC EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE NM:i:0 MD:Z:71 MC:Z:71M AS:i:71 XS:i:19 RG:Z:BT2302040049LNCBC NS500830:1150:HWNN2AFX3:2:11202:24875:15467:GGT_GGC 145 chrUn_gl000220 132375 25 71M chr7 55268915 0 TCTCTCTCTCTCTCTCTCTCTGCCTCTCTCACTGTGTCTGTCTTCTGTCTTATTCTCTTTCTCTCTCTGTC EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE NM:i:1 MD:Z:7C63 MC:Z:69M2S AS:i:66 XS:i:57 RG:Z:BT2302040049LNCBC XA:Z:chrY,-10024049,71M,4;

why 81 convert to 145?
Thx

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

No branches or pull requests

1 participant