Question: How to merge quality trimmed interlaced reads?
gravatar for O.rka
2.4 years ago by
O.rka210 wrote:

I have paired end reads from Illumina NextSeq. I've removed adapters/primers, quality trimmed, and interlaced the forward and reverse reads. I now have a single fasta (not fastq) file with the following structure:

>NS500647:208:HYKFFBGX2:1:11101:9580:1154 1:N:0:GATCAG
>NS500647:208:HYKFFBGX2:1:11101:9580:1154 2:N:0:GATCAG
>NS500647:208:HYKFFBGX2:1:11101:24675:1156 1:N:0:GATCAG
>NS500647:208:HYKFFBGX2:1:11101:24675:1156 2:N:0:GATCAG

Is it possible to merge the R1 and R2 reads into a single fragment? Or is that not how the paired end reads work?

rna-seq • 1.5k views
ADD COMMENTlink modified 20 months ago by Biostar ♦♦ 20 • written 2.4 years ago by O.rka210

You don't want to merge them together. Paired end reads are the sequences at the both ends of a DNA fragments. When you have a 400 bp DNA fragment, your reads cover only 75 or 150 bases (NextSeq) at the end of the DNA. So if you merge PE reads into one, you will lose internal sequences. Instead, you can map the reads onto reference genomes or transcriptomes using various tools such as HISAT2, bwa, bowtie2, and etc.

ADD REPLYlink written 2.4 years ago by mbk0asis560

Thank you. I think I may have had some fundamental information incorrect about paired end sequencing technology. I was under the impression that you have overlapped segments so if you had 75 bp reads you would end up with ~150 bp reads. It seems like that is not the case. So for paired end reads, do mappers count each read of the pair individually?

ADD REPLYlink written 2.4 years ago by O.rka210

It is possible for paired-end reads to overlap in the middle (if length of sequencing is longer than the size of the fragment being sequenced, e.g. 16S amplicons are designed like that) but that will not always be the case. You can easily find out if reads in your case are able to merge into a longer representation by using from link given by @Charles.

Counting programs need to be told that you have paired-end reads (so they are not double counted) since each paired-end reads come from one fragment of DNA (mappers don't do the counting but will consider paired-ness of reads when aligning them).

ADD REPLYlink modified 2.4 years ago • written 2.4 years ago by genomax87k

you dont map them individually - using the tools suggest above by mbk0asis you can map them as paired reads meaning the tool will know they are paired and take that into account when mapping.

ADD REPLYlink written 20 months ago by kristoffer.vittingseerup3.4k

Have you looked at the software listed in the post "Tools to merge overlapping paired-end reads"?

ADD REPLYlink written 2.4 years ago by Charles Plessy2.7k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1515 users visited in the last hour