Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Weird distribution of insert size

    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.



  • #2
    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.

    Comment


    • #3
      Thanks! The aligner used was bwa-mem, nothing set up to restrict the fragment size. Mapping rate was like 99.8%. I guess that I should re-do the alignment with the fastq files before trimming just for curiosity.

      Comment


      • #4
        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


        • #5
          It must be the total frag size. Since the above plot is post alignment it will correspond nicely to 150 bp inserts.

          Comment


          • #6
            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

            • seqadmin
              Strategies for Sequencing Challenging Samples
              by seqadmin


              Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
              03-22-2024, 06:39 AM
            • seqadmin
              Techniques and Challenges in Conservation Genomics
              by seqadmin



              The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

              Avian Conservation
              Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
              03-08-2024, 10:41 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 06:37 PM
            0 responses
            10 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            9 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            49 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            67 views
            0 likes
            Last Post seqadmin  
            Working...
            X