Variant calling
Variant calling of both single nucleotide variations (SNVs) and small insertion/deletions (indels) was then performed on the processed alignment files using a combination of the SAMtools mpileup