search

marbl / merqury

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

Is there a script to rerun the plot? #121

Closed Axze-rgb closed 3 months ago

Axze-rgb commented 3 months ago

Hello, my spectra CN plot cut a peak in the last third. I would like to remake it and give instructions for the x-axis. Is that possible? I mean is there a script I can rerun by changing a number? How do I have to redo the plot in R myself? Sorry for the many questions

Thanks

Axze-rgb commented 3 months ago

Yes there is, for the record : Rscript $MERQURY/plot/plot_spectra_cn.R -m 250 -n 2274025 -f H5A2.merqury.hairsplitter_final_assembly.spectra-cn.hist -o H5A2

This is inside a mamba install of merqury. Thanks

Axze-rgb commented 3 months ago

However, I didn't find how to remake the spectra asm?

arangrhie commented 3 months ago

Use the spectra-asm.hist instead for -f

Axze-rgb commented 3 months ago

Ah of course, I should have thought of it. Thanks!