Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 250bp reads in idba_ud

    Has anyone experimented with using idba_ud on 2x250bp MiSeq data? I'm getting the following error when I tried to run it using -r and -l

    terminate called after throwing an instance of 'std::logic_error'
    what(): SequenceReader::SequenceReader() istream is invalid
    Aborted (core dumped)
    The assembler works fine when the dataset is trimmed back to 150bp. Is there some configuration that I'm missing?

  • #2
    I think 250-bp is too long to assemble at least with -r option.
    In default settings, the max is set to be 128 bp.

    To change the max value, you have to re-compile IDBA from the source as follows:

    1. change the vaue 0f 128 to 250 (or more) in /src/sequence/short_sequence.h

    static const uint32_t kMaxShortSequence = 128;
    ->
    static const uint32_t kMaxShortSequence = 250;

    2. overwrite the short_sequence.h

    3. configure & make!

    In my case, it works well with 2x250bp miseq-paired-end reads by setting 512
    Last edited by akinori@jp; 04-09-2013, 02:19 AM.

    Comment


    • #3
      Thanks for the info, re-compiling IDBA with that change worked. It's strange that it doesn't work when I use the -l parameter though, isn't that the whole point of the option?

      Comment


      • #4
        I checked it under my environments.
        It seemed to be running without any errors.
        I've not changed any other points......

        Comment


        • #5
          Are you folks wokring from idba-1.1.0 or 1.0.9? I made the change above and it is altering only the max read length accepted, but not changing the maxK! When I try running it doesn't run saying my maxK is too large:

          $ idba_ud -r MiSeq_260bp.fa -o test --maxk 250


          Would either of you mind posting your full short_sequence.h files?

          Comment


          • #6
            I have been trying to run IDBA for sometimes now. I have illumina paired end reads which I merged together for testing into one file and they are 151 bases in length. I have had the same error as Kcchan and I increased the static const uint32_t kMaxShortSequence = 1600 and recompiled. It now runs till it hits k=50 and then it exits with the following error:

            kmer 50
            kmers 9815948 9845052
            merge bubble 519
            contigs: 7985 n50: 4122 max: 71069 mean: 920 total length: 7348632 n80: 919
            reads 9449880
            aligned 8473710 reads
            distance mean 229.862 sd 12592.2
            invalid insert distance
            ./command.sh: line 4: 14229 Segmentation fault /share/apps/idba-1.1.1/bin/i
            dba -r test.fa -o test_idba_out

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Current Approaches to Protein Sequencing
              by seqadmin


              Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
              04-04-2024, 04:25 PM
            • 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

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 12:08 PM
            0 responses
            11 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 10:19 PM
            0 responses
            17 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 09:21 AM
            0 responses
            14 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-04-2024, 09:00 AM
            0 responses
            43 views
            0 likes
            Last Post seqadmin  
            Working...
            X