change code to use tmap samtools