python /oak/stanford/groups/akundaje/marinovg/code/SAMstats.py Oxy_TSSseq-pooled-STAR-2.5.3a/STAR-2.5.3a-Aligned.sortedByCoord.out.bam  SAMstats-Oxy_TSSseq-pooled -bam /oak/stanford/groups/akundaje/marinovg/genomes/Ciliates/Oxytricha_trifallax_asm_v1.1/GCA_000295675.1_oxytricha_asm_v1.1_genomic.chrom.sizes  samtools -paired