- custom gatk workflow adopted from CIDR
- Modify to suit the setting at local servers
- Install bwa, picard, and gatk in $HOME/bin
- Download gatk_bundle files from GATK ftp sites in $GATK_RESOURCE_PATH
- this includes reference genome (hg19), known indels, dbsnp138
- set up bwa index
- put fastq files in $PROJECT_FOLDER/$FASTQ_DIR
- This script will create a subscript for each process in $SAMPLE/$SCRIPT
- Change the number of threads to reflect your CPUs
- Change the INPUT DATA below to match with your input