COMPARE_GENOMES.PY(1)

✨ NAME ✨

compare_genomes.py — find methylation in nanopore reads 💖

🧬 DESCRIPTION 🧬

okay bestie so here's the vibe: this lil script compares methylation between two genomes by looking at the probabilities of methylation for aligned positions. it's giving very much "compare and contrast but make it epigenetic" 🔥

here's how you call it, exactly:

usage: compare_genomes.py [-h] --bed1 BED1 --bed2 BED2 [-g GENOME_ALIGNMENT]
                          [-v]

🌸 optional arguments 🌸

• -h, --help — shows this help message and then dips, no notes 💅
• --bed1 BED1 — bed file 1, and it has to be the one with verbose output from make_bed.py — don't come in here with the quiet version, i'm begging 🙈
• --bed2 BED2 — bed file 2, same situation: verbose output from make_bed.py, no cap 🎀
• -g GENOME_ALIGNMENT, --genome_alignment GENOME_ALIGNMENT — an xmfa file from mauve 🦋 and listen, if you don't pass this one, then the alignments are assumed to be to the same reference genome. that's the whole lore, iykyk
• -v, --version — prints the version, that's it, that's the job ⭐

fr fr the two --bed files are the main character energy here — everything else is supporting cast 💜

🫶 sign-off 🫶

this iconic little page is from compare_genomes.py version 1.0.3+git202..., dated June 2021 — and yes the version string is literally cut off like that, we don't question it, we simply respect it 👑 periodt 💖✨