Merged reads in SPAdes 3.14
0
1
Entering edit mode
19 months ago
lars.snipen ▴ 30

I have paired-end Illumina reads, R1.fq and R2.fq. From these I merge overlapping reads, using BBmerge, and get the file merged.fq. In SPAdes I can input both paired-end (options --pe1-1 and --pe1-2) files and the merged.fq file (option --pe1-m). My question is: Should the paired-end files be the _original files_, containing all reads, or only the _leftover files_ with reads that did not merge? From what I read in the manual I would say the latter, but it is not entirely clear. One may ask why bother? Well, it might be a memory issue, since these samples have short insert size, and a large fraction of the reads merge .Thus, providing them twice (both in paired-end and merged files) seems like a waste of memory.

Assembly • 490 views
ADD COMMENT
0
Entering edit mode

Should the paired-end files be the _original files_, containing all reads, or only the _leftover files_ with reads that did not merge?

I would not worry about merging the files since SPAdes should be able to figure that out. Just try the assembly with the original reads as is initially.

ADD REPLY

Login before adding your answer.

Traffic: 1559 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6