Get exon regions from ref_annot.bed ...
Load BAM file ...  Done
Total read pairs  used 4095