HTSeq script to plot reads around tss