This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Vcftools Vcf-Stats Field Definitions

Does anyone know the definitions of the fields that vcftools vcf-stats returns:

  • hom_AA = ?
  • het_RA =?
  • snp_count = total count SNPs
  • count = SNP+missing+Ref?
  • hom_RR_count = ?
  • ref = reference calls?
  • missing = missing calls
  • private = private calls not found in other samples
  • het_AA_count = ?
  • ref_count = ? reference calls (same as ref?)
  • unphased = unphased snps

The AA, RA and RR confuse me. I thought they might stand for reference and alternative but why is there a hom_RR count, het_AA_count, ref and ref_count field?

Or does anyone know another tools that also provides the total SNP counts per strains and the private SNP counts per strain?

vcf statistics vcftools
• 7,707 views • link
updated by dangenet • dangenet • written by William William

1 answer

hom_AA = homozygous for a single alternate allele (eg. both alleles have the same mutation).

het_AA = both alleles are non-reference but they are not the same allele (e.g. one has the S98A mutation and the other has the L206P mutation). Think of it as het_A1A2 if that helps.

hom_RR = homozygous reference

het_RA = one reference allele, one alternate allele

"counts" appear to be raw counts at a given locus. The documentation isn't totally clear.

EDIT: My vcf-stats is broken ATM, so I'm having trouble getting more details.

• 266 views • link
written by dangenet • dangenet

Log in to answer this question.

More posts like this