Genetics · Module 14 · Genomics

Genomics Lab

Shatter a genome and assemble it, watch repeats break the assembly, run a GWAS and see which SNPs cross the genome-wide line, practise coverage, ORF, odds ratio, fold change and gene tree problems with animated worked solutions, match technologies to their data, and read real genomics results. Nothing here is graded.

Part 1: Shotgun Assembly Bench

Unique sequenceRepeat copy (R1, R2) ReadContigUncovered base / gapCoverage depth
Assembly results and Lander-Waterman predictions
MeasureThis runTheory

Simplifications: a 10 kb linear genome; reads land uniformly at random with no sequencing errors (reads may hang off either end, so every base has the same expected coverage C); two reads join when they overlap by at least 20 bp. A repeat copy is resolved only when some read covers it with at least 20 bp of unique sequence on both sides; otherwise the assembly breaks there, and a greedy assembler may misjoin the flanks. Theory uses the Lander-Waterman model: uncovered fraction e−C and expected contigs 1 + N·e−C(1−θ/L) (one, plus the expected number of gaps), with θ = 20 bp.