baozg / Potato_C88

Chromosome-scale and haplotype-resolved assmebly of tetraploid cultiavetd potato with polyploidy graph binning
23 stars 3 forks source link

Question for Figure 2A #3

Closed YvonneennovY closed 1 year ago

YvonneennovY commented 1 year ago

Hi Zhigui,

Thank you very much for your outstanding work, which has allowed me to learn a lot!

I am a bit confused about the analysis of Figure 2A in your article. You mentioned that you used the longest haplotype as the reference to detect the SNP/indel. So how you compared SNP/indel density among the six different haplotype pairs? Would it be possible to share the relevant scripts for Figure 2A? I would be immensely grateful.

Best wishes, Yvonne

baozg commented 1 year ago

If we assume the hap1 as a reference, then align the hap1 and hap2 to the hap1 and call variants. In this case, all variants in hap2 will be detected. If you compare hap2 and hap3, then only the heterozygous variants will be considered.

YvonneennovY commented 1 year ago

Thank you for your detailed response!