Queries on 10x sequencing 3' Read 1
I have seen on umi_tools that the read 1 only contains the cell barcode and UMI. However, when I look on my data, the read length is much longer than 28. Does that still mean that my read 1 contains only the cell barcode and UMI?
10x
• 1,623 views
•
link
updated
by
ATpoint
•
written
by
Cheng Wei •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to Analyze 10x Genomics Data with Full-Length R1 and R2 (100bp Each)?
written by zbidav •Dear Biostars community, I’m used to analyzing 10x Genomics single-cell RNA-seq data where R1 is typically 26bp (cell barcode + UMI) and R2 is longer …
-
Extract the true single-cell RNA sequencing reads for running SAHMI
written by 573704669 •Hello everyone! I’m currently using the SAHMI pipeline to annotate microbiome information from the single-cell data. However, I encountered two potential problems when applying SAHMI …
-
Can I use umi_tools to discard cell barcodes and extract the UMI while keeping the UMI in the read?
written by cag104 •I have paired end reads where the read structure is as follows: ``` >Read1 5'-6N-20CB-XXXXXXXXXX....XXXXGS-3' ``` The 6N corresponds to the UMI, the 20CB corresponds …
-
Improving my STARsolo results
written by Aaron •I'm pretty new to using STARsolo but I've built a STAR index, and I've just tried to align and count. I'm not quite sure what …
-
UMI extraction from 10X visium spatial transcriptome data
written by Omics data miningHello everyone I have to analyse visium spatial transcriptome (ST) sequencing data (2 x150 bp) . I want to extract Spatial barcode and UMI from …
-
Will mapping Read 1 in Single Cell 3' sequencing increase mapping accuracy or not?
written by Zhilong JiaWhy lots of pipeline of 3` sequencing (such as [zUMI][1], [scruff][2]) do not map Read 1, including barcode and UMI, as PE mode mapping? Though …
-
umi_tools using regex greedy quantifiers only returning 1 character
written by williamtmills •I am trying to use umi_tools to remove UMIs and cell barcodes and leave the remaining sequence. Unfortunately, after correctly removing the umi and cell …
-
STARsolo read1 length greater than cell barcode length and UMI length
written by from the mountainsI am trying out STARsolo (STAR 2.7.3a) for the first time and wondering if i can side-step an additional trimming step. my reads are sequenced …
-
umi_tools extract error: IndexError: string index out of range
written by Rituriya •Hi All, Background: I have completed adapter trimming and checked QC on Illumina NextSeq **miRNA** single end reads of length 75bp. I want to run …
-
How to append the cell barcode and UMI information to the fastq header in paired-end single-cell RN…
written by sc243 •Hello everyone, I have a question regarding my single-cell RNA-seq data. I have the following pair-end data in `fastq.gz` format. Read1 (contains 6bp UMI, followed …
Yes, the rest might be adapters or garbage, you can blast it if you like.
Here is an absolutely excellent resource on single cell lib structures.
https://github.com/Teichlab/scg_lib_structs
The rest of R1 would be the oligo dT primer and eventually the 3' sequence of the transcript. Illumina sequencing would have a hard time with base calling on the polyA region though due to the low complexity problem.
That having said, it is on the user to define on the sequencer how long a read is. If you sequence the library as a 2x150bp run (very common run mode on Novaseq machines, most cost effective) then R1 is 150bp, but the meaningful part is much shorter as said above. One would simply discard (many tools do that automatically) everything after the CB/UMI part as it is basically garbage.