Hello, I am visualizing RNA-seq PE reads using IGV and I see some strange but consistent patterns in the variation of depth of expression within some transcripts. I am aligning reads to a transcriptome that was assembled de novo. Is there a simple explanation for the pattern depicted in the attached image (see post directly below for image)? (There are two tracks, each representing a different individual.)
Is this alternative splicing? reads from two different paralogs aligning to the same transcript? an error in the construction of the original transcript in the denovo assembly that caused two different transcripts to merge?? If anyone out there has a lot of experience looking at these kind of alignments in RNA-seq data maybe it is obvious to you?
Many thanks!
Is this alternative splicing? reads from two different paralogs aligning to the same transcript? an error in the construction of the original transcript in the denovo assembly that caused two different transcripts to merge?? If anyone out there has a lot of experience looking at these kind of alignments in RNA-seq data maybe it is obvious to you?
Many thanks!
Comment