I got a strange distribution of insert size as analyzed by Qualimap on the BAM file. Not totally aware of the library preparation procedure yet, but this is a PE-sequencing on NextSeq500, libraries were made with NEBNext Ultra II kit and size selection beads. The median insert size was expected to be around 300bp. The problem is that we got not only a much smaller median insert size (190bp), but also a non-symmetrical distribution. What could be the reason? Something related to the trimming I used to remove the adapters? Thanks for any help.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
-
Can you tell us what aligner you used for this? Is it possible that the aligner kept reads with reasonable distance (did you provide an expected fragment size) when mapping and discarded smaller inserts (what was the alignment %?). Otherwise you are seeing a distribution of insert sizes that were in your dataset.
-
Are you sure the 300bp was the insert size and not the total fragment size? The adapters add about 150bp to the total fragment size so this looks about right in that case. The asymmetric distribution on the left hand side is caused by the size selection beads removing the shorter fragments.
Comment
-
From the fastQC results, we saw plenty of adapter sequence, which could be attributed to the wrong (smaller) insert size selection. So I cut the adapter sequence and did the alignment. The shown weird distribution was from the alignment after trimming. Now I have just finished alignment with the raw fastq files without any trimming, the distribution of insert size looks more symmetrical, at least no sudden loss of small inserts. So it seems the shape was caused by trimming (and loss of reads without minimum length which I set at 80). Anyway, since our PE read length is 151, the library was simply not prepared with the right insert size (300). Thanks for all the comments. I appreciate it!Last edited by zacshi; 05-08-2018, 08:28 PM.
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...-
Channel: Articles
Yesterday, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
51 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
45 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment