Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Possible issues for large genome (No index is available for this BAM file. and exceeding integer length in GTF file) #64

Open
JuiTse opened this issue Jan 17, 2025 · 0 comments

Comments

@JuiTse
Copy link

JuiTse commented Jan 17, 2025

Dear @mozack ,
Here I am processing the transcriptome for plant species with around 22g genome size. The bam file has already been sliced to per chromosome level rather than all-in-one bam file to match the acceptable bam file size. When processing local alignment by abra2, two issues I encountered.

  1. Using the samtools, the bai file cannot be created which arise from exceeding 536,870,912 positions in my bam file.
    As an alternative, I use -c function in samtool to create csi index for my bam file.

However, the following issue regarding "No index is available for this BAM file." still present.
Is there any way for abra2 to execute the csi index format? Or suggestions to solve this issue?

The log file is attached.
abra.log

  1. When I try on chromosome 1, the chromosome range of gtf file is too large, with the error message in the log file:
    NumberFormatException: For input string: "2147694110"

Is that possible to enlarge the acceptable "Integer data type" in java?
The log file is attached.
abra.log

I am looking forward to know more, thanks.

@JuiTse JuiTse changed the title Possible index issue for large genome (htsjdk.samtools.SAMException: No index is available for this BAM file.) Possible issues for large genome (No index is available for this BAM file. and exceeding integer length in GTF file) Jan 17, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant