Hi ! I did scRNAseq (10XChromium) of cancer cells and now I am looking for a simple way to split my FASTQ Files into celll specific FASTQ Files. Alternativley split my BAM File into cell specific BAM Files - because I need one FASTQ (BAM) File per cell in order to do perform variant calling on the single cell level. My goal is to obtain one single VCF file per cell.
Is there a build-in function in CellRanger for that? Are there any other tools to demultiplex my FASTQ/BAM Files using cell barcodes and split them into cell specific FASTQ/BAM Files?
Any help is highly appreciated! Thank you so much!
•
written
by
niklas.lang •
2 answers
If you want to call variants from 10x data, there is a specialized tool by 10x that is optimized for their data and uses the combined BAM file (no need to modify your BAM file): https://github.com/10XGenomics/vartrix
If you are still interested in splitting a BAM file by barcode, there are some approaches outlined in this blog post: https://divingintogeneticsandgenomics.rbind.io/post/split-a-10xscatac-bam-file-by-cluster/
Log in to answer this question.
-
Change sample ID in BAM file to cell barcode
written by martin.grasshoff •
-
What is the barcode and UMI in my 10x single cell sequence data?
written by Hmm •
-
5' single-cell RNA-seq Demultiplexing
written by learner-MD •
-
how to demultiplex scRNAseq fastq with barcode in read2
written by Binghong •
-
Demultiplexing 10x fastq files by GROUPS of barcodes
written by bds1217 •
-
Separating each cell in fastq file
written by sidrah.maryam •
-
From 10x single cell data to a bulk analysis
written by enho •
-
How to find percentage of a gene (eg.Ppbp) for condition A and condition B from 10X Single Cell RNA…
written by bioinforesearchquestions
-
Creating a cell by SNP matrix
written by dominicdhall •
-
Analyzing MARS-seq (single-cell RNA seq) fastq files
written by Eldad Shulman •
It is incredibly unlikely that this will prove useful. 10X reads will be enriched around the polyA site and be highly redundant.
My idea was to split the BAM File into several cell-specific BAM File by using the "CB" Tag. I don't know how this could possibly interfer the PolyA enrichment, so I would be glad to hear what your concerns are?
The question is whether single cell variant calling with this type of data is likely to be useful. The answer is "no", because you're mostly sequencing around the TES.
...meaning you're limited to detecting (exonic) variants close to the TES