![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
SOAPdenovo | vani s kulkarni | Bioinformatics | 6 | 06-17-2013 01:08 PM |
SOAPdenovo v1.05 | scami | Bioinformatics | 1 | 01-05-2013 12:04 AM |
SOAPdenovo -error | robelb4 | Bioinformatics | 0 | 08-05-2011 12:39 PM |
SOAPdenovo assembler | shailesh | Bioinformatics | 6 | 10-04-2010 07:35 PM |
SOAPdenovo | AnthonyB | Bioinformatics | 2 | 05-30-2009 07:12 PM |
![]() |
|
Thread Tools |
![]() |
#1 |
Member
Location: Virginia Join Date: Mar 2011
Posts: 72
|
![]()
Hello.
I am wondering if anyone here some some experience with SOAPdenovo. I am trying to assembly Illumina data, I have 3 lanes of 101 cycle PE data and 2 lanes of single read data. I have tried all sorts of options for running the program down to the default ./SOAPdenovo63mer all -s SOAPdenovo_config.txt -o out_ At the end of the day (ok, only 2 hours), it is telling me 'no paired reads found' and tossed a terminal error. I have used this same data in both CLC and Velvet and it works fine including scaffolding in Velvet and paired mapping in CLC. I am sure I am doing something silly. Any thoughts? Thanks, Bob ps. my config file is: #maximal read length max_rd_len=100 [LIB] #average insert size avg_ins=200 #if sequence needs to be reversed reverse_seq=0 #in which part(s) the reads are used asm_flags=1 #in which order the reads are used while scaffolding rank=1 p=/common/groups/Illx/Genomic_SIPES_lane5_1_paired_trimmed_paired.fastq p=/common/groups/Illx/Genomic_SIPES_lane6_1_paired_trimmed_paired.fastq p=/common/groups/Illx/Genomic_SIPES_lane7_1_paired_trimmed_paired.fastq q=/common/groups/Illx/Unpaired_SIPES_PF_lane5_1_trimmed.fastq q=/common/groups/Illx/Unpaired_SIPES_PF_lane7_1_trimmed.fastq |
![]() |
![]() |
![]() |
#2 |
Member
Location: Netherlands Join Date: Sep 2008
Posts: 13
|
![]()
It has been a while since I looked at SOAPdenovo, but could you try this config ?
---- #maximal read length max_rd_len=100 [LIB] #average insert size avg_ins=200 #if sequence needs to be reversed reverse_seq=0 #in which part(s) the reads are used asm_flags=1 #in which order the reads are used while scaffolding rank=1 q=/common/groups/Illx/Unpaired_SIPES_PF_lane5_1_trimmed.fastq q1=/common/groups/Illx/Unpaired_SIPES_PF_lane7_1_trimmed.fastq [LIB] #average insert size avg_ins=200 #if sequence needs to be reversed reverse_seq=0 #in which part(s) the reads are used asm_flags=3 #in which order the reads are used while scaffolding rank=2 q=/common/groups/Illx/Genomic_SIPES_lane5_1_paired_trimmed_paired.fastq q1=/common/groups/Illx/Genomic_SIPES_lane6_1_paired_trimmed_paired.fastq q2=/common/groups/Illx/Genomic_SIPES_lane7_1_paired_trimmed_paired.fastq |
![]() |
![]() |
![]() |
#3 |
Member
Location: Virginia Join Date: Mar 2011
Posts: 72
|
![]()
Thanks.
Unfortunately, it appears as though it was a case of RTFM more carefully. Paired reads have to be in fast-a format if in a single file. I only wasted 2 days because I didn't see that earlier. =/ |
![]() |
![]() |
![]() |
Thread Tools | |
|
|