I ran the sample twice on a WES platform, and there is about 3% discordance between the 2 VCFs. Most of the discordant "variants" seem to be homozygous Refs in the databases, except that the reference genome has the other allele as Ref. It seems to be a problem with the ref genome (hg19). Is there a way to filter out such variants?
Thanks
1 answer
bcftools norm in https://samtools.github.io/bcftools/bcftools.html
-c, --check-ref e|w|x|s
what to do when incorrect or missing REF allele is encountered: exit (e), warn (w), exclude (x), or set/fix (s) bad sites. The w option can be combined with x and s. Note that s can swap alleles and will update genotypes (GT) and AC counts, but will not attempt to fix PL or other fields.
Log in to answer this question.
-
Converting custom VCF to standard VCF
written by avelarbio46 •
-
SNP concordance between NGS and SNParray
written by emma.a •
-
A few questions about variant calling not clear to me
written by louisber2765 •
-
Standard approach for telling if sites were removed during calling due to poor quality or non-varia…
written by curious
-
Why is 0/0 homozygous reference call in VCF useful information?
written by lillian.thistlethwaite •
-
Should I replace heterozygous variants to make a consensus
written by deepti1rao •
-
allelic unbalance for multiplex PCR based amplicon data
written by J.F.Jiang
-
Is there a phasing algorithm or tool that first determines undeniably present haplotypes?
written by anp375
-
getting sequence by combining reference and variants (fasta and vcf files)
written by Lin •
-
Vcf Format And Filtering
written by Ashutosh Pandey
