Question: bwa MEM for illumina 1.3+
0
gravatar for Zev.Kronenberg
3.3 years ago by
United States
Zev.Kronenberg11k wrote:

I see that BWA ALN has an "-I" flag for changing the quality of the reads.  I don't see this option for BWA MEM? Does anyone know if I'm missing something?

phred bwa mapping quality • 1.3k views
ADD COMMENTlink modified 3.1 years ago by Brian Bushnell15k • written 3.3 years ago by Zev.Kronenberg11k
2

You're not missing anything: Bwa Mem: How To Specify The Fastq Phred Format? (see my link to a bwa email list message about this from Heng Li).

ADD REPLYlink written 3.3 years ago by Devon Ryan74k
1

what is the size of your reads ? See Heng's comment:Aligning Sequence Reads, Clone Sequences And Assembly Contigs With Bwa-Mem

BWA-MEM however lacks the following features:

  • Working with very short reads. The performance and the accuracy of BWA-MEM will degrade.
  • Guaranteed sensitivity to hits within a certain edit distance threshold (as with BWA and GEM).
ADD REPLYlink written 3.3 years ago by Pierre Lindenbaum102k

100bp.  I wanted to use mem for consistency across some simulations and real data.

ADD REPLYlink written 3.3 years ago by Zev.Kronenberg11k

Honestly, I was stuck between mem and aln for my 100bp data, and aln performed much better, more reads aligned.

ADD REPLYlink written 3.3 years ago by Adrian Pelin2.0k
1
gravatar for Brian Bushnell
3.1 years ago by
Walnut Creek, USA
Brian Bushnell15k wrote:

You can use my reformat tool to change ASCII-64 qualities to ASCII-33 (and most other format conversions).  It supports piping from stdin and to stdout, so you don't have to write a temp file.  And it's the fastest format convert of which I am aware.

ADD COMMENTlink written 3.1 years ago by Brian Bushnell15k
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: 647 users visited in the last hour