Question: bam-readcount values are higher than Mutect results
0
gravatar for sey
10 months ago by
sey0
sey0 wrote:

I use Mutect to find somatic mutations in multiple tumor samples (to be used in clonality analysis). I don't get the readcounts if the caller did not identify a somatic mutation in a sample. I'm using bam-readcount to fill in the empty values after merging Mutect results from different samples. Where I have values from both tools, the readcounts are lower in the Mutect results. Is there a way to use bam-readcount in a way to substitute for missing Mutect values?

bam-readcount • 323 views
ADD COMMENTlink modified 10 months ago by Chris Miller19k • written 10 months ago by sey0
0
gravatar for Chris Miller
10 months ago by
Chris Miller19k
Washington University in St. Louis, MO
Chris Miller19k wrote:

Many callers do some filtering of reads first (by base quality, by mapping quality, maybe even local realignment) and report only counts from reads that pass. bam-readcount, by default, outputs "raw" counts. You may or may not be able to mimic those Mutect settings with bam-readcount parameters.

This problem is compounded when merging the output from different variant callers, and there is no "perfect" solution.

ADD COMMENTlink written 10 months ago by Chris Miller19k
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: 1805 users visited in the last hour