(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.