Question: Mismatch penalty in BWA
0
gravatar for Juan Ledesma
3 months ago by
United Kingdom
Juan Ledesma0 wrote:

Hi

Does anyone know the highest mismatch penalty that I can set to get a more selective alignment in BWA?

I would like to map my reference sequence only to the reads that are identical and exclude those ones that can be similar but no identical to the reference sequence.

I am trying to identify the reads belonging to two genotypes in a sample even if both are very similar.

Thanks

bwa alignment mapping • 163 views
ADD COMMENTlink modified 3 months ago • written 3 months ago by Juan Ledesma0
0
gravatar for Jennifer Hillman Jackson
3 months ago by
United States
Jennifer Hillman Jackson23k wrote:

Hello,

This prior post contains setting for finding hits with no internal mismatches. https://www.biostars.org/p/132681/

The advice to remove reads with unaligned ends (not fully contained within the central alignment) through filtering is a good if you want to do that as well. There are a few tools that can filter BAM datasets - search in the tool panel with the keywords "bam filter". Your target is the CIGAR string.

Line command arguments map directly to those on the Galaxy wrapper's tool form.

Thanks, Jen, Galaxy team

ADD COMMENTlink written 3 months ago by Jennifer Hillman Jackson23k
0
gravatar for Juan Ledesma
3 months ago by
United Kingdom
Juan Ledesma0 wrote:

Thanks for your comments Jen, it is very useful for me. Cheers

Juan

ADD COMMENTlink written 3 months ago by Juan Ledesma0
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 16.09
Traffic: 94 users visited in the last hour