Currently when I ran HTSeqGenie on RNA-Seq samples. The following error occurred:

[E::bam_write1] Positional data is too large for BAM format

This caused issues in writing bam files. Has anyone else seen this error before?

sessionInfo( )

R version 4.3.2 (2023-10-31)
Platform: x86_64-conda-linux-gnu (64-bit)
Running under: Red Hat Enterprise Linux Server 7.9 (Maipo)

Matrix products: default
BLAS/LAPACK: /project/pathology/Kapur_lab/Data_processing/Results_post_062022/Conda_env/py3.11-R4.3.2/lib/libopenblasp-r0.3.24.so;  LAPACK version 3.11.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

time zone: US/Central
tzcode source: system (glibc)

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
[1] compiler_4.3.2

You will have to use SAM format if you have a chromosome that is longer than 2Gb. Or you could hypothetically cut the chromosome(s) into pieces that are <2Gb, but that's up to you.