I have been using FastX to process reads prior to de novo assembly and mapping. What I have discovered and few have pointed out is the FastX will delete reads leaving reads unpaired which changes the order of the separate paired fastq files. While it is difficult to know if this is affecting assembly with Trinity, it is definitely a problem for assembly with Velvet/Oases and mapping with Bowtie or BWA. Because the order of the paired reads has changed due to deletions of low quality reads, the reads are no longer order properly and will not map as paired.
There are some work arounds provided by sfg.stanford.edu and others to separate the reads that are still paired and place the orphaned reads in a separate file. But here is the problem, I would like to use paired reads in combination with single reads for de novo assembly. In Trinity, one designates as --right -- left or --singles, but you cannot do both.
Question: Can any assembler use both paired and single reads at the same time for de novo assembly?
Q2: Has anyone else run into this problem? Here is a related post: http://seqanswers.com/forums/showthread.php?t=24076
Q3: This issue is going to eliminate FastX from my pipeline of assembly and mapping. It seems like this should be a bigger issue but there is fairly little out there about this. Am I doing something wrong with FastX that is causing this problem?
Thanks.
Hi, have you sorted that out? Im using trinity and i'm struggling with the same problem: I want to use my single data (merged paired end reads) and left-right (unpaired reads) together for a big assembly that I want to use as a reference, if not I loose a lot of data.