How to properly map 4 lanes for differential expression using BWA?
1
0
Entering edit mode
7.5 years ago
satshil.r ▴ 50

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?

DE RNA-Seq BWA Differential Expression HT-Seq • 1.8k views
ADD COMMENT
0
Entering edit mode
7.5 years ago

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 COMMENT

Login before adding your answer.

Traffic: 2448 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6