From 554eec9f95fb83d4f7ebeafe7b068820cc7144c8 Mon Sep 17 00:00:00 2001 From: youyupei Date: Sun, 30 Jun 2024 12:52:41 +1000 Subject: [PATCH] Refactor loop to handle larger sequence lengths in BAM records (to fix the memory issue when running demuxlet on long reads) --- hts_utils.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/hts_utils.cpp b/hts_utils.cpp index 4c0411d..9e2b0e0 100644 --- a/hts_utils.cpp +++ b/hts_utils.cpp @@ -175,7 +175,7 @@ void bam_get_seq_string(bam1_t *s, kstring_t *seq) { seq->l=0; uint8_t* sq = bam_get_seq(s); - for (uint16_t i = 0; i < bam_get_l_qseq(s); ++i) + for (uint32_t i = 0; i < bam_get_l_qseq(s); ++i) { kputc("=ACMGRSVTWYHKDBN"[bam_seqi(sq, i)], seq); }