Outputting DP4 flag in GATK
1
0
Entering edit mode
6.5 years ago
Jautis ▴ 510

Hello, is there any way to output DP4 (as well as DP) while using GATK's Haplotype Caller to call variants?

Thanks!

variant calling sequencing gatk • 1.7k views
ADD COMMENT
2
Entering edit mode
6.5 years ago

Annotate the variants for StrandAlleleCountsBySample using VariantAnnotator tool.

ADD COMMENT
0
Entering edit mode

So that appears similar to DP4, but I need it in DP4 rather than the set of 8 values StrandAllelesCountsBySample would give. I'm piping the vcf into another program and I'm not familiar enough with how that one runs to modify it to use SAC instead.

ADD REPLY
0
Entering edit mode

I am not sure if you can do that directly. However, you have all the required info, and you can generate an intermediate vcf file (or pipe it directly) using shell-commands.

ADD REPLY

Login before adding your answer.

Traffic: 1003 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