![]() |
![]() |
#1 |
Senior Member
Location: Australia Join Date: Jul 2013
Posts: 119
|
![]()
I am trying to merge two illumina reads set, R1 and R2. However R1 has one particular read named read1 which can NOT be found in R2.
The file after merging only contains part of read1. So how to use PEAR to do a merge which include all the sequence of read1 in the merged file? Thanks, |
![]() |
![]() |
![]() |
#2 |
Member
Location: Brazil Join Date: Dec 2009
Posts: 24
|
![]()
I like to use flash2 to do it (https://github.com/dstreett/FLASH2).
These software assume both files contain the exact same number of reads. You should manage to get rid of unpaired reads. Please find attached a script that can solve that. -Cheers |
![]() |
![]() |
![]() |
#3 | |
Senior Member
Location: Australia Join Date: Jul 2013
Posts: 119
|
![]() Quote:
Thanks. Looks like flash2 do the similar job of PEAR. |
|
![]() |
![]() |
![]() |
#4 |
Senior Member
Location: East Coast USA Join Date: Feb 2008
Posts: 7,087
|
![]()
There are many programs that can do read merging. BBMerge.sh from BBMap is also user friendly.
|
![]() |
![]() |
![]() |
Tags |
illumina, pear |
Thread Tools | |
|
|