Question: bam-readcount tool: output zero coverage bases
0
gravatar for user230613
7 months ago by
user230613250
Europe
user230613250 wrote:

I think the question title is self explanatory. Is it possible to output bases with zero coverage using bam-readcount tool. Something similar to samtools depth:

Usage: samtools depth [options] in1.bam [in2.bam [...]]
Options:
   -a   output all positions (including zero depth)
bam-readcount • 318 views
ADD COMMENTlink modified 7 months ago by Chris Miller19k • written 7 months ago by user230613250
1
gravatar for Santosh Anand
7 months ago by
Santosh Anand3.0k
Santosh Anand3.0k wrote:

is it not? The README says otherwise: "This program reports readcounts for each base at each position requested."!

ADD COMMENTlink written 7 months ago by Santosh Anand3.0k

You're right. My mistake :/

ADD REPLYlink written 7 months ago by user230613250
0
gravatar for Chris Miller
7 months ago by
Chris Miller19k
Washington University in St. Louis, MO
Chris Miller19k wrote:

You are correct that bam-readcount does not output zero-coverage bases. It requires a little post-processing to add those back in. (such as what's done here https://github.com/genome/genome/blob/master/lib/perl/Genome/Model/Tools/Analysis/Coverage/BamReadcount.pm)

ADD COMMENTlink written 7 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: 1408 users visited in the last hour