Greetings.
I am doing some SNP-/variant-calling, mapping to a somewhat divergent reference genome (probably about 1.5-2% different), but synteny is preserved.
I have been fine-tuning my SNP-calling using SAMtools/bcftools, however I am still getting false negatives in regions of high divergence (~10-20% difference) and (often, but not always) when two SNPs are next to each other. It looks like the caller is probably inserting gaps in these regions rather than calling them as SNPs, but there are no indels in corresponding regions in the VCF.
Even looking at the unfiltered variant calls, nothing shows up in these regions, and they are the interesting ones.
All help is appreciated. Thanks!
I am doing some SNP-/variant-calling, mapping to a somewhat divergent reference genome (probably about 1.5-2% different), but synteny is preserved.
I have been fine-tuning my SNP-calling using SAMtools/bcftools, however I am still getting false negatives in regions of high divergence (~10-20% difference) and (often, but not always) when two SNPs are next to each other. It looks like the caller is probably inserting gaps in these regions rather than calling them as SNPs, but there are no indels in corresponding regions in the VCF.
Even looking at the unfiltered variant calls, nothing shows up in these regions, and they are the interesting ones.
All help is appreciated. Thanks!
Comment