this script inputs a list of genome coordinates and calculates the average depth and breadth of coverage for each genome coordinate using a samtools mpileup file