Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Why is there a big difference between kallisto and cellranger count result? #50

Open
aizhimin opened this issue Feb 16, 2023 · 6 comments

Comments

@aizhimin
Copy link

sample:
https://www.ncbi.nlm.nih.gov/sra/SRR8315735

kallisto result:
image

cellranger result:
image

@gaofan83
Copy link
Collaborator

The data you analyzed only has 25 million reads (people typically sequence 10X more reads). How many "cells" did you recover from the two pipelines?

@aizhimin
Copy link
Author

@gaofan83 About 160 cells.

But why is there a big difference between kallisto and cellranger?

@gaofan83
Copy link
Collaborator

gaofan83 commented Feb 17, 2023

@aizhimin 160 cells from kallisto? How about cellranger?
My understanding is that the dataset you analyzed is a multi-omics data with both 5' expression and immune-seq. Are you sure you chose the correct 5' barcode whitelist as the input for kallisto?

@aizhimin
Copy link
Author

@gaofan83 This is my cammand:
kb count -i /mnt/data/ref/transcriptome.idx -g /mnt/data/ref/transcripts_to_genes.txt -t 8 -m 16G --mm --h5ad --filter bustools --cellranger --gene-names --verbose --overwrite -x 10XV2 -o output/SRR8315735 SRR8315735_1.fastq.gz SRR8315735_2.fastq.gz

Where is the the correct 5' barcode whitelist ?

@gaofan83
Copy link
Collaborator

I am not familiar with kb. If you use kallisto/bustools, you will see there is an option for the whitelist. You can search 10X site and see whether 3' barcode whitelist (people typically call as 10xv2 or 10xv3) is the same as 5' barcode whitelist.
https://github.com/BUStools/bustools

@aizhimin
Copy link
Author

@gaofan83 10xv2 Contains 5' barcode whitelist. Look at this image from 10x site。

image

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants