implement the bam_csi_index parameter in atacseq pipeline

[BianyunY-ACRD]

Description of feature

I had success when running atacseq (-r 2.1.2) for samples with small genome. When I tried to run other samples with big genome and one of the chromosome size is larger than 500Mb. I set the "--bam_csi_index". However, I still got error message "....cannot be stored in a bai index. Try using a csi index". It seems the atacseq pipeline does not implement the bam_csi_index parameter. It would be greatly appreciated if this parameter can be implemented in atacseq pipeline.