Skip to content
This repository has been archived by the owner on Mar 19, 2019. It is now read-only.

query: bammarkduplicates2 - are QCfail included #60

Closed
keiranmraine opened this issue Jan 22, 2018 · 2 comments
Closed

query: bammarkduplicates2 - are QCfail included #60

keiranmraine opened this issue Jan 22, 2018 · 2 comments

Comments

@keiranmraine
Copy link

If a read-pair has one end marked as qc-fail (512 / 0x200) how would bammarkduplicates2 interpret this if another pair has the same r1/r2 mapping data but doesn't have 0x200 set?

@gt1
Copy link
Owner

gt1 commented Jan 23, 2018

bammarkduplicates2 excludes QCfail entries from processing, so if one end is marked and the other is not that will effectively turn the pair into a single end read. If there is a pair mapping to the same coordinates, then a single end at the same position should be marked as a duplicate.

@keiranmraine
Copy link
Author

thanks

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants