bwa MEM for illumina 1.3+
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?
• 3,872 views
•
link
updated
by
Ram
•
written
by
Zev.Kronenberg
•
written
by
Zev.Kronenberg
1 answer
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.
• 0 views
•
link
updated
by
Ram
•
written
by
Brian Bushnell
•
written
by
Brian Bushnell
Log in to answer this question.
More posts like this
-
BWA-MEM or BWA-ALN + BWA-SAMSE for MGI sequencing
written by khanhlpbao •
-
Fastq screen: use bwa aln and not bwa mem
written by Matteo •
-
Why ExAC used bwa aln instead of mem?
written by abascalfederico
-
bwa aln for paired-end sequencing data in one .fastq file
written by ciemanek
-
BWA-MEM Vs BWA-ALN
written by SOHAIL
-
Bwa mem is giving me alignements that are causing noise in subsequent analysis. How I can filter t…
written by shinken123
-
does bwa mem use base qualities?
written by Richard
-
bwa mem on a bam file
written by vimartin •
-
When and why is bwa aln better then bwa mem?
written by dariober
-
Bwa-Mem And Quality Trimming
written by nupurgupta0806 •
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:
100bp. I wanted to use mem for consistency across some simulations and real data.
Honestly, I was stuck between mem and aln for my 100bp data, and aln performed much better, more reads aligned.
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).