This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to explain expression level peaks(uneven gene expression) within a gene?

Hi all, I have a question about gene expression. I mapped RNA-seq reads back to the genome and noticed there are my expression peaks within a gene. Does it cause by library preparation process, sequencing problem, or a real biological thing? The blue bar on the top is a gene(single-exon gene), the red below is bam file generated by mapping reads. Has anyone got the same problem?

I uploaded the figure of a gene expression pattern on my GitHub, here is the link. Hope it works. I was thinking if there is a way to make uploading figure easier?

https://github.com/YIRAN117/Antisense/blob/master/Q%20about%20gene%20expression.png

rna-seq assembly expression

What is the length of the exon ? Though it should not matter. Do you have a screen shot from IGV ? Did you filter for uniquely mapped reads ?

I think maybe I understand what is happening with the uneven expression. FastQC showed enriched clusters at the beginning of reads which means the library is biased by random priming during library preparation. Maybe it is why...

0 answers

No answers yet.

Log in to answer this question.