Reference
Count the number of reads
Genome Variant Analysis
java -jar GenomeAnalysisTK.jar -R reference.fasta -T CountReads -I input.bam [-L input.intervals]