Address of the bookmark: https://bioinfologics.github.io/post/2018/09/17/k-mer-counting-part-i-introduction/

https://academic.oup.com/nar/article/47/11/e63/5377471

Address of the bookmark: https://github.com/linzhi2013/MitoZ

https://www.bioconductor.org/packages/devel/bioc/vignettes/maftools/inst/doc/maftools.html

Address of the bookmark: https://github.com/PoisonAlien/maftools

You can try this here, or try it later on your own data.

Get data

We will use the same Illumina data as we used above:

• illumina_R1.fastq.gz: the Illumina forward reads
• illumina_R2.fastq.gz: the Illumina reverse reads

Assemble

Run Spades:

spades.py -1 illumina_R1.fastq.gz -2 illumina_R2.fastq.gz --careful --cov-cutoff auto -o spades_assembly_all_illumina

• -1 is input file of forward reads
• -2 is input file of reverse reads
• --careful minimizes mismatches and short indels
• --cov-cutoff auto computes the coverage threshold (rather than the default setting, “off”)
• -o is the output directory

Results

Move into the output directory and look at the contigs:

infoseq contigs.fasta

Smudgeplots are computed from raw or even better from trimmed reads and show the haplotype structure using heterozygous kmer pairs. For example:

Address of the bookmark: https://github.com/KamilSJaron/smudgeplot

Address of the bookmark: https://www.ncbi.nlm.nih.gov/nuccore/ON563414

MitoHiFi was first developed to assemble the mitogenomes for a wide range of species in the Darwin Tree of Life Project (DToL)

https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-023-05385-y 

Address of the bookmark: https://github.com/marcelauliano/MitoHiFi

• AssemblyMapper : reference-guided genome assembly
• GapFiller : long-reads based gap filling
• TeloExplorer : telomere identification
• CentroMiner : centromere candidate prediction

https://academic.oup.com/hr/article/10/8/uhad127/7197191?login=false 

Address of the bookmark: http://www.atcgn.com:8080/quarTeT/home.html

nf-core/mag v5 (https://github.com/nf-core/mag) addresses this challenge by extending its reproducible Nextflow-based workflow for modern genome-resolved metagenomics.

A key addition is support for long-read-only metagenomic assembly and bin refinement, enabling long-read data to be processed through an integrated workflow. The release also introduces five additional binning tools, providing complementary strategies for recovering genomes from complex microbial communities.

The workflow has also expanded beyond conventional bacterial and archaeal MAGs, with improved classification of viruses and eukaryotes, together with enhanced genome-quality assessment.

Conceptually, the workflow brings together:

Reads → QC → Assembly → Binning → Bin refinement → Taxonomic classification → MAG quality assessment

The real strength of nf-core/mag is not any single algorithm, but the integration and reproducibility of multiple tools within a standardized workflow. This is particularly important for large-scale metagenomic studies where software versions, parameters, databases and computational environments can strongly influence results.

After seven years of development involving multiple curator teams and the wider nf-core community, v5 demonstrates how community-driven workflow development can keep metagenomic analysis aligned with rapidly evolving sequencing technologies.

The future of genome-resolved metagenomics is not simply longer reads—it is better integration of assembly, binning, refinement and quality control within reproducible computational workflows.

Read the full paper in Bioinformatics https://academic.oup.com/bioinformatics/article/42/9/btag628/8770536?

Analysis of regulatory sequences (RSAT project)
Classification and analysis of mobile genetic elements (ACLAME project).
Analysis of molecular interaction networks (NeAT project)
Inference of metabolic pathways from genomic and post-genomic data
(metabolic pathfinding, see also metabolic pathfinding in NeAT)
Critical assesment of protein interactions (CAPRI)

Lab Page http://www.bigre.ulb.ac.be/