Is there a paired end joiner which writes also the reads before merging?

rexxi · 03-07-2019, 03:54 PM

Hi, I have a weird question. But I'm looking for a software which not just merge the paired reads but also writes them in a new file.

When a software matches a pair of reads it writes them in the output fastq (or fasta), but you don't have the option of knowing what they paired. I've checked some of them but mostly they create a file with the unpaired reads as the set of the discarded ones. I've been trying to find one which would give an output like that

SAMPLE_MATCHED_FORWARD.fastq
SAMPLE_MATCHED_REVERSE.fastq
SAMPLE_JOINED.fastq

the first two files would contain the reads that are going to be merged, but have been already found their mate. Is there a software with an option like that? I know it sounds weird but it could help me to discriminate from a pool containing a lot of unwanted sequences from different sources.

thanks for your time.

torben · 03-11-2019, 02:46 PM

If it's not a problem to run the command twice, bbmerge can do this:

Code:

bbmerge.sh in1=SAMPLE_INPUT_FORWARD.fastq in2=SAMPLE_INPUT_REVERSE.fastq
out1=SAMPLE_MATCHED_FORWARD.fastq out2=SAMPLE_MATCHED_REVERSE.fastq merge=f

Code:

bbmerge.sh in1=SAMPLE_INPUT_FORWARD.fastq in2=SAMPLE_INPUT_REVERSE.fastq
out=SAMPLE_JOINED.fastq merge=t

Run bbmerge.sh without options for a list of possible options. See also this thread.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Merging paired end reads for BLAST	JJenks	Bioinformatics	9	11-05-2018 10:40 AM
A question on Illumina paired-end reads alignment. Merging from different samples	linudz	Bioinformatics	4	02-27-2017 11:16 AM
Merging non-overlapping paired end reads	karenr	Illumina/Solexa	9	12-16-2016 07:02 PM
Merging illumina V4 paired end reads	rEDI	Illumina/Solexa	8	09-09-2016 09:10 AM
Merging paired end reads (R1 and R2 files)	vectorborne5	Illumina/Solexa	2	10-09-2013 11:13 AM

03-07-2019, 03:54 PM	#1
rexxi Member Location: Santiago, Chile Join Date: Jun 2012 Posts: 20	Is there a paired end joiner which writes also the reads before merging? Hi, I have a weird question. But I'm looking for a software which not just merge the paired reads but also writes them in a new file. When a software matches a pair of reads it writes them in the output fastq (or fasta), but you don't have the option of knowing what they paired. I've checked some of them but mostly they create a file with the unpaired reads as the set of the discarded ones. I've been trying to find one which would give an output like that SAMPLE_MATCHED_FORWARD.fastq SAMPLE_MATCHED_REVERSE.fastq SAMPLE_JOINED.fastq the first two files would contain the reads that are going to be merged, but have been already found their mate. Is there a software with an option like that? I know it sounds weird but it could help me to discriminate from a pool containing a lot of unwanted sequences from different sources. thanks for your time.

03-11-2019, 02:46 PM	#2
torben Member Location: Norway Join Date: Oct 2012 Posts: 20	If it's not a problem to run the command twice, bbmerge can do this: Code: bbmerge.sh in1=SAMPLE_INPUT_FORWARD.fastq in2=SAMPLE_INPUT_REVERSE.fastq out1=SAMPLE_MATCHED_FORWARD.fastq out2=SAMPLE_MATCHED_REVERSE.fastq merge=f Code: bbmerge.sh in1=SAMPLE_INPUT_FORWARD.fastq in2=SAMPLE_INPUT_REVERSE.fastq out=SAMPLE_JOINED.fastq merge=t Run bbmerge.sh without options for a list of possible options. See also this thread.