Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Insert size != Fragment size?

    Hi All,

    This is a very simple question, or should be, but there seems to be some confusion out there.

    My understanding is that insert size is the number of bases between paired end reads (example: 300 bp fragments, 2*75 reads -> insert size 150). However, looking trough different threads here and information elsewhere the fragment is sometime referred to as insert as well, making insert size == fragment size.

    I wonder whether different programs (BWA, picards CollectInsertSizeMetrics, bfast, samtools etc) have different definitions of insert size, which might make things messy. Right now I am especially interested in picards definition.

    Any ideas?

    Thanks,
    Boel

  • #2
    I'm not sure about Picard specifically, but I found a thread that discusses insert size here, which seems to suggest that the insert size is the stretch of sequence between the adapters (so in your example, 300 would be correct). However, certain tuxedo suite programs (tophat/cufflinks/bowtie) take a --mean-inner-dist option, defined as fragment length - reads.

    Comment


    • #3
      Yes, I've seen both definitions used by different software. I know Bowtie uses insert size == fragment size for example.

      Personally I like this definition because you could sequence the same library with different read lengths, or you could have variable length reads (Ion Torrent) or you could trim your 3' read tips. Each step would vary the insert size, but the fragment size would remain constant.

      I think part of the difficulty comes from the difference between Illumina paired-end protocols (e.g. bidirectional sequencing) where insert size is always related to fragment size and the long mate-pair/jumping protocols, where the insert size relates instead to the sizing step (e.g. 8kb gel-cut) and is independent of fragment length.

      Comment


      • #4
        According to samtools help (which includes picard):
        "For Illumina paired-end data, the inferred insert size would be the difference between the 5' positions of the two reads." This translates to 300 bases in my previous example, since nucleotides are added in 5' to 3' direction.

        Thanks for replying!

        Comment


        • #5
          Yup - Glad you found an answer!

          Comment


          • #6
            Whether insert size = fragment size does vary from tool to tool, I would specifically look it up for whatever you're using.

            Regardless, the total adaptor-insert-adaptor length is useful as the best predictor of a library's amplification behavior (both in qPCR QC and on the flowcell). For example, you wouldn't want to pool together two libraries with identical 200bp inserts but very different adaptor + index lengths (unless it was acceptable for the majority of the reads to come from the smaller construct).

            Comment


            • #7
              Don't bother, use a different term. You will be misunderstood. When you mean fragment size, say "fragment size", when mean the distance between the pairs. Say, "distance between the pairs".

              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
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              8 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