More at https://www.nature.com/articles/s41467-021-22671-6

Address of the bookmark: https://github.com/dmcblab/InfoGenomeR

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

Install

pip install genomenotebook

How to use

Create a simple genome browser with a search bar. The sequence appears when zooming in.

import genomenotebook as gn

g=gn.GenomeBrowser(genome_path, gff_path, init_pos=10000)
g.show()

Tracks can be added to visualize your favorite genomics data. See Examples for more !!!!

Address of the bookmark: https://dbikard.github.io/genomenotebook/

What is PhyloHerb: PhyloHerb is a wrapper program to process genome skimming data collected from plant materials. The outcomes include the plastid genome (plastome) assemblies, mitochondrial genome assemblies, nuclear ribosomal DNAs (NTS+ETS+18S+ITS1+5.8S+ITS2+28S), alignments of gene and intergenic regions, and a species tree. It is designed to be a high throughput program dealing with lower quality data. Examples include low-coverage (5x cpDNA) plastome phylogeny, recycling plastid genes from target enrichment data, retrieving low-copy nuclear genes from medium coverage (5x nucDNA) genome skimming.

License: GNU General Public License

Citation:

• Cai, Liming, Hongrui Zhang, and Charles C. Davis. 2022. PhyloHerb: A high‐throughput phylogenomic pipeline for processing genome‐skimming data. Applications in Plant Sciences 10(3): 1–9. https://doi.org/10.1002/aps3.11475

Address of the bookmark: https://github.com/lmcai/PhyloHerb/

1. Earth BioGenome Project (EBP) - A global effort to sequence and annotate roughly 1.5 million known eukaryotic species and create a comprehensive reference library of Earth's biodiversity.
2. Vertebrate Genomes Project (VGP) - Focuses on producing high-quality reference genomes for vertebrate species, providing a major foundation for the wider Tree of Life.
3. Darwin Tree of Life - A UK and Ireland initiative aiming to sequence approximately 70,000 species, creating a detailed genomic map of regional biodiversity.
4. African BioGenome Project (AfricaBP) - A pan-African initiative focused on sequencing Africa's biodiversity while strengthening genomics and bioinformatics capacity across the continent.
5. European Reference Genome Atlas (ERGA) - A European effort to generate high-quality reference genomes for Europe's biodiversity and connect national sequencing efforts.
6. Canada BioGenome Project - Building reference genomes for Canadian biodiversity, supporting conservation, research and the study of ecosystems.
7. California Conservation Genomics Project - Uses genomics to understand and protect California's biodiversity and help inform conservation decisions.
8. Bat1K - An ambitious global effort to sequence the genomes of all known bat species, helping scientists study their evolution, longevity, immunity and unique biology.
9. Bird 10,000 (B10K) - A global project aiming to create genome resources covering the world's bird diversity and reconstruct avian evolutionary history.
10. Global Invertebrate Genomics Alliance (GIGA) - Builds genomic resources for the enormous and genetically diverse world of invertebrates.
11. GenomeIndia - India's national human-genome initiative, designed to capture the country's extraordinary population diversity and create a reference resource for Indian genomics.
12. All of Us Research Program - A US national research program combining genomic information with health and lifestyle data at population scale.
13. 1+ Million Genomes / Genome of Europe - A European effort to enable secure, cross-border use of genomic and health data and develop genome-scale resources involving more than one million people.
14. Korea National Integrated Bio Big Data Project - South Korea's large-scale national effort combining genomic and health information from a very large population cohort.
15. Saudi Genome Program - A national genomics initiative focused on genetic diseases, population genomics and precision medicine in Saudi Arabia.
16. Australian Genomics / Genomics Health Futures Mission - Australia's growing investment in genomic medicine, rare disease research and population-scale health genomics.

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/