MuTect2 annotations: number of reads at the given position
1
0
Entering edit mode
9 months ago
wiki97 • 0

Hi, I am new to exome-seq and would be grateful for any suggestions :)

I want to check coverage for variants called with MuTect2 (GATK 4.1.8.1).

I would like to filter variants basing on the number of reads at the given position and the number of reads in which alternate alleles have been called.

However, I am not sure which annotation should I analyze. I found "AS_SB_TABLE" annotation, and I am wondering - does it show total read counts for reference and alternate alleles? Or are the reads number somehow normalized for the strand bias tests?


Here is the header for "AS_SB_TABLE":

##INFO=<ID=AS_SB_TABLE,Number=1,Type=String,Description="Allele-specific forward/reverse read counts for **strand bias tests**. Includes the reference and alleles separated by |.">
MuTect2 GATK coverage • 372 views
ADD COMMENT
1
Entering edit mode
9 months ago

number of reads at the given position

look at FORMAT/DP (or INFO/DP)

number of reads in which alternate alleles have been called.

see FORMAT/AD

ADD COMMENT
0
Entering edit mode

Dear Pierre Lindenbaum, thank you very much for your answer!!

ADD REPLY

Login before adding your answer.

Traffic: 1261 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6