I have complete sequences from 10x Chromium that look like the following:
|----my 16 bp cell barcode----|---8bp 10x UMI---|--SO--|-------------------cDNA (100 bp)----------------------| SO = switch oligo
The reads are demultiplexed and don't have the Illumina barcodes in the sequence. I want to map the cDNA (sequence of interest), however I want to preserve the data provided by my barcode, 10x UMI, and switch oligo (about 40 bp). This is so I can identify duplicates and which cell each read originated from.
I can found trimming softwares (cutadapt and trimmomatic) that can delete the the 40bp. However, because I want to retain this information, I do not want to simply trim the sequences. On the other hand, keeping the sequences have lead to 70-90% not mapping because of the 40bp of non-genetic information.
I have thought about 'cutting and pasting' the information into the header, but am not aware of a program that already does this.
How to I align these reads while preserving the data provided by the 40bp barcodes?
1 answer
Have you looked at umitools (https://github.com/CGATOxford/UMI-tools )?
Log in to answer this question.
-
How to get host and virus UMI count per cell with Bacterial genome data?
written by San •
-
With 10x Visum HD : How to map the barcode id to the barcode sequence ?
written by sacha
-
Using UMI-tools on Smart-seq3 RNA-Seq data
written by Agastya •
-
ScRNA data
written by abbas.waseem.gcu •
-
10x 3' library creates R1 and R2 fastq files with the same read length
written by tomas4482 •
-
Alignment to scRNA-seq data
written by tomleee012 •
-
Tutorial: STARsolo config for 10x Chromium v1, v2, v3
written by Kevin Blighe
-
A question about drop-seq and Nextera XT Library prepare
written by 610225668 •
-
split BAM file with same cell-barcode and UMI pair
written by newbinf •
-
Tool to append Illumina Read 1 to Read 2 for downstream demultiplexing
written by achamess
