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.

How are supplementary alignments handled by HaplotypeCaller?

sp580sp580 GermanyMember

Hello,

I am curious about how HaplotypeCallers deals with secondary alignments. I know (based on the documentation and my own runs) that the filter NotSecondaryAlignmentReadFilter is applied but I do not see the filter NotSupplementaryAlignmentReadFilter (is it obsolete?, redundant?, hast it been renamed or merged with another filter maybe?)

Thanks in advance!

Best Answer

Answers

  • SChaluvadiSChaluvadi Member, Broadie, Moderator admin

    @sp580
    I believe that HaplotypeCaller filters out secondary alignments and only uniquely mapped reads are used for variant discovery. The NotSupplementaryAlignmentReadFilter is still present in GATK4 documentation with the same name. Please let us know if this is the information you were looking for!

Sign In or Register to comment.