Question: Pcr Free Variant Calls
1
gravatar for tommy.carstensen
5.3 years ago by
tommy.carstensen20 wrote:

If I have a PCR free high coverage bam and a low coverage bam and I downsample the PCR free high coverage bam to the same coverage as the low coverage bam, would I then expect the variants called from each of the two bams to be different?

variant-calling • 1.0k views
ADD COMMENTlink modified 5.3 years ago by rbagnall1.4k • written 5.3 years ago by tommy.carstensen20

Is the low covg bam also PCR free?

ADD REPLYlink written 5.3 years ago by zam.iqbal.genome1.7k
0
gravatar for rbagnall
5.3 years ago by
rbagnall1.4k
Australia
rbagnall1.4k wrote:

There probably will be some differences.

The low coverage BAM may have not have a sufficient number of reads to reliably call some heterozygous sites, or may call some heterozygous sites as homozygous if the few reads at such a position are mostly from the alternate allele. The low coverage BAM may also have some PCR errors.

The high coverage BAM should not encounter these issues, and downsampling should not affect this, unless you have downsampled too much. Differences in variant calls between the high coverage BAM and the downsampled high coverage BAM could indicate excessive downsampling.

ADD COMMENTlink written 5.3 years ago by rbagnall1.4k
Please log in to add an answer.

Help
Access

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