Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Pinal
    Member
    • Feb 2012
    • 11

    #1

    large genome assembly using paired end and mate paire reads.

    Which assembler should i go for if i want to assemble genome of 1 -1.5 GB size. i have illumina paired end and mate pair reads of 101 bp length. how can i use mate paire reads for Scaffolding?
  • sarvidsson
    Senior Member
    • Jan 2015
    • 137

    #2
    I'd try Allpaths-LG (http://www.broadinstitute.org/softwa...paths-lg/blog/) if your paired end reads mostly overlap (i.e. fragment size of ~180 bp). It will use your mate pairs for scaffolding.

    Comment

    • Pinal
      Member
      • Feb 2012
      • 11

      #3
      Thank you sarvidsson,

      I have access to 96GB memory and 24 core machine. is it sufficiant to work with Allpaths-LG?

      Comment

      • sarvidsson
        Senior Member
        • Jan 2015
        • 137

        #4
        96 GB can be a bit tight, 24 cores should be fine - I'd expect the assembly to run for up to 3 days. If you error correct and normalize the paired end reads prior to assembly (with e.g. BBNorm http://seqanswers.com/forums/showthread.php?t=49763) you typically reduce memory usage for the assembly.

        Comment

        • Brian Bushnell
          Super Moderator
          • Jan 2014
          • 2709

          #5
          We have 100-200Mbp fungal assemblies that run out of memory (with AllPaths-LG) on 128GB nodes, but complete on 256GB nodes. I'm guessing memory may be a serious problem; you probably are going to need more.

          Megahit is fast and seems to have a relatively low memory consumption, and Minia was designed for low memory consumption, so if AllPaths fails you might try those. Or, buy more memory, which will be essential if you plan to routinely assemble large genomes.

          Comment

          • lac302
            Member
            • Dec 2012
            • 64

            #6
            With that amount of memory I'd recommend SGA...Minia is great too but there isn't a scaffolding option.

            Comment

            • fahmida
              Member
              • Aug 2010
              • 54

              #7
              Allpaths-LG is a good option if you have enough RAM and CPUs. Also I wonder whether one of your PE libraries are overlapping i.e. from Allpaths-LG doc "average separation size must be slightly less than twice the read size, such that the reads from a pair will likely overlap".

              Comment

              • Pinal
                Member
                • Feb 2012
                • 11

                #8
                Thank you all.

                I think i have to go for SGA or minia due to lack of memory. is it a good option to use paired end reads for assembly and then go for scaffolding with mate pair data.?
                which tool would be suitable for 101 bp mate pair data for scaffolding?

                Comment

                • fahmida
                  Member
                  • Aug 2010
                  • 54

                  #9
                  Originally posted by Pinal View Post
                  Thank you all.

                  I think i have to go for SGA or minia due to lack of memory. is it a good option to use paired end reads for assembly and then go for scaffolding with mate pair data.?
                  which tool would be suitable for 101 bp mate pair data for scaffolding?
                  You can start with SSPACE.

                  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...
                    Yesterday, 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
                  29 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-28-2026, 09:51 AM
                  0 responses
                  20 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-22-2026, 11:13 AM
                  0 responses
                  71 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-22-2026, 10:48 AM
                  0 responses
                  48 views
                  0 reactions
                  Last Post SEQadmin2  
                  Working...