This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Recalculate QUAL and INFO fields for a subset of samples in a VCF

Dear All,

I extracted some samples from a larger VCF (one population of a species). I would like to filter the sites by quality and depth. However, the values in the new VCF still seem to reflect the data from the bigger VCF (i.e. the average depth was calculated for all samples, not just the extracted ones). Do any of the usual tools (GATK, bcftools, etc) offer a way of recalculating site quality or the INFO fields for the current subset of samples?

Many thanks

snp next-gen
• 3,697 views • link
updated by GouthamAtla GouthamAtla • written by hermathena • hermathena

Besides the issue of "software that can do that", there is the issue of whether the data to do so are supplied with the VCF samples. What format fields are present for the genotype records?

• 0 views • link
updated by Ram Ram • written by Sean Davis Sean Davis

If you have GL/PL, some versions of bcftools can do this.

• 0 views • link
written by lh3 lh3

Indeed, the FORMAT fields available are AB, AD, DP, GQ, GT, MQ0, PL.

Which version of bcftools, please?

• 0 views • link
written by hermathena • hermathena

1 answer

Enrichr, looks updated and queries many databases.

http://amp.pharm.mssm.edu/Enrichr/

• 0 views • link
written by GouthamAtla GouthamAtla

Could you elaborate on how this tool would help answer the original question?

• 0 views • link
written by Sean Davis Sean Davis

I answered a different question, by mistake.

• 0 views • link
written by GouthamAtla GouthamAtla

Log in to answer this question.

More posts like this