Question: very small insert sizes for 16S rRNA after paired-end reads overlap
0
gravatar for ElUretsky
17 months ago by
ElUretsky20
Germany
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:

GCTCAGGA-AGATCGGAAGAACACACGTCTGAACTCCAGTCACGG....

TCCTGAGC-AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTCTC...

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 • 472 views
ADD COMMENTlink modified 16 months ago by Biostar ♦♦ 20 • written 17 months ago by ElUretsky20
1

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 bbduk.sh from BBMap suite, trimmomatic or cutadapt).

ADD REPLYlink written 17 months ago by genomax71k

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 17 months ago by ElUretsky20

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

ADD REPLYlink modified 17 months ago • written 17 months ago by genomax71k
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: 1831 users visited in the last hour