Rarefaction/saturation curve based on NGS data
Hi all,
most likely a simple question, but I'm looking for a tool (software, python/Perl/R script) that produces a rarefaction curve base on an input assembly (ACE would be easiest) to assess the nr reads needed to yield all observed contigs (cfr species diversity index). This would most likely be done through resampling of reads within the ACE file. I would like to compare this for normalized and non normalized libraries.
thanks
most likely a simple question, but I'm looking for a tool (software, python/Perl/R script) that produces a rarefaction curve base on an input assembly (ACE would be easiest) to assess the nr reads needed to yield all observed contigs (cfr species diversity index). This would most likely be done through resampling of reads within the ACE file. I would like to compare this for normalized and non normalized libraries.
thanks
