Parse a BAM file and extract reads ovelapping given genomic coordinates