diff --git a/Cargo.lock b/Cargo.lock index 7c883c5..5875344 100644 --- a/Cargo.lock +++ b/Cargo.lock @@ -1699,7 +1699,7 @@ checksum = "1f3ccbac311fea05f86f61904b462b55fb3df8837a366dfc601a0161d0532f20" [[package]] name = "trgt" -version = "1.4.0" +version = "1.4.1" dependencies = [ "arrayvec", "bio", diff --git a/Cargo.toml b/Cargo.toml index 3ee5753..552639c 100644 --- a/Cargo.toml +++ b/Cargo.toml @@ -1,6 +1,6 @@ [package] name = "trgt" -version = "1.4.0" +version = "1.4.1" edition = "2021" build = "build.rs" diff --git a/README.md b/README.md index 1e8c6c8..902f212 100644 --- a/README.md +++ b/README.md @@ -146,6 +146,8 @@ tandem repeats at genome scale. 2024](https://www.nature.com/articles/s41587-023 - Parameters appropriate for targeted sequencing can now be set with `--preset targeted` option - Waterfall plots no longer panic when there are no reads in a locus - Algorithmic changes to `--genotyper cluster` allow fewer reads to be assigned to an allele; this may result in minor changes to consensus sequence and read assignment +- 1.4.1 + - Bug fix: corrected the type of the rq tag in BAM output ### DISCLAIMER diff --git a/src/trgt/writers/write_bam.rs b/src/trgt/writers/write_bam.rs index b253d97..9663b53 100644 --- a/src/trgt/writers/write_bam.rs +++ b/src/trgt/writers/write_bam.rs @@ -113,8 +113,9 @@ impl BamWriter { let tr_tag = Aux::String(&locus.id); rec.push_aux(b"TR", tr_tag).unwrap(); - rec.push_aux(b"rq", Aux::Double(read.read_qual.unwrap_or(-1.0))) + rec.push_aux(b"rq", Aux::Float(read.read_qual.unwrap_or(-1.0) as f32)) .unwrap(); + if let Some(meth) = &read.meth { let mc_tag = Aux::ArrayU8(meth.into()); rec.push_aux(b"MC", mc_tag).unwrap();