Separating a bam file into multiple bam files

How can I separate a bam file into multiple bam files such that each bam file contains information for one cell using UMI-tools?

bam scrna-seq umi-tools

• 1,815 views • link

updated 3.1 years ago by Ram

• written 3.5 years ago by Researcher •

Assuming your data was processed by cellranger there should be a BAM fie with CB tags in it. Using those to split the file would be the way to go.

Some background info: https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/output/bam#bam-bc-tags

There are some approaches outlined in the blog linked by @igor in this thread: scRNAseq: Split FASTQ into cell specific FASTQ Files

• 0 views • link

written 3.5 years ago by GenoMax

1 answer

10x provides an unsupported tool for this: https://github.com/10XGenomics/subset-bam You will need a list of cell barcodes.

• 0 views • link

written 3.5 years ago by GenoMax

Thank you. I appreciate your help. I have not tried your answer yet because I found an R package, Sierra that seems to have a function for splitting a bam file into multiple bam files for every cell, "https://rdrr.io/github/VCCRI/Sierra/man/SplitBam.html". Have you ever tried this?

• 0 views • link

written 3.5 years ago by Researcher •

More posts like this

How to get host and virus UMI count per cell with Bacterial genome data?
written 5 months ago by San •

Hello, I am working with bacterial genome fasta file. I have ran Virsorter2 tool to identify the viral regions. And i have annotated the viral …
How to Split FASTQ File by UMI Indices into Multiple Files
written 2.2 years ago by LDT •

Hello Biostars Community, I have a FASTQ file that contains sequences with UMI indices in the headers. I have already counted the UMI indices and …
counting from a single BAM file with multiple samples
written 2.9 years ago by dr-device •

I've got a BAM file that contains RNA-seq data from about 20 samples. The samples are differentiated by one another via the use of a …
Separate single cell BAM file by the cell barcode
written 3.0 years ago by zbidav •

Dear Biostar community, I am a bit new to Dropseq analysis (10x sequenced files, if not mistaken). I followed the standard CellRanger protocol and received …
Generating vcf files per individual instead of one big file from bam files
written 4.0 years ago by salman_96 •

Hi, I have multiple bam files that I used to generate vcf file for variant calling. I used the command below bcftools mpileup -O v …
Separating each cell in fastq file
written 5.8 years ago by sidrah.maryam •

Hello, I want to analyze single-cell RNA seq data that I took from a paper. The paper mentioned it has 8000 single cells and the …
from BAM to BED passing cell barcodes
written 6.3 years ago by BrunoGiotti

Hello all, It seems that different flavors of this question have been already addressed yet i couldn't find this one in particular. I have been …
scRNAseq: Split FASTQ into cell specific FASTQ Files
written 7.2 years ago by niklas.lang •

Hi ! I did scRNAseq (10XChromium) of cancer cells and now I am looking for a simple way to split my FASTQ Files into celll …
Dealing with 1000's of demultiplexed fastq files
written 7.3 years ago by volvicpellegrino •

I am trying to process single cell data from a public dataset, containing over 6000 single-end 51 bp fastq files. Each fastq file represents a …
split BAM file with same cell-barcode and UMI pair
written 7.8 years ago by newbinf •

I have single cell RNA-seq reads (from 10x Chromium) that have already been pre-processed. The cell-bacode and UMI tag were cut-and-pasted to the header (via …