How can I count the number of reads that support a variant in a vcf file?

I have seen some methods to count the total number of SNPs, or the total SNPs at a site, but how can I count the reads that support each variant at each site from either vcf, bam, or sam?

Hi @iaincambeul,
if you are interested in specific positions, this tool seems to suit your requirements: BAM-readcount.


Depending on what you’re using for variant calling, some of the keys in the INFO field of the VCF might give you what you want. e.g. AO is the count of the number of reads supporting a particular variant allele.

1 Like