Question: Identifying shared exact sequence variants between two datasets
0
gravatar for adam.sorbie
3 months 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 • 225 views
ADD COMMENTlink modified 3 months ago by genebow150 • written 3 months 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 3 months ago • written 3 months ago by geek_y9.4k

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 3 months ago by adam.sorbie0
2
gravatar for h.mon
3 months ago by
h.mon24k
Brazil
h.mon24k 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 3 months ago by h.mon24k
0
gravatar for genebow
3 months ago by
genebow150
USA/Chicago
genebow150 wrote:

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

ADD COMMENTlink written 3 months ago by genebow150
1

OP data is amplicon, probably 16S.

ADD REPLYlink written 3 months ago by h.mon24k
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: 1513 users visited in the last hour