Question: how to merge pair-end reads with maximum score 41 (J)
0
gravatar for tiantai_wing
10 weeks ago by
tiantai_wing0 wrote:

I want to merge my pair-end reads into one single reads with about 50 bp overlap. However, my fastq file is Illumina-1.8 Phred+33 format, so the maximum score is 41 (J). I tried NGmerge for merging but it can not expand to score 41 with -u 41 and editing the qual_profile. Is there any other method available to merge this kind of data? Need your help. Thank you!

next-gen assembly sequence • 136 views
ADD COMMENTlink written 10 weeks ago by tiantai_wing0
1

Have you tried bbmerge.sh from BBMap suite? A guide is available here.

ADD REPLYlink written 10 weeks ago by genomax73k
1

Any paired-end read merger (FLASH, PEAR, BBMerge, and several others) will merge Illumina-1.8 (Phred+33) fastq files without problems, as this is the current standard Illumina format, and by far the most common enconding found in the databases.

What is the problem you are having? What do you expect when you merge R1 and R2? It seems you expect the score of merged bases will be higher than 41, but this is not happening - is this the problem?

ADD REPLYlink written 10 weeks ago by h.mon27k

Thanks all! bbmerge is so fast (merging in 2 min for 7 million paired-reads)!

ADD REPLYlink written 10 weeks ago by tiantai_wing0
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: 2007 users visited in the last hour