This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Variant Calling For Ultra High Coverage Region Of Whole Exome Sequencing Data

Hi,

How does everyone solve the variant calling (SNVs, Indels) problem in ultra high coverage region of whole exome sequencing data? For example, if the ref allele count is 278 and the alt allele count is 83, how can we determine if this site is a variant, or not?

Thanks.

variant calling

2 answers

If want to make the calls by your self, you could likely calculate the alternate allele fraction for known sites intersecting with dbSNP, 1K genomes, ESP project etc and create a cut-off based on that.

I don't know if < 400 counts as ultra-high any more, but most variant calling software will take other things into account such as the quality distributions of the reference bases compared to those of the alternate, strand bias, etc.

That coverage is low enough that if you get a smallish number of candidates, you can load the bam(s) up in IGV and have a look. That usually helps to determine if there are alignment issues in that region or if things look off.

Log in to answer this question.