Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Running MetaVelvet 1.2.01

    Hello,

    I am assembling 6 datasets with MetaVelvet on a computer that has 128 GB of memory with huge pages enabled in Linux. The purpose is to draw a comparison with Ray Méta.

    I am using MetaVelvet 1.2.01 and Velvet 1.1.06, as indicated in the documentation available at http://metavelvet.dna.bio.keio.ac.jp/

    I compiled Velvet with MAXKMERLENGTH=51.


    I assembled

    ERS006494 186073978 sequences
    ERS006497 161222460 sequences
    ERS006592 26934980 sequences

    without errors in the process.


    However, I am not enable to assemble these:

    SRS011098 202487723 seq.

    Reason: velvetg: Can't calloc 18446744072008167979 Annotations totalling 18446744046487413424 bytes: Cannot allocate memory


    The range of an unsigned 64-bit integer is 0 to 18446744073709551615.
    Therefore, there is a integer overflow.


    SRS017227 139002751 seq.

    Segmentation fault (core dumped)

    (gdb) bt
    #0 0x0000000000429f71 in newPreGraph_pg ()
    #1 0x000000000040293c in main ()



    SRS018661 288475194 seq.

    Segmentation fault (core dumped)

    (gdb) bt
    #0 0x0000000000426940 in importRoadMapArray ()
    #1 0x000000000040292f in main ()


    Are there things in particular I should check in the input data ?

    Are these problems known bugs that were fixed since v1.1.06 ?

    Is MetaVelvet 1.2.01 compatible with Velvet 1.2.07 ?


    If it can help, my commands for SRS011098 are


    velveth SRS011098/velvet-block 51 -fastq -shortPaired SRS011098/block-device.fastq

    velvetg SRS011098/velvet-block -exp_cov auto -ins_length 211

    meta-velvetg SRS011098/velvet-block -ins_length 211


    Thank you for your support.


    Sébastien Boisvert

  • #2
    I am not an expert, but I think the error message suggests that you don't have enough memory to do the assembly?

    WJ

    Comment


    • #3


      Is this relevant? Also, if you have to use a page file (i.e. hard drive) for extra memory, it will probably mean the assembly takes far too long because reading into/out of a hard drive is about 100X slower.

      Comment


      • #4
        I was having a similar problem but my error message was coming up at the coverage estimation stage. It seems that MetaVelvet doesn't compile so well with a kmer limit above 64. If you recompile MetaVelvet and make sure you set your kmer limit to 64, you might solve the issue. It worked for me at least.

        Comment


        • #5
          63. Velvet doesn't like even numbered kmer values (something to do with it not dealing well with palindromes). So you'd never use 64 as an actual kmer value.

          Comment


          • #6
            Yep sorry that was a typo. I meant 63.

            Comment


            • #7
              Hi Sébastien Boisvert,

              Currently, the last version of MetaVelvet is MetaVelvet version 1.2.02 which is compatible with the last version of Velvet, version 1.2.10.
              Anyway, for the error in your case, you need more memory and also the integer type is not enough to process your data (we will revise the type). Thanks a lot for using and trying MetaVelvet.

              Best regards,
              Afia (currently continue the MetaVelvet project)

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Best Practices for Single-Cell Sequencing Analysis
                by seqadmin



                While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                06-06-2024, 07:15 AM
              • seqadmin
                Latest Developments in Precision Medicine
                by seqadmin



                Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                Somatic Genomics
                “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                05-24-2024, 01:16 PM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Today, 07:24 AM
              0 responses
              9 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 08:58 AM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-12-2024, 02:20 PM
              0 responses
              16 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-07-2024, 06:58 AM
              0 responses
              182 views
              0 likes
              Last Post seqadmin  
              Working...
              X