The current GATK version is 3.8-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Get notifications!


You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Formatting tip!


Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Jump to another community
Download the latest Picard release at https://github.com/broadinstitute/picard/releases.
GATK version 4.beta.3 (i.e. the third beta release) is out. See the GATK4 beta page for download and details.

ComputeReadSpanCoverage

Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie
edited September 2012 in GenomeSTRiP Documentation

1. Introduction

The ComputeReadSpanCoverage walker traverses a set of BAM files to generate
genome-wide statistics.

The read span coverage is the count of bases in between two paired-end reads,
not counting the lengths of the reads themselves. For fixed-length reads of
length L with ungapped alignments, this would be InsertSize - 2*L. The
read span coverage is used as an estimate of the power for detecting
breakpoints using read pairs. This estimate assumes a model where the aligner
is unlikely to align a read to a breakpoint unless the breakpoint is close to
the end of the read.

Read pairs where the ends align to different sequences are never counted.

Read span coverage is computed and reported for each readgroup, but the output
is keyed by sample and library to allow easy roll up.

See also MergeReadSpanCoverage.

2. Inputs / Arguments

  • -I <bam-file> : The set of input BAM files.

  • -md <directory> : The metadata directory. Insert size histograms are
    loaded from the default isd.hist.bin file in this directory. This
    argument is also used to load a default list of excluded read groups.

  • -maxInsertSize <n> : Read pairs with an insert size greater than n are
    not counted in span coverage.

  • -maxInsertSizeStandardDeviations <sd> : Read pairs with an insert size
    greater than the median plus sd robust standard deviations are not counted
    in span coverage.

3. Outputs

  • -O <span-coverage-file> : Tab delimited output file (default is stdout).
Post edited by Geraldine_VdAuwera on
Sign In or Register to comment.