Question: Identifying shared exact sequence variants between two datasets
gravatar for adam.sorbie
8 days ago by
adam.sorbie0 wrote:

I have two microbiome datasets which are sequenced from different sites within the environment. I would to like to identify amplicon sequence variants inferred using DADA2 which are shared between the two datasets and also those which are unique. The problem is i'm not 100% sure of the best way to do this. So far I have tried megablast and just looking at exact matches between the two datasets but perhaps there is a better way of doing this, for example another alignment tool or alignment free method but unfortunately I don't have the experience to be able to decide on the best method. Does anyone with more experience of working with microbiome data know of a better way?

microbiome ngs alignment 16s dada2 • 114 views
ADD COMMENTlink modified 4 days ago by genebow140 • written 8 days ago by adam.sorbie0

"amplicon sequence variants " - you mean the common and different variants between two datasets ? Does your organism has reference genome available ?

ADD REPLYlink modified 6 days ago • written 6 days ago by geek_y8.8k

The term is just what dada2 use to refer to each individual 16S rRNA sequence representing a "species" of bacteria, they are not variants in the genetics sense. Essentially yes, I wish to find the common and differing sequences between the two datasets.

ADD REPLYlink written 5 days ago by adam.sorbie0
gravatar for h.mon
4 days ago by
h.mon22k wrote:

You can use cd-hit-est-2d for this task. The VSEARCH wiki also has an example on how to use it to (among other things) compare datasets.

ADD COMMENTlink written 4 days ago by h.mon22k
gravatar for genebow
4 days ago by
genebow140 wrote:

You may use sliding window approach to align windowed sequence from one genome onto the other genome.

ADD COMMENTlink written 4 days ago by genebow140

OP data is amplicon, probably 16S.

ADD REPLYlink written 4 days ago by h.mon22k
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: 1201 users visited in the last hour