/usr/local/bin/bwa mem -R "@RG\tID:${bamfile}\tSM:${bamfile}\tPL:ILLUMINA\tLB:NexteraXT" -t4 -v2 ${h37}bwai.fasta ${dir}$file1 ${dir}$file2 | /usr/local/bin/samtools import ${h37}sami.fasta.fai - - | /usr/local/bin/samtools sort -@4 -m5G - ${dir}processed/$bamfile 
/usr/local/bin/samtools index ${dir}processed/${bamfile}.bam 
cd /usr/share/qualimap_v0.7.1; ./qualimap bamqc --java-mem-size=4G -bam ${dir}processed/${bamfile}.bam > /dev/null 
/usr/local/bin/samtools rmdup ${dir}processed/${bamfile}.bam ${dir}processed/${bamfile}_nodup.bam; mv ${dir}processed/${bamfile}_nodup.bam ${dir}processed/${bamfile}.bam; /usr/local/bin/samtools index ${dir}processed/${bamfile}.bam 
/usr/bin/java -Xmx24g -jar /usr/local/share/GenomeAnalysisTK.jar -T RealignerTargetCreator -R ${h37}sami.fasta -I ${dir}processed/${bamfile}.bam -o ${dir}processed/${bamfile}.bam.intervals -rf BadCigar -nt 12 
/usr/bin/java -Xmx4g -jar /usr/local/share/GenomeAnalysisTK.jar -T IndelRealigner -R ${h37}sami.fasta -I ${dir}processed/${bamfile}.bam -targetIntervals ${dir}processed/${bamfile}.bam.intervals -o ${dir}processed/${bamfile}_realigned.bam --defaultBaseQualities 12 
/usr/bin/java -Xmx4g -jar /usr/local/share/GenomeAnalysisTK.jar -T BaseRecalibrator -I ${dir}processed/${bamfile}_realigned.bam -R ${h37}sami.fasta --knownSites ${h37}Resi-List-16-12-2013.vcf -o ${dir}processed/${bamfile}_Resilist.grp -nct 8 -rf BadCigar 
rm ${dir}processed/${bamfile}.bam ${dir}processed/${bamfile}.bam.bai 
/usr/bin/java -Xmx4g -jar /usr/local/share/GenomeAnalysisTK.jar -T PrintReads -R ${h37}sami.fasta -I ${dir}processed/${bamfile}_realigned.bam -BQSR ${dir}processed/${bamfile}_Resilist.grp -o ${dir}processed/${bamfile}.bam -rf BadCigar -nct 8 
rm ${dir}processed/${bamfile}_realigned.bam ${dir}processed/${bamfile}_realigned.bai; cp ${dir}processed/${bamfile}.bai ${dir}processed/${bamfile}.bam.bai