Question: How to properly map 4 lanes for differential expression using BWA?
gravatar for satshil.r
3.7 years ago by
United States
satshil.r50 wrote:

Hey guys,

I have 8 samples, each sample having 4 lanes. How do I go about mapping the sample to my reference when there's multiple lanes?

I'm not sure how cating the 4 lanes together would work out for counts. Here's my thought process:

If I cat the 4 lanes together before mapping (BWA), and then run BWA mem, won't that give me 4x the expression/counts since the mapping consists of 4 duplicate reads?

If this is actually what would happen, or do I have the process all wrong?

How would you properly map/align these 4-lane reads?

ADD COMMENTlink modified 3.7 years ago by Devon Ryan88k • written 3.7 years ago by satshil.r50
gravatar for Devon Ryan
3.7 years ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k wrote:

Splitting a library over 4 lanes doesn't result in 4 copies of every read, no one would spend the extra money to run the additional 3 lanes then. Just concatenate the files and map them, though I wouldn't recommend bwa for RNAseq (it's not splice aware, try STAR or hisat or even BBMap, which quite flexible).

ADD COMMENTlink written 3.7 years ago by Devon Ryan88k
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: 1612 users visited in the last hour