16S metagenomics data normalization
1
0
Entering edit mode
4.0 years ago
junaidkori • 0

I have 16S amplicon metagenomics data i want to normalize my reads because number of reads varies between different samples. Can any body help me.

sequencing genome sequence R • 2.1k views
ADD COMMENT
0
Entering edit mode
4.0 years ago

If you use QIIME then you could summarize your biom table, and choose appropriate depth through '-e' parameter in core diversity analyses. This accounts for differences in depth between samples. To choose this value you could simply create a rarefaction plot again in QIIME.

Just read on:

http://nbviewer.jupyter.org/github/biocore/qiime/blob/1.9.1/examples/ipynb/illumina_overview_tutorial.ipynb

and

http://qiime.org/scripts/make_rarefaction_plots.html

I usually sumamrize table and take the "sensible" lowest value. I.e. if the values are more or less for ex. 1000, and there is one around 150, I choose 1000, therefore excluding a particularly shallow sample and ensuring further normalization.

There may be other methods though.

ADD COMMENT
2
Entering edit mode
ADD REPLY

Login before adding your answer.

Traffic: 1606 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6