QuickFix: Refactor markovprobability function to assert DNAAlphabet t…

…ype consistency in model and sequence
BioJulia · Dec 9, 2024 · d918766 · d918766
1 parent 15778eb
commit d918766
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/transitions.jl b/src/transitions.jl
@@ -259,8 +259,8 @@ function markovprobability(
     seq::NucleicSeqOrView{A};
     model::BioMarkovChain{A}=ECOLICDS,
     logscale::Bool=true
-) where {A<:NucleicAcidAlphabet}
-    @assert Alphabet(model) == Alphabet(seq) "Sequence and model state space are inconsistent."
+) where {A<:DNAAlphabet}
+    @assert eltype(model) == eltype(seq) "Sequence and model state space are inconsistent."
 
     init = logscale ? log2(model.inits[_dna_to_int(seq[1])]) : model.inits[_dna_to_int(seq[1])]
     probability = init