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

How to only mark duplicates without removing them? #11

Open
Dries-B opened this issue Apr 21, 2020 · 1 comment
Open

How to only mark duplicates without removing them? #11

Dries-B opened this issue Apr 21, 2020 · 1 comment

Comments

@Dries-B
Copy link

Dries-B commented Apr 21, 2020

For my workflow, I want to combine all information, including duplicates, into one BAM file. I am looking for similar behaviour to Picard's MarkDuplicates.

I was hoping --annotated_output_bam would give the desired output, but it does not contain the duplicate flag 1024. This file does contain tags though; could these be used for subsequent filtering?

@cgates
Copy link
Member

cgates commented Apr 21, 2020

Hi @Dries-B

Thanks for you interest in Connor. Takara Bio is the primary support for Connor and I've forwarded your question to them.

For my own part, I believe I understand your question and while Connor does not set the duplicate bit in the SAM flag, the annotated BAM file contains an undocumented tag, X6, which indicates that specific alignment would be the template for the consensus alignment sequence. So the absence of the X6 tag is roughly analogous to the duplicate SAM flag bit. I hope that helps.

Best,
cg

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

2 participants