Question: bwa mem -T (alignment score) not doing anything
gravatar for chris.bird
2.6 years ago by
chris.bird0 wrote:

I'm trying to control which reads are mapped to a reference genome by adjusting the the -T parameter in bwa mem, which is the minimum acceptable alignment score. The only problem is that reads with an alignment score below the threshold are being mapped anyway. I can see this when using samtools view to visualize the bam file, reads with an AS value less than the threshold are there. My data is paired end, so I tried doubling the threshold to account for both reads, but that didn't work either.

What's going on here? Do some settings in bwa mem interfere with -T?

snp sequence genome • 925 views
ADD COMMENTlink written 2.6 years ago by chris.bird0

For things like this, please give the version of BWA, full command line, the reference genome and an example of the problematic reads.

ADD REPLYlink written 2.6 years ago by ATpoint44k
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: 1934 users visited in the last hour