Holiday Notice:
The Frontline Support team will be slow to respond December 17-18 due to an institute-wide retreat and offline December 22- January 1, while the institute is closed. Thank you for your patience during these next few weeks. Happy Holidays!

How does gatk deal with bam XA tag when calling variations?

Hello,

When a read has mutli-hits on a reference genome, the bam file would record mutli positions in XA tag. So, would GATK consider the XA tag when calling variations?
After getting filtered VCF using gatk, I found that some sites were covered by many reads (for example, more than ten). But when I check the bam file in the IGV viewer, these sites didn't have any reads but the consensus sequence was intact. So, I wonder whether GATK also count the reads in the XA tag when calling variations?

Thanks a lot!

Tagged:

Issue · Github
by Sheila

Issue Number
2471
State
closed
Last Updated
Assignee
Array
Milestone
Array
Closed By
sooheelee

Best Answer

Answers

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    @liujing251614
    Hi,

    I have asked another team mate to respond, as she knows more. But, you may find this thread useful in the meantime.

    -Sheila

  • liujing251614liujing251614 ChinaMember

    Yeah, thank you~

  • liujing251614liujing251614 ChinaMember

    I get it, thanks a lot!

Sign In or Register to comment.