added code to output read counts at the first step, ie, bamtofastq step, and bamsort step, ie, the step after bwa alignment