Cell ranger reads do not match Seurat reads
1
0
Entering edit mode
2.9 years ago
firestar ★ 1.6k

This is the summary from the CellRanger run showing ~150,000 reads per cell.

cellranger-report

But, when I read the results into R, I get only ~5000 reads per cell. This is before any pre-filtering in Seurat.

x <- Seurat::Read10X_h5("data/PBMC_1_3_1/count/sample_feature_bc_matrix.h5")$`Gene Expression`
hist(colSums(x))

enter image description here

median(colSums(x))
[1] 5050.5

This is a huge difference and I wonder if I am making a mistake somewhere.

10x cellranger seurat • 1.4k views
ADD COMMENT
1
Entering edit mode
2.9 years ago

Your Seurat histogram is actually of unique UMIs per cell, not raw read counts per cell. If you look at the cell ranger results under "Median UMI counts per cell" it matches (almost exactly) to your results from Seurat.

ADD COMMENT
0
Entering edit mode

Ah right! I knew I was overlooking something obvious. Thanks!

ADD REPLY
1
Entering edit mode

Yeah, 10xGenomics cheats a little by displaying the total reads far more prominently than the value you really want, the # of UMIs.

ADD REPLY

Login before adding your answer.

Traffic: 890 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