Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • lh3
    Senior Member
    • Feb 2008
    • 686

    #16
    Different assemblers respond to errors and trimming in different ways. For example, if you read the SGA paper, the authors recommend not to trim reads. Allpath-lg does not trim reads, either, as I remember. Other assemblers, such as SPAdes, may be less sensitive to trimming as they trim reads by default. Also, I have read somewhere (could be wrong) that SOAPdenovo developers recommend not to correct reads if you have enough RAM, but SGA/Allpath-lg etc always include error correction as a necessary step. At the end of day, which trimming/error correct approach to use is assembler dependent. If it were me, I would just use the tools/pipelines recommended by the developers. If I had time, I would combine different strategies/correctors and see what I would get. Probably the result is data dependent.

    K-mer based error correctors typically use short k-mers. I think that is fine. With shorter k-mers, we more often collapse segmental duplications/repeats and will not be able to correct errors when they occur right at the sites differentiating repeats. However, only a small fraction of errors are not correctable due to repeats. If such errors can be corrected with long k-mers, assemblers can usually handle them well. I would not worry to much about the k-mer length in error correction, unless it is too short.

    Comment

    • jwag
      Member
      • Apr 2013
      • 42

      #17
      Originally posted by lh3 View Post
      Different assemblers respond to errors and trimming in different ways. For example, if you read the SGA paper, the authors recommend not to trim reads. Allpath-lg does not trim reads, either, as I remember. Other assemblers, such as SPAdes, may be less sensitive to trimming as they trim reads by default. Also, I have read somewhere (could be wrong) that SOAPdenovo developers recommend not to correct reads if you have enough RAM, but SGA/Allpath-lg etc always include error correction as a necessary step. At the end of day, which trimming/error correct approach to use is assembler dependent. If it were me, I would just use the tools/pipelines recommended by the developers. If I had time, I would combine different strategies/correctors and see what I would get. Probably the result is data dependent.
      I've been trying to follow the pipeline of the folks that assembled the Giant Panda, as they used SOAPdenovo and exclusively short Illumina reads (I think the group that developed SOAPdenovo is the same that put out the Giant Panda assembly). Unfortunately, documentation with SOAPdenovo isn't the most thorough. According to their Nature paper (supplemental section), they did trim for low quality bases at the 3' end before error correction, though they don't specify their threshold.

      After looking at the parameters in SOAPec, it seems that they might have trimmed for Q<2 because that is the only threshold available for trimming during error correction in that program.

      I'm going to first try letting SOAPec do the trimming and error correction from the full raw data. Right now I have all paired-end data, and I don't think SOAPec can handle both paired-end and single end data while still keeping pairs together in the output.

      If that doesn't work too well, I'll try trimming first with something like Q<10 before error correction.

      Comment

      • Wallysb01
        Senior Member
        • Feb 2011
        • 286

        #18
        Originally posted by jwag View Post
        I'm going to first try letting SOAPec do the trimming and error correction from the full raw data. Right now I have all paired-end data, and I don't think SOAPec can handle both paired-end and single end data while still keeping pairs together in the output.
        You just have to do two seperate runs with the same output from KmerFreq.

        Comment

        • jwag
          Member
          • Apr 2013
          • 42

          #19
          Originally posted by Wallysb01 View Post
          You just have to do two seperate runs with the same output from KmerFreq.
          Ah that makes sense. So I can just initially put in all my data into KmerFreq, so that it counts all instances of K, and I can use any of my data against that same frequency distribution (even in smaller chunks). That will definitely save me some time. Thanks for the tip.

          Comment

          Latest Articles

          Collapse

          • SEQadmin2
            16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
            by SEQadmin2




            NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
            Today, 06:46 AM
          • SEQadmin2
            Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
            by SEQadmin2



            Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

            Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
            ...
            10-14-2026, 09:55 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 10-28-2026, 09:57 AM
          0 responses
          22 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-28-2026, 09:51 AM
          0 responses
          14 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-22-2026, 11:13 AM
          0 responses
          67 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-22-2026, 10:48 AM
          0 responses
          47 views
          0 reactions
          Last Post SEQadmin2  
          Working...