Question: Proper use of BWA-backtrack and BWA-MEM
0
gravatar for cfarmeri
13 months ago by
cfarmeri130
Japan
cfarmeri130 wrote:

Hi, Biostars.

I have a set of Illumina GA2 paired-end fastq files.

These reads of the fastqs has 76b length.

Before mapping, I trimmed barcode sequence from 5' end of both 1st and 2nd reads, finally, the reads have 63-67bp length.

BWA's author says following.

・reads < 70bp should be mapped by BWA-backtrack (BWA-backtrack can map 70bp< and <100bp reads, but they mention than these reads should be BWA-MEM algorithm.)

・reads > 100bp should be mapped by BWA-MEM. BWA-MEM algorithm is recommeded.

Any one has an idea which BWA-backtrack or BWA-MEM I should use?? Thanks!!!

genome • 965 views
ADD COMMENTlink modified 13 months ago by dariober8.0k • written 13 months ago by cfarmeri130
1
gravatar for Brian Bushnell
13 months ago by
Walnut Creek, USA
Brian Bushnell13k wrote:

If the reads are < 70 bp after trimming, which it sounds like they are, then the official advice is to use BWA-backtrack. It does not matter how long they were initially.

Personally, if I were to use a BWA variant, I would choose BWA-mem over BWA-backtrack because it is so much faster and, in my testing, more accurate. But I have not tested it on < 70bp reads.

ADD COMMENTlink written 13 months ago by Brian Bushnell13k

Thanks, Brian.

I confused the 63-67bp is a value close to 70bp. First of all, I try BWA-backtrack algorithm.

ADD REPLYlink written 13 months ago by cfarmeri130
1
gravatar for dariober
13 months ago by
dariober8.0k
Glasgow - UK
dariober8.0k wrote:

Hi- Have a look also at this post When and why is bwa aln better then bwa mem?

Personally, these days I use bwa mem even when read length is < 70bp. I haven't tested it properly, but I find that for shorter reads, especially if paired-end, bwa mem becomes slower and more memory hungry than bwa aln.

ADD COMMENTlink written 13 months ago by dariober8.0k
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: 590 users visited in the last hour