added ira's script for extracting split reads from a pos sorted BAM