This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Synchronization Of Pair-End Reads

After adapter trimming and quality filtering, I found some reads are too short (eg. 10bp),
so I delete all reads shorter than 25bp.
However, after this manipulation, some reads are missing from one PE library,
(eg. read1 was presented in 5'library, but missed in 3'library)
I try to re sync the PE library, but i did not find efficient methods.
I wrote a perl script to do this job, using hash table to store reads in each library, but the memory requirement is too high.
do anyone have some memory efficient methods to do this job?

paired
• 4,810 views • link
updated by brentp brentp • written by Plantae Plantae

2 answers

Have you tried using trimmomatic? I takes care about the read pairing.

• 1 views • link
written by Frenkiboy Frenkiboy

Usually, this is done by the trimmer. However, with your constraints, you could use a bloom filter (in perl) backed by a hash.

There's an example of how you'd do this here.

But, the simplest way is to use something like Sickle that does this for you as @frenkiboy suggests.

• 1 views • link
written by brentp brentp

Log in to answer this question.

More posts like this