Question: BWA inexact matching
0
gravatar for chirag.agrawal.91
9 months ago by
chirag.agrawal.910 wrote:

Hi, I have been trying to understand BWA algorithm. I tried running the tutorial as in -

https://icb.med.cornell.edu/wiki/index.php/Elementolab/BWA_tutorial

bwa aln -t 4 hg19bwaidx s_3_sequence.txt.gz >  s_3_sequence.txt.bwa

bwa samse hg19bwaidx s_3_sequence.txt.bwa s_3_sequence.txt.gz > s_3_sequence.txt.sam

I think these steps only give me exact matches for the queries. Can someone guide me as to how to run for inexact matches. How to set the threshold for inexact matches?

bwa alignment • 288 views
ADD COMMENTlink modified 9 months ago by genomax69k • written 9 months ago by chirag.agrawal.910

Take a look at the `bwa aln options on the help page.

Can someone guide me as to how to run for inexact matches

I assume that means you want to allow for mismatches in alignment?

ADD REPLYlink written 9 months ago by genomax69k

Hello,

to understand the algorithm by BWA I guess you have to first take a look at the code :)

BTW: I think there is no reason to use bwa aln. Use bwa mem instead.

fin swimmer

ADD REPLYlink written 9 months ago by finswimmer11k
Please log in to add an answer.

Help
Access

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