added Rsamtools script