search

marbl / merqury

k-mer based assembly evaluation
Other
272 stars 19 forks source link

hap_blot #34

Closed duhuipeng closed 3 years ago

duhuipeng commented 3 years ago

dear author I finished hap_blot this step, why is it like this, what caused it? image I found out in the log that this step was not wrong

arangrhie commented 3 years ago

This plot is showing how many hapmers are found in each contig/scaffold in the assemblies. Have you switched your father_purged and purged_mo assembly labelings, or the parental db naming?

duhuipeng commented 3 years ago

Dear author The code I run is as follows: ./merqury.sh Son.meryl Mother.hapmer.meryl Father.hapmer.meryl purged_mo.fasta father_purged.fasta test

But I didn't found error in Hap_blot.sh in log file,Looking forward to your advice,What's wrong with me Here is my log file for this step test.blob.log

arangrhie commented 3 years ago

Double check the assembly was properly named, as well as the kmer dbs. The Merqury plot is naming the axises following the input file. The blob plot is telling me your assembly is completely swapped, when the hapmers are correct.

arangrhie commented 3 years ago

Closing this issue as this is not a bug in merqury.

duhuipeng commented 3 years ago

Dear author image Why is the picture like this?Hope to get your advice