fix equal short sequences being 1|2 for short alleles

wdecoster · Jul 2, 2024 · 1fa9fcc · 1fa9fcc
1 parent 25c43e7
commit 1fa9fcc
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/Cargo.toml b/Cargo.toml
@@ -1,6 +1,6 @@
 [package]
 name = "STRdust"
-version = "0.8.3"
+version = "0.8.4"
 edition = "2021"
 
 

diff --git a/src/vcf.rs b/src/vcf.rs
@@ -106,7 +106,7 @@ impl VCFRecord {
             "."
         } else if levenshtein(&allele2.seq, &repeat_ref_sequence) < repeat_ref_sequence.len() / 20 {
             "0"
-        } else if levenshtein(&allele2.seq, &allele1.seq) < allele1.seq.len() / 20 {
+        } else if allele2.seq == allele1.seq || levenshtein(&allele2.seq, &allele1.seq) < allele1.seq.len() / 20 {
             genotype1
         } else {
             next_alt