Multiple alleles in REF and ALT in VCF file

Hi all, I am looking at some data in my VCF file and came across this line below:

20      62855516        20:62855516:GC:AC       GC      AC      .       .       PR;AC=8;AN=70   GT      0/0     0/0     0/0     0/0     0/1     0/1     0/1   0/0      0/0     0/0     0/0     0/0     0/0     0/0     0/0     0/0     0/0     0/1     0/0     0/0     0/0     0/0     0/1     0/1     0/1     0/0     0/0   0/0      0/1     0/0     0/0     0/0     0/0     0/0     0/0

My question is, is this considered a multiallelic site? How should I handle this entry? I am also looking to compare sites with another VCF, and the equivalent position in the other VCF has G in REF and A in ALT, is there a way to clean the data to consider only the first allele for this entry?

I've tried to remove/fix these entries with bcftools view -m2 -M2 -v snps and bcftools norm -m -any but it doesn't seem to catch it.

Thanks.

vcf allele indel bcftools vcftools

• 1,369 views • link

updated 5.4 years ago by Biostar

• written 5.5 years ago by Volka

Instead of removing, try decomposing vcfs with Vt. Refer to the decompose biallelic block substitutions section here: https://genome.sph.umich.edu/wiki/Vt

• 0 views • link

written 5.4 years ago by cpad0112

0 answers

No answers yet.

More posts like this

Multiallelic variants when merging VCF's with GLnexus
written 2.6 years ago by Michal Frenkel •

I'm attempting to combine around 140 .g.vcf files into a single file using GLnexus on the DNAnexus platform. To examine multiallelic variants, I'm normalizing the …
Can not get bcftools norm to join biallelics into a multiallelic.
written 3.4 years ago by a615ebfb •

is this the right way to use bcftools to join/merge biallelic records into a multiallelic? If so, it is not working. No errors but it …
How to remove lines from a vcf file
written 3.6 years ago by Khaleesi95 •

Hi guys, I'm really new to bash. I have a vcf file, and these are the first 37 rows of the file: ##fileformat=VCFv4.2 ##FILTER=<ID=PASS,Description="All filters …
How to remove shared SNP from a VCF with multiple individuals
written 3.9 years ago by cuamatzi •

Hi! I have a VCF produced by **MafFilter** with 29 samples, with the next format (trimmed to 5 strains for easier reading): ##fileformat=VCFv4.0 ##fileDate=202291 ##source=Bio++ …
bcftools multiallelic split not working
written 5.0 years ago by from the mountains

I am attempting to split multiallelic sites using `bcftools norm` with the following command: ``` zcat ${inputVcf} | \ sed 's/AD,Number=./AD,Number=R/g' | \ sed 's/ADR,Number=./ADR,Number=R/g' …
What is the correct way of setting the genotype after splitting multi-allelic sites in a VCF file?
written 5.2 years ago by terdon

Consider this minimal VCF file which contains a multi-allelic call made by a somatic variant caller on a human tumor sample: ``` ##fileformat=VCFv4.3 ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> ##FORMAT=<ID=AD,Number=R,Type=Integer,Description="Allelic …
bcftools to remove non-variants sites
written 5.9 years ago by shiy05 •

Dear all, I am used to using vcftools to remove non-variants sites. The command I used: ``` vcftools --vcf input.vcf --mac 1 --recode --out input.mac …
VCF files - how to handle entries with identical ID but different positions and alleles?
written 6.4 years ago by Volka

Hi all, I am working on VCF files on imputed data from the Sanger imputation server, and ran into a problem where there are entries …
Removing both duplicates from a VCF
written 6.6 years ago by john.michel.rouhana •

Hi, I'm aware of ways to merge duplicates in a vcf, like [this one][1]. With my vcf, however, there appear to be duplicate variants with …
bcftools annotate picks one attribute but not another
written 8.7 years ago by Ram

Hello, I've run into a weird error and would like to figure out what's going on. I'm annotating a VCF file from a library VCF …