Question about --collapse

[Antonia-Chalka]

Question about --collapse

[Antonia-Chalka]

From a quick view at the code in the methods script, it seems the correlation has to be perfect, but there's also a mention of having a 'softer' mention so I'm not 100% sure 😅

[AdmiralenOla]

Thanks for your question, and sorry about the wait.

As you have already figured out, the genotypes need to be 100% correlated to be collapsed. You may also have seen from the code that I thought about using a softer threshold, but I have never gotten around to implementing that.

I'm also a bit uncertain how the distribution of the collapsed variant should be counted, i.e. should it be present in all isolates with either of the original variants? I'm uncertain how that would impact other assumptions that are made.

Another thing I'm not sure about is whether the collapsed genes should then go through subsequent rounds of correlation -> collapse. That is, when we collapse two genes into one, this will have a new distribution pattern, and there is a chance that this new pattern will fall within the correlation threshold of being collapsed with yet another gene.