Hi all,

I'm interested in learning more about best practices for merging RNA-seq technical replicates. I've read many Biostars posts on the matter, but I have a somewhat special case.

**Background:**

I sent RNA samples for sequencing, which were split across 4 lanes each as per normal practice. However, the sequencing depth for this run was much lower than expected and the sequencing core re-sequenced the same samples. Thus, I now essentially have 8x technical replicates per sample, 4x lanes from each of two runs. The repeat runs achieved much better depth.

**Questions:**

1. In this case, should the lower depth data be discarded, or can these data still be used in combination with the updated runs? 
2. If combined, is there a need to mitigate batch effect? Aside from read depth, the samples have nearly identical statistics with respect to FastQC analysis, % genome alignment individually with STAR, etc.
3. For combining, what stage is most appropriate? Aside from file sizes, is there any difference between merged .fastq files and merged .bam files? What about at the level of raw counts? In the past, I have merged .bam files from different lanes and found that the effect was summing the raw reads per gene between replicates.
It's very very unlikely that the library prep was redone.  You almost certainly should just put everything together; **running the same library on different lanes or different days does not introduce any technical artifacts.**

The only exception would be if the first run had a serious technical problem with the instrument itself causing the reads to be totally unusable, but I doubt they would have given you the reads at all if that were the case.