Try   HackMD

nf-core/genomeassembler assembly validation and curation discussion

Desired components:

  • K-mer copy number spectra profile - Merqury
  • Contamination profile - Blobtools
  • Gene space checks - Busco ( Dogma too?)
  • Assembly comparison and stats - Quast
  • Other accuracy metrics ( Inspector ? FRCBam ? )
  • Pairwise genome plots
  • Assembly graph plots ( Bandage )
  • Read coverage ( Qualimap, Blobtools )
  • Output genome browser tracks with useful information, like polished areas, missassembly signals, coverage, etc

Workflow diagram draft:

Assemblies

Quast

BUSCO

on-demand:Bakta

Merqury

PacBioHifi

Implementations:

  • Pairwise genome comparison workflow
  • Generate list of pairwise inputs:
    ​​workflow {
    ​​  Channel.of( 'A', 'B', 'C' )
    ​​      .toSortedList()
    ​​      .flatMap {
    ​​          it.subsequences()
    ​​              .findAll { it.size() == 2 }
    ​​      }.view()
    ​​  // Output is [B, C], [A, B], [A, C]
    ​​}