an overlooked fix for multichrom gBGC

MesserLab · Jan 16, 2025 · eed96b8 · eed96b8
1 parent 4b22f8d
commit eed96b8
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/core/population.cpp b/core/population.cpp
@@ -4466,7 +4466,7 @@ void Population::DoHeteroduplexRepair(std::vector<slim_position_t> &p_heterodupl
 	Chromosome *chromosome = p_child_haplosome->AssociatedChromosome();
 	double gBGC_coeff_scaled = (chromosome->mismatch_repair_bias_ + 1.0) / 2.0;
 	bool repairs_biased = (species_.IsNucleotideBased() && (gBGC_coeff_scaled != 0.5));
-	NucleotideArray *ancestral_sequence = (repairs_biased ? species_.TheChromosome().AncestralSequence() : nullptr);
+	NucleotideArray *ancestral_sequence = (repairs_biased ? chromosome->AncestralSequence() : nullptr);
 	int heteroduplex_tract_count = (int)(p_heteroduplex.size() / 2);
 
 	if (heteroduplex_tract_count * 2 != (int)p_heteroduplex.size())