Walks over the input data set, calculating the number of reads seen for diagnostic purposes.
Can also count the number of reads matching a given criterion using read filters (see the --read-filter command line argument). Simplest example of a read-backed analysis.
One or more BAM files.
Number of reads seen.
java -Xmx2g -jar GenomeAnalysisTK.jar \ -R ref.fasta \ -T CountReads \ -I input.bam \ [-L input.intervals]
This Read Filter is automatically applied to the data by the Engine before processing by CountReads.
This tool can be run in multi-threaded mode using this option.
This tool does not apply any downsampling by default.
The arguments described in the entries below can be supplied to this tool to modify its behavior. For example, the -L argument directs the GATK engine restricts processing to specific genomic intervals (this is an Engine capability and is therefore available to all GATK walkers).
GATK version 3.0-0-g6bad1c6 built at 2014/03/06 06:38:04. GTD: NA