Using dbSNPs for mapping

kellermackellermac Baton Rouge, LouisianaMember

Hi. So I want to take a dbSNP file and highlight SNPs that are found in some BAM files. Instead of filtering these I want to search for these SNPs only. I will then use these SNP positions as markers to trace recombination events that have occured between some strains I am studying. So far I havn't found a good strategy to accomplish this. I was hoping someone over here knows how to parse the dbSNP, and could guide me through it. Alternatively you could use this idea to generate a new tool! Please let me know if you have thoughts, or questions. Thanks! -Keller



  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    Hi Keller, I'm not sure I understand exactly what you want to do, but in general it sounds like something that should be possible with GATK. For example if you have a VCF file containing variants of interest, and you just want to select those that are present in dbSNP, that is very easy using the GATK tool SelectVariants.

    Can you please clarify what data you have available and exactly what you would like to produce?

Sign In or Register to comment.