Entering edit mode
9.1 years ago
rebzi87
•
0
Hi -
I'm trying to do a de novo assembly using Illumina paired-end reads and 3 mate-pair libraries of different insert length. I use the command:
abyss-pe np=32 k=61 name=wagtail lib='pe1' mp='mp1 mp2 mp3 mp4 mp5 mp6 mp7 mp8 mp9' \
pe1='cutadapt1.fq.gz cutadapt2.fq.gz' mp1='2kb-out_A_R1.fasta.gz 2kb-out_A_R2.fasta.gz' \
mp2='2kb-out_B_R1.fasta.gz 2kb-out_B_R2.fasta.gz' mp3='2kb-out_C_R1.fasta.gz 2kb-out_C_R2.fasta.gz' \
mp4='4kb-out_A_R1.fasta.gz 4kb-out_A_R2.fasta.gz' mp5='4kb-out_B_R1.fasta.gz 4kb-out_B_R2.fasta.gz' \
mp6='4kb-out_C_R1.fasta.gz 4kb-out_C_R2.fasta.gz' mp7='9kb-RCout_A_R1.fasta.gz 9kb-RCout_A_R2.fasta.gz' \
mp8='9kb-RCout_B_R1.fasta.gz 9kb-RCout_B_R2.fasta.gz' mp9='9kb-RCout_C_R1.fasta.gz 9kb-RCout_C_R2.fasta.gz'
where each mate-pair library has group A, B, and C as output by the quality filtering program NextTrim. I keep getting the following error:
Building the suffix array...
terminate called after throwing an instance of 'std::bad_alloc'
what(): std::bad_alloc
Mateless 0
Unaligned 0
Singleton 0
FR 0
RF 0
FF 0
Different 0
Total 0
abyss-fixmate: error: All reads are mateless. This can happen when first and second read IDs do not match.
error: `pe1-3.hist': No such file or directory
make: *** [pe1-3.dist] Error 1
make: *** Deleting file `pe1-3.dist'
I saw a past post about the naming of paired reads and changed the files to have the /1 /2 naming format (instead of a space) so I changed my fasta file names from:
>D00723:25:C6C71ANXX:4:1101:1748:2147 1:N:0:14
>D00723:25:C6C71ANXX:4:1101:1748:2147/1
I still get the same errors and am at a loss of what to do. Any help would be greatly appreciated! Thanks.
It may be because your left and right reads are very different sizes. I got the same error after Trimmomatic discarded most of my left reads.