Question: Identifying shared exact sequence variants between two datasets
0
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
2
gravatar for h.mon
4 days ago by
h.mon22k
Brazil
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
0
gravatar for genebow
4 days ago by
genebow140
USA/Chicago
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
1

OP data is amplicon, probably 16S.

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