Question: How to get an estimate the reads mapping to different regions of the genome and plot to get the peak
0
gravatar for Biocode_user
3.3 years ago by
Biocode_user10
Sweden
Biocode_user10 wrote:

How can I get a density plot to get an idea of reads mapping to different regions of the genome.

I see that only 40% of the reads maps to CDS and around 60% mapping to UTR's or some other non-coding RNA. I have generated a bam file after mapping the reads to genome. I could like to visualize this by generating a plot kind. Since I am new-bie to this, it would be helpful if someone could help me out with what tools could be used to get this estimate and how to plot them.

I see from other posts that picard tools can be used. Also wiggle format can be used to solve this problem. Is there any reccomendation or how to get this done?

 

Thanks,

Bio

density plot rna-seq assembly • 1.2k views
ADD COMMENTlink modified 3.3 years ago by Sukhdeep Singh9.7k • written 3.3 years ago by Biocode_user10
0
gravatar for 5heikki
3.3 years ago by
5heikki8.4k
Finland
5heikki8.4k wrote:

ggplot2 is IMO the best tool for plotting. Here's a nice resource.. From the same source you can find out e.g. how to read tables into R. I very much recommend that you install RStudio. Also, don't expect to master it in minutes. This will take hours/days/weeks depending on your level of expertise. Good news is that Google will know the answer to pretty much any question you might have..

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by 5heikki8.4k
0
gravatar for Sukhdeep Singh
3.3 years ago by
Sukhdeep Singh9.7k
Netherlands
Sukhdeep Singh9.7k wrote:

If you are comfortable with R, there are several packages which will make this task easy.

ChIPseeker: an R package for ChIP peak Annotation, Comparison and Visualization

ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip data

 


 

ADD COMMENTlink written 3.3 years ago by Sukhdeep Singh9.7k

I have not worked much with R before. And its not just about plotting  . But I would like to know from scratch what should be done with the bam alignment file to get how the reads are mapped to different regions. I have no idea of how to get the read counts at different regions?

ADD REPLYlink written 3.3 years ago by Biocode_user10

You can convert Bam file (which contains information about mapped reads) to Bed file (which retains read location information per read per line) using bamToBed and then can use one of the above tools. Otherwise, have look at bedTools which can be used to intersect bam files with genomic files (ex UTR, promoters etc) to give what you are essentially asking for.

ADD REPLYlink written 3.3 years ago by Sukhdeep Singh9.7k
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: 792 users visited in the last hour