Question: Plot normalized coverage across genes/transcripts in single cell rna seq data
gravatar for kevin.stachelek
10 months ago by
kevin.stachelek10 wrote:

I'd like to generate a plot similar to the included image. What is the name of such a plot and does there exist a simple tool to generate, preferably compatible with multiqc and/or R/bioconductor?

enter image description here

rna-seq single cell R • 447 views
ADD COMMENTlink modified 9 months ago by kristoffer.vittingseerup3.4k • written 10 months ago by kevin.stachelek10
gravatar for kristoffer.vittingseerup
9 months ago by
European Union
kristoffer.vittingseerup3.4k wrote:

Maybe CoverageView? Looks like it was made for other experiment types but since it uses bam files as input it might work. Else Bioconductors bamsignals have a bamCoverage() function

Else I know RSeQC have a gene-coverage tool which you can find here.

Just be aware that generally in single cell data (especially for 5' end data such as 10X) you have way to few reads per gene to look at coverage. Most genes have very few reads so how much sense it makes to calculate this I'm note sure.

ADD COMMENTlink written 9 months ago by kristoffer.vittingseerup3.4k

thanks. I'm working with Smart-seq data so I hope I have enough reads to make an intelligible plot.

ADD REPLYlink written 9 months ago by kevin.stachelek10
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: 1444 users visited in the last hour