freebayes calls presumably homozygous variants as heterozygotes
Entering edit mode
7 months ago

Hi Everybody,

I ran freebayes to detect variants in multiple samples and now I saw that sometimes at a site there is only one reference observation and 10 or more alternate observations but still the sample is called heterozygous, even though it is more likely to be homozygous for the alternative allele in fact. Often also the reference observation has quite low quality in comparison to the alternative observation.

I used the --min-alternate-total and --min-alternate-fraction parameters to set a threshhold for when to consider an alternative allele, but is there a way to set a threshhold like this for the reference allele as well or any other way to reduce this behaviour?

SNP freebayes variant-call genotyping • 228 views

Login before adding your answer.

Traffic: 1759 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6