SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
200-250bp read length cumulonix Illumina/Solexa 22 05-05-2016 06:46 AM
splitting longer reads into two smaller reads for bowtie input a_mt Bioinformatics 0 02-19-2013 09:14 PM
250bp paired end Miseq reads alignment empyrean Bioinformatics 3 12-12-2012 12:36 PM
duplicate reads in Illumina short, single end reads of RNAseq data inbarpl Bioinformatics 4 05-22-2012 08:36 AM
Alignment of ABI solid reads and 454 reads baohua100 Bioinformatics 2 02-23-2009 04:58 PM

Reply
 
Thread Tools
Old 04-08-2013, 01:21 PM   #1
kcchan
Senior Member
 
Location: USA

Join Date: Jul 2012
Posts: 182
Default 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

Quote:
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?
kcchan is offline   Reply With Quote
Old 04-09-2013, 01:36 AM   #2
akinori@jp
Junior Member
 
Location: Tokyo, Japan

Join Date: Apr 2013
Posts: 2
Smile

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 at 02:19 AM.
akinori@jp is offline   Reply With Quote
Old 04-09-2013, 09:12 AM   #3
kcchan
Senior Member
 
Location: USA

Join Date: Jul 2012
Posts: 182
Default

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?
kcchan is offline   Reply With Quote
Old 04-09-2013, 09:59 AM   #4
akinori@jp
Junior Member
 
Location: Tokyo, Japan

Join Date: Apr 2013
Posts: 2
Default

I checked it under my environments.
It seemed to be running without any errors.
I've not changed any other points......
akinori@jp is offline   Reply With Quote
Old 04-10-2013, 05:26 PM   #5
ewilbanks
Member
 
Location: Davis, CA

Join Date: Mar 2009
Posts: 82
Default

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?
ewilbanks is offline   Reply With Quote
Old 09-02-2014, 08:11 PM   #6
tsucheta
Member
 
Location: Arlington, TX

Join Date: Nov 2009
Posts: 17
Default

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
tsucheta is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:57 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO