Hey all,
I'm trying to merge my reads (300bp, demultiplexed) using BBMerge and for my samples, sequencing results of the fungal ITS2 region, BBMerge yields much lower merged reads compared to PEAR (~50% vs ~90%, respectively). Most of the unmerged reads by BBMerge are ambiguous. Few questions:
Is there something I can try with BBMerge?
I've read here that PEAR results in a lot of FP, but can it be up to 50% of the sample?
Many of my reads are longer than the actual length of the fragment sequenced (ITS2 region vary in size, 100-500bp), is BBMerge sensitive to post-adapters "garbage" bases?
I'd rather not to trim my reads based on quality \ primers' sequence since it interferes with follow steps of the analysis.
Thanks,
Omer
I'm trying to merge my reads (300bp, demultiplexed) using BBMerge and for my samples, sequencing results of the fungal ITS2 region, BBMerge yields much lower merged reads compared to PEAR (~50% vs ~90%, respectively). Most of the unmerged reads by BBMerge are ambiguous. Few questions:
Is there something I can try with BBMerge?
I've read here that PEAR results in a lot of FP, but can it be up to 50% of the sample?
Many of my reads are longer than the actual length of the fragment sequenced (ITS2 region vary in size, 100-500bp), is BBMerge sensitive to post-adapters "garbage" bases?
I'd rather not to trim my reads based on quality \ primers' sequence since it interferes with follow steps of the analysis.
Thanks,
Omer
Comment