search

marbl / merqury

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

NG block #81

Closed duhuipeng closed 1 year ago

duhuipeng commented 1 year ago

dear author Now I want to draw this plot separately image My input file format is as follows，but it appear error image Could you please provide the input file format Best

arangrhie commented 1 year ago

Hello, it should be the -b ....phased_block.sizes. See more details with

Rscript $MERQURY/plot/plot_block_N.R --help

For example:

block   mat.hapmer      6635314
block   mat.hapmer      6213435
block   mat.hapmer      5738677
block   mat.hapmer      5696453
block   mat.hapmer      5627365
block   mat.hapmer      5209980

The blocks are colored by col2; so extract the lines that matches what you want on col2. Not sure what you intend to measure though..? The blocks are identified on one assembly, so pulling only one type of block wouldn't make much sense on the N* scale.

duhuipeng commented 1 year ago

thanks , got it