this prints out read positions
TODO: get each gene's exons

because too many single exon genes have too low expression level,
it's seem necessary to analyze all the genes and their corresponding exons