Detecting duplicated UMI in BAM file
Dear all,
I'm trying to remove reads in BAM file resulted by Cell Ranger that have the same UMI per CB and keep only one, since I need to filter by some condition and count uniquely satisfied UMI. I'm using pysam and function is_duplicate() to detect those reads. I'm wondering if it is enough to detect duplicated reads? And if sam flag is updated after correcting UMI and cell barcode in cell ranger?
And also, if someone know how sam flag recognizes duplicated reads?
Thanks for your help.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Spliced and unspliced counts from HIVE scRNA-seq data
written by anna_o •
-
Editing/ Adding new tags to BAM file to match 10X format
written by ashwini •
-
CellRanger v. 2.0 and v. 7.0
written by Tadeoye •
-
How to find UMI counts from cell ranger VDJ output
written by Kevin •
-
How Cell ranger count UMI from intronic region in the latest version
written by tien •
-
Read is confidently mapped to the transcriptome in Cell Ranger
written by tien •
-
Count number of unique cell barcode (CB) in bam file
written by tien •
-
Cellranger output .bam file does Not contain barcodes ( CB and UB )
written by Gabriel
-
Editing BAM file using pysam
written by sam •
-
Demultiplexing Conceptual Issue
written by asyndeton17 •
You'll need to check whether cell ranger adds the duplicate flag for that to work. If not you would need to filter by the cell barcode and UMI tags.
Hi, thanks for your respond. Cell ranger has such flag (at least, there exists reads with samflag = 1024). But I'm not sure if it's enough? Because I observe something like in this post UMI and CB for the same sample but different lane by 10X single-cell RNA sequencing