HTSeq (RNA-seq counting)

htseq-count -r pos -t CDS -f bam S13.map.sorted.bam S13.gft > S13.count