Comments (1)
On digging into this issue a bit deeper: some samples in the 5' dataset have 26 bp on the first mate while other samples have 150 bp.
For samples that have 26 bp, the first mate contains just the cell barcode and UMI. I had initially made the assumption that STAR would just map these as single-end reads even when using the above 5' parameters. However, when looking at the SAM flags I realized that all reads in these samples were not properly paired. I mapped these as single-end reads, and the counts improved.
For the samples that have 150 bp, I verified that the first 39 bases contains the cell barcode, UMI, and the 13 bp TSO. When these were aligned using the 5' commands above, STAR had marked all reads as properly paired but the counts were low as indicated above. However, if these samples were aligned as single-end reads, the counts improved and were close to the Cellranger.
The figure below compares the counts when aligned as single-end reads:
Please could you let me know any insight on why the samples samples with 150 bp on the first mare are getting low counts when aligned as paired-end reads?
from star.
Related Issues (20)
- Segmentation fault (core dumped)
- STAR alignments generate different counts: Full genome vs Split genome
- fatal error, cannot create output file for unmapped.out.mate1.thread0
- Possible to implement special treatment of multi-gene reads for BCR/TCR genes?
- [Version 2.7.11b] - Error running 'make' command HOT 6
- killed after "...starting to generate Genome files"
- Problem with reading FASTQ files
- secondary alignment bit in flag not set for unmapped mate of a pair with the other mate mapped
- Chimeric reads very high
- soloUMIfiltering distribution tracking?
- --quantMode GeneCounts output lacks two T-cell receptor genes compared to htseq-count
- STAR hangs at 300% CPU
- Some reads do not seem to be assigned to genes correctly by starsolo
- Alignement of viral genome with reads from infected host
- Fatal INPUT FILE error, no exon lines in the GTF file
- [2.7.10b] Read length in FATAL ERROR message does not match actual read length HOT 2
- STARSolo, alignment of Scell from different lanes
- FATAL ERROR in input read file: the total length of barcode sequence is 151 not equal to expected 28, and also not mapped reads HOT 1
- STARsolo for chimeric reads
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from star.