Use of phase_block.sh

[zhidaobuzhidao]

Hello,
I would like to ask about the use of phase_block.sh.
I now have two haplotype assembly results, hap1.fasta and hap2.fasta. If I want to evaluate the assembled genome's swicherror, how do I use the phase_block.sh script?
In the script description, I need to input asm.fasta, hap1.meryl, hap2.meryl. What genome does asm.fasta refer to? Does it mean the genome I need to evaluate (hap1.fasta or hap2.fasta)? How do I get hap*.meryl? Is it the kmer result of using meryl to process hap1.fasta and hap2.fasta?
Look forward to a reply!

[arangrhie]

Hello, please read the wiki.
For phasing evaluation, Merqury relies on hapmers - inherited parental specific kmers, usually obtained from parental illumina/short reads. The hap1 hap2 meryl dbs needs to be built from these reads, independent from the assembly.

Best, Arang