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
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
Comment