Bimodal insert size distribution

Pepe · 12-17-2009, 07:30 PM

Hi,

we have used the RNAseq Illumina protocol to make about a dozen paired-end libraries.

Once sequenced and aligned to the reference genome I plot the distribution of insert sizes and I see a bimodal distribution. One peak corresponding to the insert size I'm expecting from the band size we cut (around 220 bp) and the other peak at a smaller size (around 100bp). See attached.

I see different sizes of that peak at 100bp in some of the libraries and not at all in others, which makes me think it happened during the library preparation.

I look at the reads under both peaks and they look very similar to me (many perfect matches, evenly distributed across the genome...)

I was expecting bigger insert sized due to the alignment of the (RNA) reads to the genomic sequence (no introns), but not smaller.

Any ideas on what's going on?

Thanks in advance

natstreet · 03-03-2010, 05:10 AM

I realise this is a long time past the original post - but could you tell me how you produced the graph you attached?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Insert size != Fragment size?	Boel	Bioinformatics	6	12-12-2013 09:28 AM
About Insert, Insert size and MIRA mates.file	aarthi.talla	454 Pyrosequencing	1	08-01-2011 02:37 PM
Bimodal Distribution of Mismatch Frequencies in Illumina	mdjones66	Bioinformatics	0	10-27-2010 08:41 PM
insert size	polystone	Sample Prep / Library Generation	0	05-04-2010 11:07 AM
insert size	adrian	Bioinformatics	1	03-18-2010 05:55 PM

03-03-2010, 05:10 AM	#2
natstreet Member Location: Sweden Join Date: Nov 2009 Posts: 83	I realise this is a long time past the original post - but could you tell me how you produced the graph you attached?