I have reference mapped paired end illumina reads and called variants using BWA and Samtools respectively. The resulting vcf was treated to remove high...