At present we use only the ancestor's sequencing data when inferring the ancestral genotype in DescribeMutant. However, due to the nature of MA experiments, all the sequencing data tells us something about the probability of ancestral genotype (since most descendants will have the ancestral allele at most sites). Thus we shuold use the over-all genotype probabilities (i.e. anc_genotypes in TetMAProbability ) when determining the ancestral genotype.
NOTE: This only effects the from->to column in our output, not the mutation probabilities which already incorporate this information.
At present we use only the ancestor's sequencing data when inferring the ancestral genotype in
DescribeMutant
. However, due to the nature of MA experiments, all the sequencing data tells us something about the probability of ancestral genotype (since most descendants will have the ancestral allele at most sites). Thus we shuold use the over-all genotype probabilities (i.e.anc_genotypes
inTetMAProbability
) when determining the ancestral genotype.NOTE: This only effects the
from->to
column in our output, not the mutation probabilities which already incorporate this information.