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
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
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?
Leave a comment:
-
I checked it under my environments.
It seemed to be running without any errors.
I've not changed any other points......
Leave a comment:
-
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?
Leave a comment:
-
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 512Last edited by akinori@jp; 04-09-2013, 02:19 AM.
Leave a comment:
-
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)Tags: None
Latest Articles
Collapse
-
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...-
Channel: Articles
04-04-2024, 04:25 PM -
-
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...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
24 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
52 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Leave a comment: