Hi,
I was asked to provide tools in a local galaxy instance to handle the following research question:
Two groups of microbial cultures taken from case and control animals, were 16S sequences on 2x250bp miseq. The aim is to identify species present in case and absent in control.
So we need to merely identify species present, and optionally some sort of quantification, no clustering and so on.
As my background is human exome seq, I would welcome advice on how to approach this. Right now, I have a blastn, followed by selecting the best hit for each read in mind.
Thanks!
I was asked to provide tools in a local galaxy instance to handle the following research question:
Two groups of microbial cultures taken from case and control animals, were 16S sequences on 2x250bp miseq. The aim is to identify species present in case and absent in control.
So we need to merely identify species present, and optionally some sort of quantification, no clustering and so on.
As my background is human exome seq, I would welcome advice on how to approach this. Right now, I have a blastn, followed by selecting the best hit for each read in mind.
Thanks!
Comment