Input for cartesian plot about merfin

[ASLeonard]

I tried it anyway with with cut ... $illum.dump | paste $hifi.dump - | ..., so the axes may be flipped from the labels.

This was using the merged hap1 + hap2 fasta file with hifi and short reads, but the short reads had fairly lower coverage (~16x).

There is an approximate R of -0.03, but the top three values below accounted for ~ 61% of all points, and so probably bias that heavily.

3185301570      0.00    0.00
348440924       0.00    -1.00
143715045       0.00    1.00

It is interesting that the two axes are pretty heavily populated, but not the diagonal. I guess this may demonstrates that kmer bias for hifi is pretty independent of kmer bias for short reads?

[arangrhie]

Hi @ASLeonard , just saw this now. Sorry for the silence!

Yes, as far as I can tell, the k-mer bias was independent, so to speak.
The different error modes in HiFi and Illumina seem to be the cause of this;
we found homopolymer and microsatellite contraction in HiFi reads and the long-known GC biases in Illumina reads as shown here in T2T-CHM13.