Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Illumina fragment length distribution

    I'm trying to characterize the fragment length distribution of some Illumina data from the 1000 Genomes Project.

    To do this, I downloaded the supplied BAM files (reads aligned with MAQ) for the mtDNA of the CEU trio daughter (to make sure rearrangements wouldn't be a problem.) Then I checked the distance between the two farthest ends of the mapped pairs.

    I thought this would be fairly straightforward, but the mean and median aren't even close to the fragment size listed for some of the libraries. They're much lower than the listed fragment size. The distributions aren't even remotely normal, either.

    For example, SRR001139 is supposed to have a fragment length of 676 with mate paired reads of size 47. The mean fragment length I estimate, however, is 234 with a standard deviation of 148.

    Here is a histogram of the length distribution:


    Is expected? Or am I making a mistake somewhere along the way in my reasoning? Does anyone know any references about the fragment lengths that result from Solexa sequencing?


    Any help or suggestions would be much appreciated.

  • #2
    Sorry. I said "fragment size" here somewhat carelessly. I really mean "insert size".

    Comment


    • #3
      randomly came across this - kinda sad that you've had no answers - doesn't bode well for the community.

      In any case, could this have been coming from contaminated "mate-pair" libraries (from large circularized fragments that are then fragmented themselves - fragments that contain the circ. point lead to outward facing read pairs that map R/F at long distances, while fragments that don't contain the circ. point are F/R pairs that map at short distances. I.e. all "mate-pair" libraries are contaminated to a certain extent with shorter-range "paired-end" reads. Now, given the opposite orientations, normally one or the other set won't map, but I think MAQ was able to align both, IIRC.

      HTH.
      ~J

      Comment


      • #4
        Does Illumina library's fragment length of 676 include adapter sequences? How big were the adaptor sequences?

        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
        11 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        51 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