Question: BWA aligner parameter tweaking to show primer mapping
gravatar for morovatunc
15 months ago by
morovatunc400 wrote:

Hi guys,

I generated some primers by mutating amino acids in a DNA sequence of 1553 bases.. My primer structure is as follows;

19bp sequence(flanking sequence%100 match)+ 3 bp (mutated) + 19 bpsequence(flanking sequence %100 match)

Say i have some sequence that start with start codon. Following amino acids are mutated into other amino acids by changing codons. To show that I did this correct. I thought mapping to a reference (1553 bp DNA sequence) will show me those primers.

I followed these steps;

bwa index ref.fa
bwa aln tolc.fa primers.FASTA | bwa samse ref.fa - primers.FASTA | samtools view -bu - | samtools sort - > aln_sorted_output.bam && samtools index aln_sorted_output.bam

I am aware that I am swimming against stream by trying to decrease the specificity of that algorithm but couldnt think of a better way.

So my questions are;

1) Is my method practical ? if not how would you achieve this ? 2) (If my method is ok) What options in bwa should I tweak. I tried to lower the scores of gap penalty and mismatch penalties but did not work well.

Thank you very much for your help,

Best regards,


bwa primer-design • 628 views
ADD COMMENTlink modified 15 months ago • written 15 months ago by morovatunc400
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 691 users visited in the last hour