Question: Assigning gene identifiers to SNPs with bcftools annotate
0
gravatar for maxa
6 weeks 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 • 78 views
ADD COMMENTlink modified 6 weeks ago • written 6 weeks ago by maxa0

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

ADD REPLYlink written 6 weeks ago by Pierre Lindenbaum129k
0
gravatar for maxa
6 weeks ago by
maxa0
maxa0 wrote:

I found the solution in another post here.

ADD COMMENTlink written 6 weeks 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: 1702 users visited in the last hour