Next changeset 1:3169a38c2656 (2021-06-26) |
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/genomescope commit aa87b7b1713b749328c5a710f32631aab2acaa3a" |
added:
genomescope.xml test-data/genomescope-in1.tab test-data/genomescope-out1-1.png test-data/genomescope-out1-1.tab test-data/genomescope-out1-1.txt test-data/genomescope-out1-2.png test-data/genomescope-out1-2.txt test-data/genomescope-out1-3.png test-data/genomescope-out1-3.txt test-data/genomescope-out1-4.png |
b |
diff -r 000000000000 -r b2f674562a18 genomescope.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/genomescope.xml Fri Apr 30 20:21:25 2021 +0000 |
[ |
@@ -0,0 +1,123 @@ +<tool id="genomescope" name="GenomeScope" version="@VERSION@" profile="20.01"> + <description>Analyze unassembled short reads</description> + <macros> + <token name="@VERSION@">2.0</token> + </macros> + <requirements> + <requirement type="package" version="@VERSION@">genomescope2</requirement> + </requirements> + <version_command>genomescope2 --version</version_command> + <command detect_errors="exit_code"><![CDATA[ + genomescope2 + --input '$input' + --output . + --kmer_length $kmer_length + $no_unique_sequence $testing $trace_flag + #if $ploidy: + --ploidy $ploidy + #end if + #if $lambda: + --lambda $lambda + #end if + #if $max_kmercov: + --max_kmercov $max_kmercov + #end if + #if $topology: + --topology $topology + #end if + #if $initial_repetitiveness: + --initial_repetitiveness $initial_repetitiveness + #end if + #if $initial_heterozygosities: + --initial_heterozygosities $initial_heterozygosities + #end if + #if $transform_exp: + --transform_exp $transform_exp + #end if + #if $true_params: + --true_params $true_params + #end if + #if $num_rounds: + --num_rounds $num_rounds + #end if + ]]> + </command> + <inputs> + <param argument="--input" type="data" format="tabular" label="Input histogram file" help="This file is a two column tabular file for example generated with the histo function of Jellyfish."/> + <param name="model_output" type="boolean" label="Add the model parameters to your history"/> + <param name="summary_output" type="boolean" label="Output a summary of the analysis"/> + <param name="progress_output" type="boolean" label="Additional information for each optimization round"/> + <param argument="--ploidy" type="integer" optional="true" label="Ploidy for model to use" help="Default: 2"/> + <param argument="--kmer_length" type="integer" value="21" optional="false" label="K-mer length used to calculate k-mer spectra"/> + <param argument="--lambda" type="integer" optional="true" label="Optional initial kmercov estimate for model to use"/> + <param argument="--max_kmercov" type="integer" optional="true" label="Optional maximum k-mer coverage threshold" help="K-mers with coverage greater than max_kmercov are ignored by the model"/> + <param argument="--no_unique_sequence" type="boolean" truevalue="--no_unique_sequence" falsevalue="" label="Turn off yellow unique sequence line in plots"/> + <param argument="--topology" type="integer" optional="true" label="Flag for topology for model to use"/> + <param argument="--initial_repetitiveness" type="integer" optional="true" label="Initial value for repetitiveness"/> + <param argument="--initial_heterozygosities" type="integer" optional="true" label="Initial values for nucleotide heterozygosity rates"/> + <param argument="--transform_exp" type="integer" optional="true" label="Parameter for the exponent when fitting a transformed (x**transform_exp*y vs. x) k-mer histogram" help="Default: 1"/> + <param argument="--testing" type="boolean" truevalue="--testing" falsevalue="" label="Create testing.tsv file with model parameters"/> + <param argument="--true_params" type="integer" optional="true" label="Flag to state true simulated parameters for testing mode"/> + <param argument="--trace_flag" type="boolean" truevalue="--trace_flag" falsevalue="" label="Turn on printing of iteration progress of nlsLM function"/> + <param argument="--num_rounds" type="integer" min="1" optional="true" label="Number of optimization rounds"/> + </inputs> + <outputs> + <data name="linear_plot" format="png" from_work_dir="linear_plot.png" label="${tool.name} on ${on_string} Linear plot"/> + <data name="log_plot" format="png" from_work_dir="log_plot.png" label="${tool.name} on ${on_string} Log plot"/> + <data name="transformed_linear_plot" format="png" from_work_dir="transformed_linear_plot.png" label="${tool.name} on ${on_string} Transformed linear plot"/> + <data name="transformed_log_plot" format="png" from_work_dir="transformed_log_plot.png" label="${tool.name} on ${on_string} Transformed log plot"/> + <data name="model" format="txt" from_work_dir="model.txt" label="${tool.name} on ${on_string} Model"> + <filter>model_output</filter> + </data> + <data name="summary" format="txt" from_work_dir="summary.txt" label="${tool.name} on ${on_string} Summary"> + <filter>summary_output</filter> + </data> + <data name="progress" format="txt" from_work_dir="progress.txt" label="${tool.name} on ${on_string} Progress"> + <filter>progress_output</filter> + </data> + <data name="model_params" format="tabular" from_work_dir="SIMULATED_testing.tsv" label="${tool.name} on ${on_string} Model parameters"> + <filter>testing</filter> + </data> + </outputs> + <tests> + <test expect_num_outputs="8"> + <param name="input" value="genomescope-in1.tab"/> + <param name="kmer_length" value="21"/> + <param name="testing" value="true"/> + <param name="model_output" value="true"/> + <param name="summary_output" value="true"/> + <param name="progress_output" value="true"/> + <output name="linear_plot" file="genomescope-out1-1.png" ftype="png"/> + <output name="log_plot" file="genomescope-out1-2.png" ftype="png"/> + <output name="transformed_linear_plot" file="genomescope-out1-3.png" ftype="png"/> + <output name="transformed_log_plot" file="genomescope-out1-4.png" ftype="png" compare="sim_size"/> + <output name="model" file="genomescope-out1-1.txt" ftype="txt" lines_diff="2"/> + <output name="summary" file="genomescope-out1-2.txt" ftype="txt" lines_diff="2"/> + <output name="progress" file="genomescope-out1-3.txt" ftype="txt" lines_diff="2"/> + <output name="testing" file="genomescope-out1-1.tab" ftype="tabular"/> + </test> + </tests> + <help><![CDATA[ + +GenomeScope 2.0: Reference-free profiling of polyploid genomes +============================================================== + +GenomeScope 2.0 applies classical insights from combinatorial theory to establish +a detailed mathematical model of how k-mer frequencies will be distributed in +heterozygous and polyploid genomes. GenomeScope 2.0 employs a polyploid-aware +mixture model that, within seconds, accurately infers genome properties from +unassembled sequencing data. GenomeScope 2.0 uses the k-mer count distribution, +e.g. from KMC or Jellyfish, and produces a report and several informative plots +describing the genome properties. We validate the approach on simulated polyploid +data created using a generative model with parameters for genome size, heterozygosity, +repetitiveness, ploidy, and sequencing coverage, and find GenomeScope 2.0 retains +accuracy across a broad range of realistic and extreme parameter values. We also +validate GenomeScope 2.0 by analyzing genuine sequence data from 11 diverse +polyploid genomes with known genome characteristics. + +]]></help> + <citations> + <citation type="doi">10.1093/bioinformatics/btx153</citation> + <citation type="doi">10.1038/s41467-020-14998-3</citation> + </citations> +</tool> |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-in1.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/genomescope-in1.tab Fri Apr 30 20:21:25 2021 +0000 |
b |
b'@@ -0,0 +1,24742 @@\n+1 398558798\n+2 20510187\n+3 4448867\n+4 2021116\n+5 1095317\n+6 684526\n+7 488031\n+8 403773\n+9 396767\n+10 453765\n+11 579419\n+12 769072\n+13 1021830\n+14 1333141\n+15 1682372\n+16 2043149\n+17 2380557\n+18 2682538\n+19 2915377\n+20 3060490\n+21 3113837\n+22 3076739\n+23 2962785\n+24 2783624\n+25 2576341\n+26 2362401\n+27 2155614\n+28 1985594\n+29 1880343\n+30 1834668\n+31 1861054\n+32 1949634\n+33 2102849\n+34 2301698\n+35 2551528\n+36 2824079\n+37 3101356\n+38 3378532\n+39 3628757\n+40 3845823\n+41 4009857\n+42 4130915\n+43 4179976\n+44 4165381\n+45 4080784\n+46 3945785\n+47 3753722\n+48 3534244\n+49 3273383\n+50 2998681\n+51 2718765\n+52 2429130\n+53 2145537\n+54 1868844\n+55 1613769\n+56 1381725\n+57 1171762\n+58 987274\n+59 825001\n+60 685604\n+61 569315\n+62 469637\n+63 387661\n+64 320226\n+65 267381\n+66 224627\n+67 190073\n+68 161812\n+69 139410\n+70 121982\n+71 110522\n+72 101679\n+73 93848\n+74 88622\n+75 83801\n+76 80197\n+77 76740\n+78 74593\n+79 72817\n+80 70641\n+81 69342\n+82 67272\n+83 66487\n+84 65950\n+85 64597\n+86 63045\n+87 62202\n+88 60132\n+89 58933\n+90 57572\n+91 55739\n+92 53404\n+93 51857\n+94 49666\n+95 47076\n+96 45242\n+97 43492\n+98 41569\n+99 39848\n+100 37681\n+101 35564\n+102 33867\n+103 31977\n+104 30786\n+105 29551\n+106 28330\n+107 27059\n+108 26038\n+109 24894\n+110 23574\n+111 23007\n+112 21998\n+113 21083\n+114 20562\n+115 19833\n+116 19124\n+117 18571\n+118 17887\n+119 17682\n+120 17185\n+121 16790\n+122 16665\n+123 16041\n+124 15647\n+125 15767\n+126 15179\n+127 14966\n+128 14579\n+129 14605\n+130 14322\n+131 13781\n+132 13673\n+133 13472\n+134 13369\n+135 12757\n+136 12691\n+137 12573\n+138 12329\n+139 12097\n+140 11883\n+141 11833\n+142 11447\n+143 11496\n+144 11158\n+145 11103\n+146 10791\n+147 10815\n+148 10596\n+149 10367\n+150 9809\n+151 9892\n+152 9646\n+153 9486\n+154 9342\n+155 9302\n+156 9055\n+157 8708\n+158 8739\n+159 8430\n+160 8441\n+161 8231\n+162 8219\n+163 7992\n+164 7816\n+165 7679\n+166 7556\n+167 7462\n+168 7183\n+169 7141\n+170 7023\n+171 6927\n+172 6772\n+173 6786\n+174 6695\n+175 6719\n+176 6505\n+177 6431\n+178 6440\n+179 6210\n+180 6013\n+181 6120\n+182 5834\n+183 5808\n+184 5721\n+185 5677\n+186 5814\n+187 5797\n+188 5783\n+189 5586\n+190 5644\n+191 5300\n+192 5531\n+193 5386\n+194 5559\n+195 5405\n+196 5448\n+197 5336\n+198 5215\n+199 5192\n+200 5018\n+201 5159\n+202 5335\n+203 5100\n+204 5245\n+205 5283\n+206 5261\n+207 5056\n+208 4913\n+209 5021\n+210 5038\n+211 5044\n+212 5189\n+213 5308\n+214 5290\n+215 5289\n+216 5269\n+217 5277\n+218 5444\n+219 5507\n+220 5272\n+221 5271\n+222 5268\n+223 5165\n+224 5362\n+225 5269\n+226 5412\n+227 5454\n+228 5455\n+229 5485\n+230 5477\n+231 5614\n+232 5416\n+233 5529\n+234 5522\n+235 5422\n+236 5610\n+237 5488\n+238 5462\n+239 5657\n+240 5827\n+241 5699\n+242 5557\n+243 5774\n+244 5703\n+245 5570\n+246 5742\n+247 5653\n+248 5704\n+249 5790\n+250 5683\n+251 5812\n+252 5950\n+253 5998\n+254 5964\n+255 5957\n+256 5976\n+257 5947\n+258 5781\n+259 6023\n+260 5726\n+261 5625\n+262 5498\n+263 5623\n+264 5702\n+265 5692\n+266 5569\n+267 5689\n+268 5643\n+269 5482\n+270 5526\n+271 5584\n+272 5506\n+273 5393\n+274 5291\n+275 5387\n+276 5309\n+277 5198\n+278 5184\n+279 5228\n+280 5066\n+281 4917\n+282 4671\n+283 4570\n+284 4574\n+285 4544\n+286 4505\n+287 4485\n+288 4356\n+289 4444\n+290 4278\n+291 4321\n+292 4166\n+293 4069\n+294 3876\n+295 3962\n+296 3836\n+297 3625\n+298 3660\n+299 3669\n+300 3460\n+301 3387\n+302 3325\n+303 3285\n+304 3239\n+305 3134\n+306 2922\n+307 2903\n+308 2642\n+309 2625\n+310 2446\n+311 2504\n+312 2362\n+313 2349\n+314 2382\n+315 2381\n+316 2252\n+317 2195\n+318 2075\n+319 1972\n+320 1989\n+321 1924\n+322 1825\n+323 1894\n+324 1828\n+325 1693\n+326 1724\n+327 1612\n+328 1650\n+329 1623\n+330 1612\n+331 1660\n+332 1567\n+333 1589\n+334 1529\n+335 1432\n+336 1490\n+337 1481\n+338 1503\n+339 1403\n+340 1392\n+341 1376\n+342 1376\n+343 1352\n+344 1284\n+345 1207\n+346 1309\n+347 1247\n+348 1182\n+349 1224\n+350 1150\n+351 1110\n+352 1140\n+353 1108\n+354 1204\n+355 1126\n+356 1086\n+357 1115\n+358 1194\n+359 1144\n+360 1110\n+361 1139\n+362 1062\n+363 1069\n+364 1080\n+365 1055\n+366 1018\n+367 1053\n+368 1044\n+369 1009\n+370 1034\n+371 1001\n+372 1000\n+373 1021\n+374 1052\n+375 927\n+376 983\n+377 901\n+378 908\n+379 945\n+380 965\n+381 912\n+382 917\n+'..b'726 1\n+118815 1\n+118886 1\n+118984 1\n+119133 1\n+119285 1\n+119416 1\n+119544 1\n+119611 1\n+120137 1\n+120179 1\n+120233 1\n+120254 1\n+120885 1\n+121817 1\n+121999 1\n+122844 1\n+123248 1\n+123770 1\n+124157 1\n+124349 1\n+124788 1\n+124860 1\n+124899 1\n+124919 1\n+125007 1\n+125198 1\n+125439 1\n+125595 1\n+125602 1\n+125619 1\n+125636 1\n+125706 1\n+126188 1\n+126515 1\n+126892 1\n+127439 1\n+127542 1\n+127643 1\n+127976 1\n+128027 1\n+128282 1\n+128400 1\n+128415 1\n+128620 1\n+130382 1\n+130625 1\n+130736 1\n+130798 1\n+131045 1\n+131060 1\n+131067 1\n+131361 1\n+131627 1\n+131738 1\n+132072 1\n+132242 1\n+133089 1\n+133434 1\n+133549 1\n+133579 1\n+133848 1\n+134025 1\n+134614 1\n+134727 1\n+134812 1\n+135203 1\n+135380 1\n+135947 1\n+137266 1\n+138112 1\n+138617 1\n+140619 1\n+141092 1\n+141693 1\n+142247 1\n+142544 1\n+142644 1\n+143553 1\n+143920 1\n+144051 1\n+144066 1\n+144368 1\n+146225 1\n+147719 1\n+147836 1\n+147924 1\n+148878 1\n+150351 1\n+150928 1\n+151071 1\n+152152 1\n+152460 1\n+153039 1\n+153527 1\n+153545 1\n+153841 1\n+154029 1\n+154358 1\n+154515 1\n+154679 1\n+155068 1\n+155624 1\n+156347 1\n+157436 1\n+157442 1\n+158094 1\n+158228 1\n+158490 1\n+158606 1\n+159165 1\n+159772 1\n+160040 1\n+160057 1\n+163239 1\n+164600 1\n+165103 1\n+165762 1\n+166301 1\n+166366 1\n+166941 1\n+167042 1\n+169241 1\n+169854 1\n+170284 1\n+170327 1\n+170677 1\n+171393 1\n+180619 1\n+181990 1\n+183353 1\n+183368 1\n+183496 1\n+186025 1\n+188096 1\n+188128 1\n+188239 1\n+189251 1\n+189684 1\n+189746 1\n+192281 1\n+193058 1\n+193406 1\n+194460 1\n+194486 1\n+194587 1\n+194932 1\n+195372 1\n+195383 1\n+195522 1\n+195774 1\n+195871 1\n+197251 1\n+197273 1\n+197466 1\n+197994 1\n+198209 1\n+198797 1\n+199073 1\n+199444 1\n+200758 1\n+202059 1\n+202452 1\n+202953 1\n+203095 1\n+204391 1\n+204419 1\n+204834 1\n+205284 1\n+205623 1\n+206031 1\n+206262 1\n+207238 1\n+207941 1\n+208922 1\n+208992 1\n+209957 1\n+211230 1\n+211391 1\n+212159 1\n+212914 1\n+213722 1\n+214997 1\n+215255 1\n+215440 1\n+215576 1\n+215687 1\n+215814 1\n+215938 1\n+215997 1\n+216139 1\n+216643 1\n+216926 1\n+217096 1\n+217238 1\n+217361 1\n+217808 1\n+218253 1\n+218437 1\n+218618 1\n+218661 1\n+219661 1\n+220448 1\n+221452 1\n+224504 1\n+224946 1\n+226113 1\n+226341 1\n+228254 1\n+228599 1\n+228901 1\n+229043 1\n+231018 1\n+232106 1\n+233370 1\n+233957 1\n+235102 1\n+237155 1\n+237317 1\n+238775 1\n+240168 1\n+241109 1\n+260899 1\n+261929 1\n+264767 1\n+265764 1\n+270290 1\n+275221 1\n+282111 1\n+288233 1\n+289131 1\n+289221 1\n+294131 1\n+294225 1\n+295541 1\n+295774 1\n+299266 1\n+299622 1\n+299964 1\n+312382 1\n+313090 1\n+326306 1\n+326884 1\n+327519 1\n+327529 1\n+328801 1\n+329719 1\n+329747 1\n+329777 1\n+330756 1\n+331458 1\n+334788 1\n+335587 1\n+336667 1\n+337799 1\n+337971 1\n+339038 1\n+339941 1\n+340370 1\n+340703 1\n+341448 1\n+341991 1\n+343677 1\n+347797 1\n+349650 1\n+349724 1\n+350764 1\n+351157 1\n+351374 1\n+351401 1\n+352249 1\n+352392 1\n+352507 1\n+352665 1\n+352923 1\n+353088 1\n+353921 1\n+354365 1\n+354424 1\n+354735 1\n+355199 1\n+355707 1\n+356901 1\n+359866 1\n+360393 1\n+360779 1\n+361127 1\n+361521 1\n+364242 1\n+365372 1\n+366298 1\n+369013 1\n+369574 1\n+369678 1\n+370800 1\n+372553 1\n+374601 1\n+379438 1\n+381815 1\n+382199 1\n+383448 1\n+384752 1\n+385882 1\n+387084 1\n+387183 1\n+387320 1\n+389460 1\n+389566 1\n+393259 1\n+395184 1\n+395558 1\n+396392 1\n+401927 1\n+403274 1\n+404960 1\n+407700 1\n+413321 1\n+413762 1\n+437868 1\n+451830 1\n+456464 1\n+457637 1\n+460164 1\n+463355 1\n+464228 1\n+466550 1\n+469431 1\n+470274 1\n+470415 1\n+470765 1\n+474703 1\n+475576 1\n+475931 1\n+477419 1\n+481103 1\n+481842 1\n+484647 1\n+485159 1\n+488865 1\n+489766 1\n+492807 1\n+496142 1\n+496746 1\n+497670 1\n+502878 1\n+507149 1\n+508470 1\n+508487 1\n+508520 1\n+508765 1\n+511694 1\n+512041 1\n+512595 1\n+512627 1\n+513872 1\n+514358 1\n+514425 1\n+514473 1\n+516732 1\n+516769 1\n+517897 1\n+521053 1\n+521529 1\n+523037 1\n+541291 1\n+547000 1\n+565310 1\n+565466 1\n+577583 1\n+581379 1\n+598225 1\n+600139 1\n+601177 1\n+605676 1\n+608386 1\n+610384 1\n+613154 1\n+616398 1\n+616774 1\n+618082 1\n+621162 1\n+621693 1\n+624929 1\n+627300 1\n+636335 1\n+637144 1\n+663780 1\n+730833 1\n+737802 1\n+738025 1\n+738828 1\n+742502 1\n+742775 1\n+745763 1\n+748008 1\n+965147 1\n+967487 1\n+985751 1\n+987948 1\n+1000001 134\n' |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-1.png |
b |
Binary file test-data/genomescope-out1-1.png has changed |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-1.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/genomescope-out1-1.tab Fri Apr 30 20:21:25 2021 +0000 |
b |
@@ -0,0 +1,1 @@ +0.0286284442637493 0.0103961720120016 22.1589612104347 0.458078023574561 151946613.878237 0 |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/genomescope-out1-1.txt Fri Apr 30 20:21:25 2021 +0000 |
b |
@@ -0,0 +1,19 @@ + +Formula: y_transform ~ x^transform_exp * length * predict2_0(r1, k, d, + kmercov, bias, x) + +Parameters: + Estimate Std. Error t value Pr(>|t|) +d 2.863e-02 4.312e-04 66.39 <2e-16 *** +r1 1.040e-02 3.291e-05 315.91 <2e-16 *** +kmercov 2.216e+01 4.424e-03 5008.85 <2e-16 *** +bias 4.581e-01 2.786e-03 164.42 <2e-16 *** +length 1.065e+08 9.361e+04 1137.81 <2e-16 *** +--- +Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 + +Residual standard error: 489000 on 1995 degrees of freedom + +Number of iterations to convergence: 5 +Achieved convergence tolerance: 1.49e-08 + |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-2.png |
b |
Binary file test-data/genomescope-out1-2.png has changed |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-2.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/genomescope-out1-2.txt Fri Apr 30 20:21:25 2021 +0000 |
b |
@@ -0,0 +1,15 @@ +GenomeScope version 2.0 +input file = /tmp/tmp32yh4464/files/1/e/2/dataset_1e2ac011-2bc8-45cf-971a-75abc3cda6d6.dat +output directory = . +p = 2 +k = 21 +TESTING set to TRUE + +property min max +Homozygous (aa) 98.9538% 98.967% +Heterozygous (ab) 1.03304% 1.0462% +Genome Haploid Length 151,885,967 bp 152,007,309 bp +Genome Repeat Length 48,469,025 bp 48,507,747 bp +Genome Unique Length 103,416,942 bp 103,499,562 bp +Model Fit 69.6762% 98.569% +Read Error Rate 0.349793% 0.349793% |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-3.png |
b |
Binary file test-data/genomescope-out1-3.png has changed |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-3.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/genomescope-out1-3.txt Fri Apr 30 20:21:25 2021 +0000 |
b |
@@ -0,0 +1,6 @@ +starting +round 0 trimming to 8 trying 2p peak model... converged. score: 759966810268782 +round 1 trimming to 13 trying 2p peak model... converged. score: 758881495868533 +round 2 trimming to 18 trying 2p peak model... converged. score: 752379509919454 +round 3 trimming to 23 trying 2p peak model... converged. score: 747542116457350 +done |
b |
diff -r 000000000000 -r b2f674562a18 test-data/genomescope-out1-4.png |
b |
Binary file test-data/genomescope-out1-4.png has changed |