Comments (4)
my experience is that mapping in that region can be quite difficult and that is one of the reasons IgBlast exists.
One of the major issues is that the parameters you'd need to use for mapping to that region are different from the preferred ones for mapping to the rest of the genome which contains much less genetic variation.
In any case, I would advise consulting how it is is done in the many published workflows for assembling BCR and TCR from sequencing reads see BALDR, BASIC, BRACER etc.
from star.
Ok thank you, I wanted to see if it was relevant to look at the coverage from the STARsolo BAM output on IGV to better see from where the output quality (incomplete reconstruction of certain chains) of some of these TCR/BCR reconstruction methods on a sample stem.
from star.
that makes sense as a line of investigation -- if you're looking at why some assemblies fail I wonder if it could also be informative to look at how many reads map to the constant regions of the chains using the BAMs from STAR. I've found that is a more reliable quantification than the variable region
from star.
Dear @michael-swift, thank you for your advice. I am using a 5' biased approach which makes it so that the constant region is less covered. I will try to see if IgBlast gets me a really different BAM profile than STARsolo at the TCR / BCR loci. Thank you
from star.
Related Issues (20)
- 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
- Parameters for 10x 5' alignment HOT 1
- 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
- STARsolo UMI calling problem (UMI position shift) after cell barcode correction by barcode whitelist
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.