Question: How to merge quality trimmed interlaced reads?
gravatar for O.rka
15 months ago by
O.rka110 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 • 758 views
ADD COMMENTlink modified 6 months ago by Biostar ♦♦ 20 • written 15 months ago by O.rka110

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 15 months ago by mbk0asis420

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 15 months ago by O.rka110

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 15 months ago • written 15 months ago by genomax68k

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 6 months ago by kristoffer.vittingseerup1.9k

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

ADD REPLYlink written 15 months 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: 830 users visited in the last hour