Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
Principle of removing duplicated reads in Picard
MarkDuplicates of Picard is a useful function to remove duplicated reads. However, after reading the introduction of Picard (https://broadinstitute.github.io/picard/command-line-overview.html#MarkDuplicates), I still have some questions about filtering out the Duplicated reads. Here are the questions:
The source of PCR duplicated includes "library/PCR-generated duplicates (LB)" and "sequencing-platform artifact duplicates (SQ)". How does Picard identify LB and SQ from reads?
In default setting with REMOVE_DUPLICATES=true, which type of duplicated reads will be removed, SQ, LB, or both?
The reads A, B, C are considered as duplicated reads, and their quality scores are equal. If these reads are mapped to the same position in genome, which reads will be removed after filtering by Picard?
And if these reads are mapped to the different position in genome, which reads will be removed after filtering by Picard?
Continue to the previous question, but the qualities of read A, B, C are not equal, which reads will be removed after filtering by Picard?