how to calculate aligned base count rna seq code example

Example: how to calculate aligned base count rna seq

bedtools genomecov -ibam input.bam -bga > output_aln.bam
awk '{s+=$4}END{print s}' output_aln.bam