This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!
(howto) Compress read data with ReduceReads - DEPRECATED
Please note that this article refers to a method that is no longer recommended as part of the Best Practices!
Compress the read data in order to minimize file sizes, which facilitates massively multisample processing.
- Compress your sequence data
1. Compress your sequence data
Run the following GATK command:
java -jar GenomeAnalysisTK.jar \ -T ReduceReads \ -R reference.fa \ -I recal_reads.bam \ -L 20 \ -o reduced_reads.bam
This creates a file called
reduced_reads.bam containing only the sequence information that is essential for calling variants.
Note that ReduceReads is not meant to be run on multiple samples at once. If you plan on merging your sample bam files, you should run ReduceReads on individual samples before doing so.