RNAseq, reads outside exons

Hello. Why do readings outside of genes appear in an RNAseq experiment in IGV? The readings are supposed to match exons, right?
Thank you

1 Like

It is possible to have reads align to regions of the genome not specifically targeted by the library construction/sequencing method.

BAM datasets can filtered to reduce spurious hits, for example, by mapping quality, concordantly aligned paired reads. And can be ignored in some protocols as they are excluded by the analysis tools in various ways.

Some Galaxy tutorials that might be helpful if you are just getting started: https://galaxyproject.github.io/training-material/topics/introduction/

With all here, including those that cover RNA-seq:

Yes, but what is the cause of these alignments in the intergenic regions.
Thank you.

At least two possibilities:

  1. These are real RNA molecules that were in your library. poly-a selection in RNA-seq is not 100% efficient, so these could be ncRNAs, or just spurious transcription

  2. Alignment error: alignment is also not perfect, some reads may be placed in intergenic regions incorrectly.

This is typically not something to worry about. Downstream, if you assemble transcripts or quantify using an existing transcriptome, these will be ignored.

1 Like

Ok Jxtx, I understand.
Thank you very much for your answer.