Question: How to merge quality trimmed interlaced reads?
0
gravatar for O.rka
12 weeks ago by
O.rka40
O.rka40 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
GTGGTCAGCAGACGTTTAGCTTCGTCAACCAGGTCAGCTTCGTACAGGGATTTACAGATTTCAGGG
>NS500647:208:HYKFFBGX2:1:11101:9580:1154 2:N:0:GATCAG
ACTGGNNNTTCTGGAAAANNNNNGNNTCANC
>NS500647:208:HYKFFBGX2:1:11101:24675:1156 1:N:0:GATCAG
AGCCGATATTCACTACCTGCTCGCCTTTAACGTTCGCAATGTCTTTTAGTTGCGGCACCGCATTAACCAGA
>NS500647:208:HYKFFBGX2:1:11101:24675:1156 2:N:0:GATCAG
CGCACNNNATATGGGTTTTANTGGNGCAGCAT

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 • 197 views
ADD COMMENTlink written 12 weeks ago by O.rka40
1

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

ADD REPLYlink written 12 weeks ago by Charles Plessy2.5k
1

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 12 weeks ago by mbk0asis350

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 12 weeks ago by O.rka40
1

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 bbmerge.sh 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 12 weeks ago • written 12 weeks ago by genomax50k
Please log in to add an answer.

Help
Access

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