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

  • #76
    Hi Brian,

    Thanks for your reply. I tried using Spades and got similar results. After speaking with my experimental colleague, I think I may have a mixture of viruses rather than a single one. I ran the program on the input fastq file I used for tadpole. in=ecco.fq hist=histogram.txt

    There is a single, strong peak of singleton kmers. Maybe I'm a little slow on the uptake, but if the majority of the kmers are seen only once, this would indicate it would be hard to produce longer contigs, correct?


    • #77
      Originally posted by indapa View Post
      There is a single, strong peak of singleton kmers. Maybe I'm a little slow on the uptake, but if the majority of the kmers are seen only once, this would indicate it would be hard to produce longer contigs, correct?
      That's correct. However, most datasets have a strong peak for singleton kmers, due to sequencing error. It's just that typically, for isolates, there is also an obvious higher peak (at, say, 40x when you have 40-fold genomic coverage).

      In your case, you could have a mix of viruses with sequencing depth sufficiently low that none of them will assemble, or a single rapidly-mutating virus, or very low-quality data due to a problem with the sequencing machine... though I'd still suggest that host genome contamination is a possibility.


      • #78
        Hi Brian,

        Does Tadpole have a checkpoint option, or this a possible feature to add? I ask because I am de novo assembling the axolotl genome using the available raw paired end shotgun sequencing libraries (mate paired libraries have not been released to the public yet nor the shotgun and mate pair assembly until publication). My Tadpole assembly does not finish within 72 hours, which is the maximum walltime for the bigmem queue I am using (48 cores and 1.5TB RAM).

        The goal of the assembly is not to have the largest N50 possible, rather to map candidate baits (for a sequence capture experiment) developed from transcriptome transcripts and eliminate candidate baits that are non-specific and hybridize to multiple targets.

        I first performed adapter and quality trimming on each set of the 15 sets of raw paired end reads with BBDuk 37.56
        ~/bin/bbmap-37.56/ in1=SRR2027504_1.fastq.gz in2=SRR2027504_2.fastq.gz out1=SRR2027504_1_clean.fastq.gz out2=SRR2027504_2_clean.fastq.gz ref=~/bin/bbmap-37.56/resources/truseq.fa.gz ktrim=r k=23 mink=11 hdist=1 tpe tbo qtrim=rl trimq=15 threads=8
        Then used Tapole 37.56 to de novo assemble them
        ~/bin/bbmap-37.56/ -Xmx1400g threads=48 prealloc=t \
        in1=SRR2027504_1_clean.fastq.gz,SRR2027505_1_clean.fastq.gz,SRR2027506_1_clean.fastq.gz,SRR2027507_1_clean.fastq.gz,SRR2027508_1_clean.fastq.gz,SRR2027509_1_clean.fastq.gz,SRR2027510_1_clean.fastq.gz,SRR2027511_1_clean.fastq.gz,SRR2027512_1_clean.fastq.gz,SRR2027513_1_clean.fastq.gz,SRR2027514_1_clean.fastq.gz,SRR2027515_1_clean.fastq.gz,SRR2027516_1_clean.fastq.gz,SRR2027517_1_clean.fastq.gz,SRR2027518_1_clean.fastq.gz \
        in2=SRR2027504_2_clean.fastq.gz,SRR2027505_2_clean.fastq.gz,SRR2027506_2_clean.fastq.gz,SRR2027507_2_clean.fastq.gz,SRR2027508_2_clean.fastq.gz,SRR2027509_2_clean.fastq.gz,SRR2027510_2_clean.fastq.gz,SRR2027511_2_clean.fastq.gz,SRR2027512_2_clean.fastq.gz,SRR2027513_2_clean.fastq.gz,SRR2027514_2_clean.fastq.gz,SRR2027515_2_clean.fastq.gz,SRR2027516_2_clean.fastq.gz,SRR2027517_2_clean.fastq.gz,SRR2027518_2_clean.fastq.gz \
        out=axolotl-contigs-k63.fasta mode=contig k=63
        If a checkpoint option is not possible with Tadpole, can you recommend a de novo assembler that supports checkpoint?

        Thank you,


        • #79
          The axolotl genome consists of 14 chromosome pairs (2N = 28)14 and estimates of its physical size range from 21–48 gigabases
          I am not sure if Tadpole is designed to assemble a genome of this size. Curious to see what @Brian has to say.

          Have you tried other large genome assemblers? ALLPATHS-LG?


          • #80
            Originally posted by GenoMax View Post
            Have you tried other large genome assemblers? ALLPATHS-LG?
            No I haven't, but based on the manual- I am out of luck as ALLPATHS-LG requires shotgun sequencing and mate pair libraries.


            • #81
              Hi Gopo,

              I don't particularly recommend Tadpole for diploid (or higher) genomes, as it has absolutely no capability of dealing with heterozygous sites. However, it's really fast, so even with a huge genome 72 hours would be unusual (though possible; that one is pretty large after all) unless something went wrong. Are you certain that it not crash? Typically, if it crashed (due to running out of memory, for example) it would indicate that in the stderr output.

              You may find it helpful to perform error-correction with K=31 and add the flag "prefilter=2" to get rid of erroneous kmers and conserve memory with a Bloom filter. But as for finishing a massive assembly in 72 hours, I don't think that will help. Tadpole does not support checkpointing. I don't know what the best diploid eukaryotic assembler for Illumina reads is currently, but it's safe to bet that it's not Tadpole (unless all you care about is avoiding misassemblies and very low continuity is acceptable). There are some assemblers, though, like Ray and Hipmer, that can run distributed on a cluster to reduce the overall time as well as per-node memory requirements. Those might be worth trying in this case to fit into the 72-hour window.

              If your read pairs are mostly overlapping, you can also merge them first with BBMerge to reduce your data volume somewhat and increase quality, which will reduce both time and memory usage. Ray, for example, appears to benefit from merged reads, and I've been told by one of the developers that HipMer does as well.
              Last edited by Brian Bushnell; 10-13-2017, 10:17 AM.


              • #82
                Axolotl genome paper has used SOAPdenovo2.


                • #83
                  Originally posted by Brian Bushnell View Post
                  Hi Gopo,
                  Are you certain that it not crash? Typically, if it crashed (due to running out of memory, for example) it would indicate that in the stderr output.
                  Hi Brian, no it did not crash- unfortunately the job exceeded the allowed walltime. I'll try what you suggested for Tadpole first.

                  @GenoMax - Thank you. Unfortunately, they were unable to finish the assembly with SOAPDenovo2 (see

                  From ""
                  This assembly represents a single individual from the AGSC and was generated using 600 Gb of HiSeq paired end reads and 640Gb of HiSeq mate pair reads. Reads were assembled using a modified version of SparseAssembler [Ye C, et al. 2012].

                  I might give SparseAssembler a try.


                  • #84
                    Hi Brian,

                    I used bbmerge and am now trying to error correct my paired and merged reads with tadpole at the same, but I can't seem to get the right syntax for the input

                    I tried the following like the example,

                     ~/bin/bbmap-37.56/ in=SRR2027504_1.fq.gz,SRR2027504_merged.fq.gz in2=SRR2027504_2.fq.gz,null out=ecc_SRR2027504_1.fq.gz,ecc_SRR2027504_merged.fq.gz out2=ecc_SRR2027504_2.fq.gz,null mode=correct
                    but get,

                    Tadpole version 37.56
                    Exception in thread "main" java.lang.RuntimeException: Can't read file 'null'
                            at shared.Tools.testInputFiles(
                            at shared.Tools.testInputFiles(
                            at assemble.Tadpole.<init>(
                            at assemble.Tadpole1.<init>(
                            at assemble.Tadpole.makeTadpole(
                            at assemble.Tadpole.main(
                    What am I doing wrong?


                    • #85
                      Appears that tadpole is not able to take SE,PE reads at the same time.


                      • #86
                        I've a question.

                        Originally posted by Brian Bushnell View Post
               in=reads.fq out=merged.fq outu=unmerged.fq ihist=ihist.txt extend2=20 iterations=10 k=31 ecct qtrim2=r trimq=12 strict
                        The unmerged pairs, are they trimmed or not? Do I have to do a quality trimming again on the unmerged pairs? The same question when using adapter removing during merge.


                        • #87
                          Hi Brian,

                          I recreated the paired-end FASTQ files, performed adapter and quality trimming with bbduk, then used tadpole for error correction, and finally used tadpole in contig mode to de novo assemble the contigs within 72 hours. I used prefilter=2 and evaluated various values of K.

                          Thank you for your help. This was the only de novo assembler that I tried that could finish within 72 hours, and yes it needed up to 1.4TB RAM and 48 cores.


                          • #88
                            Hi Brian,

                            After reading this whole thread, I still have some doubts about how the mode=extend works in tadpole.

                            My understanding is: kmers of size k are extracted from the reads and, upon overlap, reads are extended. My expectation was that they were also merged together. Instead, I get the same number of reads in input and output, just extended.

                            I am ok with the output but I would like some rationale to justify it in my workflow. Such extended reads should only be used in the context of assembly, right? Because I won't try to extract a positional coverage from them, as they are an extended version of themselves.


                            • #89

                              I am getting a memory error with the following:
                     in=../fasta/2.fasta out=tad2.fa k=96 merge=t overwrite=t

                              The fasta is interleaved and has 57767162 reads. This is a metagenome file. Reads are 150 bp paired illumina novaseq, qc'd and clipped.

                              I have tried the -Xmx50g but it made no difference. I have 64GB RAM (about 62GB available, on Ubuntu 16.04) and 64GB swap, but the program doe snot seem to use the swap at all.

                              Thanks for any help.



                              Executing assemble.Tadpole2 [in=../fasta/2.fasta, out=tad2.fa, k=96, merge=t, overwrite=t, -Xmx50g]
                              Version 37.88 [in=../fasta/2.fasta, out=tad2.fa, k=96, merge=t, overwrite=t, -Xmx50g]

                              Using 8 threads.
                              Executing ukmer.KmerTableSetU [in=../fasta/2.fasta, out=tad2.fa, k=96, merge=t, overwrite=t, -Xmx50g]

                              Ways=31, initialSize=128000, prefilter=f, prealloc=f
                              Memory: max=51450m, free=50913m, used=537m

                              Initialization Time: 0.032 seconds.

                              Loading kmers.

                              Estimated kmer capacity: 585441055
                              After table allocation:
                              Memory: max=51450m, free=50376m, used=1074m

                              java.lang.OutOfMemoryError: Java heap space
                              at shared.KillSwitch.allocLong2D(
                              at ukmer.AbstractKmerTableU.allocLong2D(
                              at ukmer.HashArrayU1D.resize(
                              at ukmer.HashArrayU1D.incrementAndReturnNumCreated(
                              at ukmer.HashBufferU.dumpBuffer_inner(
                              at ukmer.HashBufferU.dumpBuffer(
                              at ukmer.HashBufferU.incrementAndReturnNumCreated(
                              at ukmer.KmerTableSetU$LoadThread.addKmersToTable(
                              at ukmer.KmerTableSetU$

                              This program ran out of memory.
                              Try increasing the -Xmx flag and using tool-specific memory-related parameters.


                              • #90
                                My mistake - I should have increased Xmx.. not decreased it. It still seems to be stalling though RAM stays full but cpu activity drops to zero.. assembly never finishes.



                                Latest Articles


                                • seqadmin
                                  Advanced Tools Transforming the Field of Cytogenomics
                                  by seqadmin

                                  At the intersection of cytogenetics and genomics lies the exciting field of cytogenomics. It focuses on studying chromosomes at a molecular scale, involving techniques that analyze either the whole genome or particular DNA sequences to examine variations in structure and behavior at the chromosomal or subchromosomal level. By integrating cytogenetic techniques with genomic analysis, researchers can effectively investigate chromosomal abnormalities related to diseases, particularly...
                                  09-26-2023, 06:26 AM
                                • seqadmin
                                  How RNA-Seq is Transforming Cancer Studies
                                  by seqadmin

                                  Cancer research has been transformed through numerous molecular techniques, with RNA sequencing (RNA-seq) playing a crucial role in understanding the complexity of the disease. Maša Ivin, Ph.D., Scientific Writer at Lexogen, and Yvonne Goepel Ph.D., Product Manager at Lexogen, remarked that “The high-throughput nature of RNA-seq allows for rapid profiling and deep exploration of the transcriptome.” They emphasized its indispensable role in cancer research, aiding in biomarker...
                                  09-07-2023, 11:15 PM





                                Topics Statistics Last Post
                                Started by seqadmin, Yesterday, 09:38 AM
                                0 responses
                                Last Post seqadmin  
                                Started by seqadmin, 09-27-2023, 06:57 AM
                                0 responses
                                Last Post seqadmin  
                                Started by seqadmin, 09-26-2023, 07:53 AM
                                1 response
                                Last Post seed_phrase_metal_storage  
                                Started by seqadmin, 09-25-2023, 07:42 AM
                                0 responses
                                Last Post seqadmin