Primer design from bam file ?
Hi Biostars,
I have 5 samples that I mapped against a genome.
My goal is to design set of primers (final amplicon size around 100bp) from this, that will allow to amplify regions to look for SNPs between individuals.
Do you know a software/method that take bam files, looking for conserved regions with SNPs inside (that are present in my 5 samples) and produce set a primers ?
Thanks a lot for your help.
• 3,840 views
•
link
updated
by
WouterDeCoster
•
written
by
Picasa
•
written
by
Picasa
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Finding out if a query sequence is present in a read library / bam file (blast against these?)
written by chronotope •
-
How to find conserved regions in a genome
written by Nicolas •
-
Adding flanking regions to a gene for primer design
written by eli_bayat •
-
How to filter .vcf based on .gbk file to remove SNP calls in non-CDS regions?
written by goatsrunfaster •
-
Finding non-polymorphic region of genome for PCR primer design
written by joneill4x
-
How to extract the same mapped region by multiple samples ?
written by Picasa
-
How to design primers when I only know the position of the mutation given in the VCF file?
written by jpxaviersantos90 •
-
primer design to discriminate sequences in intron-exon MSA
written by Anand Rao
-
DNA multialignment and conservation profile
written by Kevin D •
-
primer design using BLASTn
written by J.F.Jiang
Can you tell the purpose of primers ?
To amplify regions and looking for SNPs between individuals
You should probably just call SNPs in the samples and then see which SNPs could distinguish/identify these individuals. Not everything has to (or should be) automated, of course. Once you know the SNPs that distinguish/identify them, you can then design primers accordingly.
I wrote a SOP for designing primers a few years ago, from which you may get a few tips: Designing a single set of primers and probe for a genomic region of interest
Just a suggestion: First, create consensus sequence from bam files for each sample (Generate Consensus sequence from BAM file ). Design primers to your reference genome sequence (Primer3, https://www.ncbi.nlm.nih.gov/tools/primer-blast/ ), and then mach them to your consensuses, to see if they mach correctly (Primer Map (http://www.bioinformatics.org/sms2/primer_map.html ). Best, Agata
Hey Agata, please check your URLs / links. There is a bug in the system whereby, if a parenthesis appears after a URL, it is interpreted as part of it. You just need to add a space between these.
Corrected, thanks! Agata
Thanks Agata, but not sure to understand what the consensus sequence is ?
Consensus sequence is a linear sequence from bam file (mapped reads). For example, if you have sample reads mapped to reference genome, you can generate the linear representation of sequence for this sample. It will include all differ variants (SNPs).
In case you don't want all SNPs, but only the unique for each sample, you can compare consensus linear sequences from 5 samples, and then design primers in desire places. Agata
Suggestion: if by "mapping file" you mean "bam file" then just write the latter. Explicit is better than implicit. I have adapted your title.