Hi,
I work on plant species. I would like to get flanking sequences of SNPs from the reference genome and the different bam files.
Any help would be appreciated.
Thanks
•
written
by
myprogramming2016 •
4 answers
How can I get snp flanking sequence from the genome file?
Thanks
Do you have a VCF file listing all the SNPs in the reference? If so, you can use the SNiPlay application: http://sniplay.southgreen.fr/cgi-bin/analysis_v3.cgi Load your VCF file and select the plant genome used as reference. And then export flanking sequences of SNPs...
You could use the FastaVariant class in pyfaidx:
To get read pairs surrounding a SNP site from a BAM, you might try VariantBam. You can input any number of sites, and add a padding to the SNP if you like: https://github.com/jwalabroad/VariantBam
Log in to answer this question.
-
how to map every SNP flanking sequences to a reference genome?
written by blacktomato27
-
Extract flanking regions from bam files based on coverage threshold
written by laura.bertola •
-
retrieve multiple bacterial genes and flanking regions
written by benys •
-
XP-CLR on plant data
written by yp •
-
Extract SNPs flanking sequences based on VCF and genome Fasta files
written by Denis
-
Extract sequences flanking SNPs from organism without reference genome
written by konanamenanjacky •
-
RNAseq reads mapped to the gene
written by myprogramming2016 •
-
retrieve orthologous sequences from human fasta file
written by jfertaj
-
SNP calling with fasta files
written by joycewang83 •
-
Web resources for downloading intergenic regions from many plant species?
written by biolab
So you have bam files and a reference genome, and a list of snps (their positions), right? Or do you still need to generate the SNP calls based on the bam files? Do you want the flanking sequences as in the reference genome or as sequenced for your bam files?
Thanks for your reply. I do have SNPs and their positions. I would like to get the flanking sequences in reference as well as in bams.
Thanks
How can I get snp flanking sequence from the genome file?
Thanks
I'm not sure it makes sense to get the flanking sequence from the BAM files. This would be either the flanking sequence from all the reads at all the variant positions, or a consensus variant call which you presumably have in your VCF/FASTA combination.