Question: Assigning gene identifiers to SNPs with bcftools annotate
0
gravatar for maxa
5 months ago by
maxa0
maxa0 wrote:

Hello,

I have a VCF file with SNPs. It contains the chromosome and the position of the SNPs, but the ID column is empty (.). I want to fill the ID column with gene identifiers which i have in a tab-deliminated file with four columns: Chromosome, start position, end position, identifier. Let's call it coordinates.gff.

I used

bcftools annotate -a coordinates.gff.bgz -c CHROM,FROM,TO,ID -o variants_annotated.vcf -Ov --threads 2 variants.bcf

A tabix index exists for the coordinates file and there are no error messages. Still, the ID column of the output file is empty. I'm struggling with this for a while now and I have the feeling that I'm missing something very obvious. Do I have to specifiy that I want the IDs to be replaced in some way?

Any help would be highly appreciated!

snp • 167 views
ADD COMMENTlink modified 5 weeks ago by Biostar ♦♦ 20 • written 5 months ago by maxa0

same chromosome notation ? no +1 problem with the positions ?

ADD REPLYlink written 5 months ago by Pierre Lindenbaum131k
0
gravatar for maxa
5 months ago by
maxa0
maxa0 wrote:

I found the solution in another post here.

ADD COMMENTlink written 5 months ago by maxa0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 909 users visited in the last hour