First set up your Zorn/Bascet work directory as before.
Assembly can then performed using SKESA, via the MapCell system:
### Assemble all genomes
BascetMapCell(
bascetRoot,
withfunction = "_skesa",
inputName = "filtered",
outputName = "skesa"
)
This produce contigs for each cell, all grouped together in “skesa”. Note that both skesa and spades only assembles cells for which there are sufficient number of reads.
#TODO extract an individual genome
#TODO QUAST for QC etc