Comparing exome sequencing data sample with vastly different read coverage · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Comparing exome sequencing data sample with vastly different read coverage

I have a question about exome sequencing data. I have sample of data from the same individual taken from blood and tumour samples. The read coverage of the tumour sample is ~500X but for the blood sample it's ~100X. This coverage is fine for calling somatic mutations but in comparing the sample for CNV calls, is it a concern that the difference between read coverages is so vast? While the coverage seems adequate for both samples individually, is it possible the difference in read count will distort the estimated log ratios between the samples when they differ by several magnitudes.

Thank you for any advice you can give.

Matthew

sequencing next-gen

• 2,651 views • link

updated 4.7 years ago by Ram

• written 12.1 years ago by matthew.1.davies •

1 answer

The answer may depend on the algorithm you are using for CNV calling but, in general, most algorithms account for differing coverage between samples by "centering" the resulting copy number estimates in some way. In other words, the differences in coverage is likely not going to be problematic.

• 2 views • link

updated 4.7 years ago by Ram

• written 12.1 years ago by Sean Davis

Thank you, Sean.

• 3 views • link

written 12.1 years ago by matthew.1.davies •

Log in to answer this question.

More posts like this

Spades Genome Assembly
written 1.8 years ago by SomeOne •

Hi Reader, I am using `SPAdes 4.0.0` to assemble some fungal genomes from Illumina short read PE150 data. My data is at 100x coverage. and …
Marking Duplicate reads in deep sequencing data
written 4.2 years ago by Tom •

Hello Bio-wizards, I'm currently aligning (ultra) deep whole-exome sequencing data with a target depth of coverage of 500x on bulk (FFPE) samples. I'm going through …
Normalizing RNA-seq read counts within a single sample - gut check
written 5.5 years ago by plberry •

I'm doing a pilot analysis on some public RNA-seq datasets that include mRNA-seq, Ribo-seq, and RNC-seq for the same samples. The analysis is part of …
Seurat on both pooled and unpooled samples
written 7.5 years ago by Kristin Muench

Hello, I'm analyzing an scRNA-Seq dataset that includes both pooled and unpooled samples. Additionally, the number of cells per sample is different between the unpooled …
Bam files coverage and library size
written 7.7 years ago by David

Hi, I have aligned 10 metagenomic samples(WGS --> 2x150bp) to a bacterial reference genome to get the coverage for that specific bacterial reference. My 10 …
Difference in pipeling designing for data analysis and output files generated by Sequencing by synt…
written 8.0 years ago by shuksi1984 •

If I sequence a tissue/blood sample from the same subject (human), how will the outputs from Illumina's Genome analyzer, Roche 454's Genome Sequencer, and Applied …
CNV-Seq read count normalisation
written 8.5 years ago by nanana

I've been using [CNV-Seq][1] to detect CNV in a tumour normal pair. CNV-Seq produces 2 different files: A `.count` file, e.g.: chromosome start end test …
Multiplex In Exome Seq And "Customized/Capture" Seq Kits
written 13.1 years ago by CrazyB

<p>Does anyone have the experience of using Nimblegen exome seq kit or customized seq kit to enrich their samples before sequencing the DNA?</p> <p>My question …
What Coverage For Genome Re-Sequencing By Illumina ?
written 13.8 years ago by helene.badouin •

<p>Hello,</p> <p>I was wondering was coverage you need to do genome re-sequencing in illumina (Illumina HighSeq 2000) ?</p> <p>I was told 100x, which seems high, …
Gene Coverage In Music
written 13.8 years ago by yliubu •

<p>Hi, all,</p> <p>I am interested in MuSiC. When I use CalcWigCovg function, there is a directory called gene_covg, which contain per-gene covered base counts for …