Question: How to compare coverage eveness of two whole genome sequencing runs of the same sample?
gravatar for xuang
2.8 years ago by
xuang0 wrote:

We are comparing two library prep kits: Nextera DNA Flex Library Prep and TruSeq DNA PCR-Free Library Prep for whole genome sequencing. The purpose of this comparison is to predict how will they affect the result of CNV calling.

So we whole-genome-sequenced the same sample prepped with either of the two preps, mapped them to the reference genome with bwa mem.

Now we have the general mapping statistics but that is not enough. To evaluate how CNV calling will be affected, our fist idea is to evaluate the eveness of coverage, but don't know to evaluate that properly.

Any ideas or existing tools to compare the coverage?

sequencing cnv library prep wgs • 867 views
ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by xuang0
gravatar for Devon Ryan
2.8 years ago by
Devon Ryan98k
Freiburg, Germany
Devon Ryan98k wrote:

One simple method would be to run multiBamSummary (or make normalized bigWig files with bamCoverage and then multiBigwigSummary), after which you can use tools like plotCorrelation and plotPCA to see how well the (possibly normalized) coverage compares across bins in your dataset.

ADD COMMENTlink written 2.8 years ago by Devon Ryan98k

In this case which normalization method would you recommend for WGS

ADD REPLYlink written 23 months ago by rrdavis20

RPGC (aka 1X) normalization would be the simplest.

ADD REPLYlink written 23 months ago by Devon Ryan98k
gravatar for Carlo Yague
2.8 years ago by
Carlo Yague5.5k
Carlo Yague5.5k wrote:

There is a rather simple tool designed to do exactly what you want: cnv-seq. Its a bit old but rather straightforward. I used it in the past and as far as I remember, it was easy to use and it made nice figures.

ADD COMMENTlink written 2.8 years ago by Carlo Yague5.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1597 users visited in the last hour