I have some RNA-seq samples from mouse, 2 conditions with 4 replicates each, read quality is good and for each sample 85-90% of reads align. The number of aligned reads in millions are:
Condition 1 replicate 1: 100
Condition 1 replicate 2: 79
Condition 1 replicate 3: 52
Condition 1 replicate 4: 37
Condition 2 replicate 1: 56
Condition 2 replicate 2: 59
Condition 2 replicate 3: 31
Condition 2 replicate 4: 26
I know DESeq2 normalizes based on the number of reads but I'm wondering if the variability here will cause any problems?
Thanks that's good to know.
>This looks like a reasonable, expected variation in total read counts.
We sequenced on a NextSeq and expected the reads evenly distributed between the 8 samples and I was surprised by the variability, I expected some but not this much. I was planning to follow up with the sequencing facility because I was a bit concerned that the variability could be due to errors in quantification or poor quality control by them (we submit RNA and they do the rest). But you think the numbers seen here are normal?