python /oak/stanford/groups/akundaje/marinovg/code/SAMstats.py HARMI-1-STAR-2.5.3a/STAR-2.5.3a-Aligned.sortedByCoord.out.bam  SAMstats-HARMI-1 -bam /oak/stanford/groups/akundaje/marinovg/genomes/Metazoa/Insects/Helicoverpa_armigera-Harm_1.0/GCA_002156985.1_Harm_1.0_genomic.chrom.sizes  samtools -paired 
python /oak/stanford/groups/akundaje/marinovg/code/SAMstats.py HARMI-2-STAR-2.5.3a/STAR-2.5.3a-Aligned.sortedByCoord.out.bam  SAMstats-HARMI-2 -bam /oak/stanford/groups/akundaje/marinovg/genomes/Metazoa/Insects/Helicoverpa_armigera-Harm_1.0/GCA_002156985.1_Harm_1.0_genomic.chrom.sizes  samtools -paired