Question: Downstream analysis of bam files from scRNA-Seq Data
gravatar for rajeev.vikram
2.4 years ago by
rajeev.vikram30 wrote:


I am new to single cell RNA-Seq, and am trying to find a workflow related question. I have aligned about 1500 cells from two different experiments (but using same experimental procedure). I am wondering how do I go ahead with differential expression analysis. Should I merge all the bam files and then generate a count table? or should I use string tie on each cell (use the hisat-stringtie pipeline). I am a bit confused because the data set is too large. I understand that ultimately I should have an expression matrix (Genes/Expression level/cell). What should be the next step after alignment (i,e creation of individual bam files) ? Please advice.


ADD COMMENTlink modified 2.4 years ago by Devon Ryan96k • written 2.4 years ago by rajeev.vikram30
gravatar for Devon Ryan
2.4 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

Normally one would use featureCounts or something to get the per-gene/per-cell counts. In single-cell studies this tends to involve using UMIs for deduplication, so be sure to account for that if possible. If you happen to know what groups the various cells belong to then you can go ahead with differential expression then. Otherwise you'll need to come up with your own groups (via tSNE or what not) and then do the DE between those groups.

ADD COMMENTlink written 2.4 years ago by Devon Ryan96k

Thanks Ryan, This makes things a lot more clear.

ADD REPLYlink modified 2.4 years ago • written 2.4 years ago by rajeev.vikram30
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1164 users visited in the last hour