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.

Sign up for GitHub

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

Jump to bottom

Are duplicated records in BAM file counted or not? #118

Open
zztin opened this issue Feb 17, 2023 · 1 comment
Open

Are duplicated records in BAM file counted or not? #118

zztin opened this issue Feb 17, 2023 · 1 comment

Comments

@zztin
Copy link

zztin commented Feb 17, 2023 •
edited
Loading

Hi,
I am wondering if ichorCNA (actually readCounter in HMMcopy) count the records in bam files if they contains flag indicating they are (PCR or optical) duplicated record (flag 0x400) by picard (MarkDuplicates)?
In theory I suppose it should not count these molecules, however, I could not find clear documentation if this check is implemented or not. Could you point me to the correct forum, or perhaps you have a clue to this?
Thank you very much!
@ilykos
Copy link

ilykos commented Dec 8, 2023

@zztin I took the liberty of replicating your question in the readCounter repository. HMMCopyUtils Issue 15. As ichorCNA works on downstream .wig files, I do not believe it has any awareness of the 0x400 flag.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@ilykos @zztin