Question: Comparing vcf with dbSNP
0
gravatar for cvu
5.7 years ago by
cvu150
India
cvu150 wrote:

I've generated vcf file using samtools. now how do i compare my SNPs with dbSNP ?

alignment snp next-gen assembly • 4.5k views
ADD COMMENTlink modified 5.7 years ago by cristianrohr76830 • written 5.7 years ago by cvu150
3
gravatar for poisonAlien
5.7 years ago by
poisonAlien2.8k
Asgard
poisonAlien2.8k wrote:

Assuming you have dbSNP in VCF format you can use vcf-compare

ADD COMMENTlink modified 4 months ago by RamRS25k • written 5.7 years ago by poisonAlien2.8k

can we use samtools or GATK for this?

ADD REPLYlink written 5.7 years ago by cvu150

Samtools and GATK can't be used to add dbSNP rsID information to vcf file. Either you use vcf-compare as suggested above or your own script to add rsids to the vcf file. You will have to download dbSNP data from NCBI.

Check these links:

Snp Rsid By Location?

http://www.ncbi.nlm.nih.gov/books/NBK9777/#FTP.how_do_i_extract_information_from_db

ADD REPLYlink modified 5.7 years ago • written 5.7 years ago by Ashutosh Pandey12k

You can add dbSNP annotations with GATK using VariantAnnotator with the --dbsnp flag:

ADD REPLYlink modified 4 months ago by RamRS25k • written 5.7 years ago by Brad Chapman9.5k
2
gravatar for cristianrohr768
5.7 years ago by
Spain
cristianrohr76830 wrote:

You can annotate with dbSNP ids using SnpSift, and then with the effects of the variants using SnpEff.

ADD COMMENTlink modified 4 months ago by RamRS25k • written 5.7 years ago by cristianrohr76830
1
gravatar for Bert Overduin
5.7 years ago by
Bert Overduin3.6k
Edinburgh Genomics, The University of Edinburgh
Bert Overduin3.6k wrote:

If your goal is to find out whether your SNPs are known or novel, you can also think of using the Ensembl Variant Effect Predictor (VEP).

ADD COMMENTlink written 5.7 years ago by Bert Overduin3.6k
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: 1784 users visited in the last hour