Hello,
I tried to plot the distribution of the distance between mates among mapped reads (RNASEQ) and I noticed many different distributions.
I was expecting the distance being centered at x, where x is the insert size.
Is this true for all chromosomes? How about those reads that map mitochondria DNA?
THANKS!
•
written
by
Pfs
1 answer
Let's say your insert size is 100 bp starting out. You get your reads from the RNA, then map those reads back to the reference genome. You'll quickly figure out that many of the paired reads have ends mapping to two different exons. This means that if you take the mapping distance between the reads, your distance is going to be 100 bp (your insert size) PLUS the size of the intron between the exons.
The median intron size is about 1300bp, but is highly variable, so you're not going to get a clean peak around your insert size.
Log in to answer this question.
-
Hisat2 disable spliced alignment or not for DNA shotgun reads??
written by Jon •
-
warning[vg::giraffe]: Cannot cluster reads with a fragment distance smaller than read distance
written by cassiwatt •
-
biological vs. technical definition of "insert size"
written by sovrappensiero •
-
Bowtie: How do I map files to reference genome?
written by melissachua90 •
-
Bowtie not aligning shorter fragments from ATAC-seq
written by eyleen.corrales •
-
Insert size on PE 2*250
written by David
-
Set maximum insert size with BWA MEM?
written by James Ashmore
-
Does it make sense to limit insert size when mapping on transcriptome?
written by Matteo Schiavinato
-
the defination of insert size in RNA-seq, including the length of read, or not?
written by zju.whw •
-
Splice-aware insert size for RNA-seq
written by shuelga •
If your experiment worked, would the distances not merely reflect the distances between exons? If so, that would almost certainly be a multi-modal distribution.
your saying the distances varied between different chromosomes?
you're saying the distances varied between different chromosomes?