This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to map every SNP flanking sequences to a reference genome?

Hi to all Good Afternoon I selected some 1000 snps based on MAF, PIC and their coverage in genome based on physical distance. Now i would like to know where and how many places are my SNP flanking sequences match or hit in reference genome? here my goal is to further select snp with only 1 hit, it means it can bind at only one place in genome and will be good results in genotyping. Can i any help me how can i proceed in this process? i did normal NCBI blast but i have more snps to match and it is time taking. Is there any tool or script to make it fast and simple.? any help in this regard is highly appreciated. This is my example snp flanking sequence with [ ] brackets ATGCAATGCATGACGTGCATGCATGCATGCAGTCGCAGTCGTACATGCAGTA and then bracket Thanks in advance

snp alignment sequence blast

You can use any of the standard NGS aligners to map reads in fasta format to a reference genome such as bowtie2.

0 answers

No answers yet.

Log in to answer this question.