How to Create a snippet of reads corresponding to a genomic interval