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 to get a variant's abundance from calling file?
I'm wondering how to get a mutation's abundance based on variant calling result file. I noticed there is a 'AF' field in vcf file, but the value is always 1 or 0.5. It's different from the ratio of ref/alt reads in bam file. I also found that 'AD' filed in FORMAT column records the ref and alt reads count, but it can be different with corresponding reads number in bam file as well. So which field should I use to get variants' abundance? Thank you!