A similar routine as for RNA-seq data required for BWA alignment for the DNA-seq path (when input is BAM file)