## Load modules for feature counts
module load gcc/6.2.0 subread/1.6.2

## feature counts script
featureCounts -a /n/groups/churchman/djm46/Genomes/hg38/allPC_Genes.saf -g gene_id -s 2 -F SAF -o counts.txt sample.bam

##change parameters to adjust for specific samples