Question: very small insert sizes for 16S rRNA after paired-end reads overlap
gravatar for ElUretsky
2.4 years ago by
ElUretsky20 wrote:

I have 16s rRNA reads sequenced on an Illumina MiSeq 2x250bp. Our target region is about 300bp so pairs are expected to overlap. We used leeHom to merge paired-reads that overlap. Up to now it seems that we have good results but we are also getting many pairs where there is only a barcode combined with adapters, for example:



The dash separates the adapter and the barcode. This is normally a minority of cases but there are some libraries where this is more frequent. Has else anyone experienced this? How can we prevent this for future runs?

ngs miseq 16s rrna • 738 views
ADD COMMENTlink modified 2.4 years ago by Biostar ♦♦ 20 • written 2.4 years ago by ElUretsky20

Are you referring to primer dimers where there is no real insert? You would want to control that at the library prep/purification step. Bioinformatically you can scan and trim the data to remove the dimers (use from BBMap suite, trimmomatic or cutadapt).

ADD REPLYlink written 2.4 years ago by genomax89k

they are very small so they do not make in the analysis pipeline. They question is more about how to avoid this maybe from groups doing 16s sequencing? sorry for the confusion.

ADD REPLYlink written 2.4 years ago by ElUretsky20

You would generally do bead washing to remove those. This question is more appropriate for if you want assistance with experimental strategy.

ADD REPLYlink modified 2.4 years ago • written 2.4 years ago by genomax89k
Please log in to add an answer.


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