Skip to content

Conversation

@joshfactorial
Copy link
Collaborator

No description provided.

@joshfactorial joshfactorial changed the base branch from main to develop November 7, 2025 21:38
@joshfactorial
Copy link
Collaborator Author

I'm tweaking when we use bgzip, which is quite slow. But I stumbled upon another issue, which is that the bam alignments are all showing up on the first 500,000 bases in the bam (when viewing in the IGV genome browser), so there must be an index off somwhere. I checkd chr1 of rice and the vcf index goes up to 44.3 million, which is roughly the size of chr1, so I think we are good there. Something still in the coordinates of the bam output.

@joshfactorial joshfactorial marked this pull request as ready for review November 9, 2025 07:17
@joshfactorial
Copy link
Collaborator Author

I made a small change to how the files were written, to try to speed up intermediate file creation. While testing that, I found an error in the way bams were indexing reads, causing only 500,000 bases worth of reads to get written. Now it spans the chromosome, but the coverage looks off.

@joshfactorial joshfactorial merged commit c3eb220 into develop Nov 9, 2025
@joshfactorial joshfactorial deleted the check-slow-splits branch November 9, 2025 07:19
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

Successfully merging this pull request may close these issues.

2 participants