This is a test version of Biostars. For the public version, visit https://www.biostars.org.
updating SNP names in VCF file

Hello How can I update SNP names from position (chr1:xxxxxx) to rs name in VCF file ? I have a reference VCF file with rs names and want to update my own data VCF file.

Thanks for your help

snp rs names vcf
• 3,950 views • link
updated by Jimbou Jimbou • written by mel22 mel22

intersect sample vcf with reference vcf with bedtools or annotate with bcftools

• 0 views • link
written by cpad0112 cpad0112

2 answers

dbSNP annotation can be added in the ID filed using the VariantAnnotator tool from GATK.

gatk  --java-options '-Xmx4G' \
   VariantAnnotator \
   -R /userdata/arup/hpcc_test/genome/hg38/hg38.fa \
   -V input.vcf.gz \
   -O output.vcf.gz \
   --dbsnp dbsnp.vcf.gz
• 0 views • link
written by Arup Ghosh Arup Ghosh

Thank you very much

• 0 views • link
written by mel22 mel22

you can also use SnpSift from the SnpEff package http://snpeff.sourceforge.net/SnpSift.html#annotate

 java -jar SnpSift.jar annotate dbSnp132.vcf input.vcf > output_annotated.vcf

check ftp://ftp.ncbi.nih.gov/snp/organisms/ for other dbSNP builds

• 0 views • link
written by Jimbou Jimbou

Thank you very much

• 0 views • link
written by mel22 mel22

Log in to answer this question.

More posts like this