Metagene for first last and intermediate exons?

Hi all,

I was able to do the metagene plot for 5'utr exon and 3' utr. My question is, is there any way to do the metagene plot for first, last and intermediate exons.

ggplot rna-seq deeptools r

• 1,428 views • link

updated 3.1 years ago by Ram

• written 3.6 years ago by PK •

You can make a metagene plot for any shared feature among genes. Where are you stuck? What have you tried?

• 0 views • link

written 3.6 years ago by seidel

Hi, I used following command to produce the metageneplot

cvgs.utr5 <- coverageDepth("bamfilename", gtf = txdb, region="utr5")

and i did for CDS and utr3

metaPlot(cvgs.utr5, cvgs.CDS, cvgs.utr3, sample=1,xaxis = "mRNA")

which produced this plot enter image description here

But if you look at towards stop there is clear spike which is present in all the KO compare to wild type. Additionally, i did for each files is there any way to combine replicates together?

• 0 views • link

written 3.6 years ago by PK •

0 answers

No answers yet.

More posts like this

Sequence Homology locator
written 3.5 years ago by himalayan_yeti •

Hello helpful folks, I have several human genes of interest (with their gene name, gene id, sequences, and accession number information) that are homologous to …
How to normalize the counts when you are looking at 3'UTR?
written 3.5 years ago by PK •

I want to know whether I have elevation at 3'UTR between two conditions. What I did was took all the transcripts and extracted the coverage …
how to create a metagene plot for different features of genes?
written 3.6 years ago by PK •

Hi all, I wanted to do metagene plot for different features of genes. I have RNA-Seq sample. So would like to create a plot which …
Exons Distribution from single cell data.
written 5.1 years ago by Cheng Wei •

I was looking at the read counts from my single cell sequencing data through IGV and notice that there are very little reads that mapped …
genomic features 3'utr and 5'utr
written 6.1 years ago by ginlucks •

Hello, I am having some troubles using genomicFeatures R's package. Basically I am trying to retrieve for a specific gene the longest isoform and, for …
Bulk download introns, exons, and UTR regions from Ensembl for gene prediction training set
written 8.6 years ago by Katherine Huang •

Hi, I would like to download labeled FASTA sequences of introns, exons, 5' UTR regions, and 3' UTR regions from a nonredundant set of human …
Alternative for GeneBody coverage from RseQC (Python) in R
written 9.5 years ago by Walter Baumann •

Hi, There is a function called geneBody_coverage.py in the Python tool RseQC, which determines whether the coverage of Reads is uniform along the transcripts, o …
RNA-seq: bias in coverage
written 9.5 years ago by angela.sparago •

When I observe mapped reads of RNA-seq data, I realize that many genes have the first exon not covered. I just read this seems to …
Collecting 3' UTR and neighboring exonic sequence from GENCODE GTF files
written 10.3 years ago by mdwain.uw •

Hi, I performed RNA sequencing using a poly A 3' tagging/sequencing approach. I therefore expect the sensible reads to map to only the 3' end …
Question about how fetch intron location
written 11.6 years ago by jesselee516

Merry Christmas and Happy New Year, Hi,I am new to this research.I am trying to understand how to fetch intron location.I have get gene location …