Question: Plot normalized coverage across genes/transcripts in single cell rna seq data
1
gravatar for kevin.stachelek
10 days 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 • 67 views
ADD COMMENTlink modified 10 days ago by kristoffer.vittingseerup2.9k • written 10 days ago by kevin.stachelek10
3
gravatar for kristoffer.vittingseerup
10 days ago by
European Union
kristoffer.vittingseerup2.9k 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 10 days ago by kristoffer.vittingseerup2.9k

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

ADD REPLYlink written 5 days ago by kevin.stachelek10
Please log in to add an answer.

Help
Access

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