BAM write error in HTSeqGenie
1
0
Entering edit mode
zh9118 • 0
@zh9118-21668
Last seen 7 months ago
United States

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
HTSeqGenie BAM • 338 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 13 hours ago
United States

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.

Login before adding your answer.

Traffic: 643 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6