Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Novoalign options

    Hello Group,

    I've been trying to find the optimal parameters for aligned bisulfite treated 76bp PE reads to a small reference that includes repeats. I'm confused with the option to set the fragment length and standard deviation. My question is: what exactly does the fragment length refer to? Is it the distance between mapped mates or does it include the reads.

    The library I am using isolated adapter fragments between 250-300bp (minus adapter = 131-181). Given 76 bp reads, some mates could overlap. If the fragment length refers to distance between aligned reads than this value could be negative. Anyone want to clarify?
    Thanks a lot!

  • #2
    HI, Just found your post. The fragment length mean and standard deviation are outer coordinates, it doesn't matter if you set them a bit higher than the real values.
    Colin

    Comment


    • #3
      So by outer coordinates you mean the most distal coordinates of an adapter-ligated molecule, i.e. average length of your library that you would see after PCR amplification on a gel?

      For the library I mentioned, I used 150bp as the mean fragment length since I assumed it was the distance between reads, and the alignment worked pretty well. Is it more computationally strenuous to have a larger fragment length setting?

      Comment


      • #4
        It's the length of the DNA fragment as mapped by the aligner and hence doesn't include any adapters. It's not the length of the gap between the two read alignments.
        If your gel includes PCR adapters then you need to adjust for this.

        It doesn't really affect computation as Novoalign adjusts to the fragment lengths it sees, but setting it too short may mean reads don't get paired properly.
        It's not a major issue as range for pairs is from 0 to mean + 6 standard deviations so it's usually enough.
        You can also run Novoalign on a few K reads and check the reported fragment length distribution. Use the -# option to limit the number of reads processed. e.g. -# 2K will map 2000 reads and stop.

        Comment


        • #5
          Great, that makes a lot of sense. I actually meant "reads" instead of "adapters." Of course in Illumina sequencing the read begins before the adapter since the sequencing primer anneals to the adapter... I'm glad you caught that. I appreciate all the help. Thanks a lot!

          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
          12 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
          68 views
          0 likes
          Last Post seqadmin  
          Working...
          X