Easy: predicting road flooding using open sensor data & measured datums (levelslevels.uk) Almost impossible: remotely ascertaining whether a given Somerset Council flood gate is open or closed, given their seldom updated website (currently 9 days since update) www.somerset.gov.uk/roads-travel...

Some bioinformatics papers: this lightweight method requires 4 gigabytes of RAM. Other bioinformatics papers: this lightweight method requires 4 terabytes of RAM.

We are recruiting for a number of PhDs in the HPRU in Public Health Genomics, based at @unibirmingham.bsky.social, in collaboration with UKHSA. 🦠🧬 Apply below to work on cutting-edge science and public health priorities! (NB deadline of 9th January / UK only due to funding restrictions) (1/7)

This is Bam, shall I bring him in? Bites more than Sam.

Not the first paper I've been a part of, but the first one I've ever written, organized, and executed! I learned a lot getting this one done. Also, Mojo! 🔥 👇

Amazing novel mpox recombination between epidemic clade 1 and 2 in returning traveller detected in UK: virological.org/t/inter-clad...

'Can you stop an outbreak becoming a pandemic?' ARTIC Network recently delivered a workshop at the PHA4GE conference in Cape Town. Read more: www.linkedin.com/posts/artic-... @pha4ge.bsky.social

"Prolonged exposure in shared, poorly ventilated spaces… drives respiratory virus transmission more than close contact." www.nature.com/articles/s41...

We're half a decade into studies finding that improving airflow in classrooms will reduce disease transmission enormously, and that bleaching surfaces etc. does very little. And yet nothing changes. Waves of flu and colds wash over schools, and the schools pretend it's an act of God.

And if you have a .fastq.gz.mim already present, you get another up to 2x speedup (or more with >6 cores).

TIL Deacon decompresses fastq.gz 3x faster than GNU coreutils gzip/zcat. You can decompress, dehost and recompress fastq.gz using Deacon in less time than only decompressing using coreutils gzip/zcat or pigz.

github.com/bede/deacon For anyone still using Bowtie2 for filtering or depletion of host sequences or specifics, I can recommend Deacon from @bedec.bsky.social . It is so much faster and easier than Bowtie2, and its performance is equal or better (tested with metagenomes and mitogenomes).🧬 & 🖥️

Minimizer Density revisited: Models and Multiminimizers https://www.biorxiv.org/content/10.1101/2025.11.21.689688v1

579 high-quality human genomes from @humanpangenome.bsky.social, Arab Pangenome and individual papers (CHM13, CN1, KSA001, I002C, YAO and KOREF1). Sequences available in the AGC format (3.7GB) and FM-index in the ropebwt3 format (20.3GB). For details, see github.com/lh3/human-asm

Ok; mim (github.com/COMBINE-lab/...) preprint submitted! Excited for folks to see it and share thoughts. The key takeaway; mim allows the quick, one-time, building of a small auxiliary index that then allows scaling gzipped FASTQ parsing linearly in # of threads. 1/2

@wytamma.bsky.social : so, it took a little bit of extra time (not the flight back from the CZI meeting), but I decided to just f#&$ing do it, and the basic code to build and parse with the auxiliary fastq index is working (github.com/COMBINE-lab/...). 1/2

New preprint: we looked into production of the bacterial toxin colibactin and found that MDR E. coli from the global north have co-evolved with endemic colibactin producers, acquiring colibactin resistance genes before undergoing clonal expansions. www.biorxiv.org/content/10.1...

I want to spell this out in case the implications aren't clear: This means all public tools/webapps of GISAID data (all the ones you've been used to seeing thru the pandemic, as far as we can tell) are prohibited. The file allowed this. Cut that - cut off all tools the public & others were using.

My account's upload and bulk download access were terminated permanently in 2021 without explanation after I published *checksums* of GISAID genomes. GISAID and its SAB have since ignored a dozen emails seeking explanation. 4 yrs on, even Nextstrain has lost access. GISAID has rotted from its core.

I was on Last Word, the Radio 4 obituary programme, trying to sum up Jim Watson’s near-century long life.

Long term, good software requires lazy users: complain and file issues as soon as things don't work first try. If nothing else, it means documentation should be improved.

🔸️Early data suggests we could be in for a worse than normal flu season, brought on by a cluster of escape mutations in H3N2 this year that may lift the Re from 1.2 to 1.4. We are starting to see an uptick in the US. Data from BIOFIRE.

New post from me, for UK folks only, on how you need to start preparing for Apple to switch off Advanced Data Protection and the end-to-end encryption of the data you store on it. Like I said, UK only. #SunlitUplands heatherburns.tech/2025/11/10/t...

As expected, unfortunately. If ever you needed a reason for never using GISAID ever again (as a data producer or data user - we're both), look no further. Time to move on to more trusted and transparent solutions.

Our method for genome size estimation from long-read overlaps is now published 🥳 academic.oup.com/bioinformati...

“my brain is open” users.monash.edu/~normd/docum...

Really exciting that the preprint on Barbell, a new demultiplexer, is finally out! It's the first tool that builds on Sassy, the approximate-DNA-searching tool that @rickbitloo.bsky.social and myself developed earlier this year, specifically with this application in mind.

RIFs at CDC. Destroying the Epidemic Intelligence Service, NCIRD, NCIPC & a dozen other areas/divisions/branches will cause enormous harm and suffering to America, and indeed to the whole world. Call your Reps about it. Then call them again.

Our recent paper on rifampicin resistant subpopulations in M. tuberculosis (M. tb) has been published at JAC-antimicrobial resistance. I am really happy to see this work published just hours before submitting my DPhil thesis! 🔗👇 doi.org/10.1093/jaca...

For more information about the Friday night massacre at CDC, I wrote up an analysis of who got terminated and what that means for public health. Grateful to @saveamericamvmt.bsky.social for supporting and amplifying. We are in really terrible trouble. rasmussenretorts.substack.com/p/the-death-...

So what's the equivalent of `perf record && perf report` on a MacBook? I want to see the generated assembly and which lines are hot.

Last week we were in the Washington Post for our characterization of Robertsonian chromosomes. This week we are entering our 10th day of being shut down and all of our research is on hold. To help me feel not-so-bad, here is a thread of some studies we released right before the shutdown 🧵 [1/n]...

Funny story, though, we found this gene in NCBI databases, but it was annotated in Streptococcus pneumoniae! This is surely human contamination in a bacterial strep sample that was not properly filtered. Lesson: use CHM13, or better yet a pangenome, when filtering for human contamination...

Just published an interactive article about a magical algorithm known as the Burrows-Wheeler Transform, which powers sequence alignment tools like bowtie and bwa: sandbox.bio/concepts/bwt It's also notoriously unintuitive so I'm hoping this article helps you build that intuition.

I've added 7 videos to my Burrows-Wheeler indexing playlist (www.youtube.com/playlist?lis...), rounding out the r-index series and adding a 5-part series on the move structure. Now 27 videos in that playlist. I aim to add videos on prefix-free parsing, PBWT, Wheeler languages/automata in the future.

Deacon 0.11.0: - Local server mode - Ultra-careful handling of non-ACGT - Faster indexing & index loading - Denser index now stores k-mers not hashes - xxHash & FxHash replaced with rapidhash::fast - Bug fixes Thanks @curiouscoding.nl (and others!) for contributions github.com/bede/deacon/...

"OpenZL is our answer to the tension between the performance of format-specific compressors and the maintenance simplicity of a single executable binary." engineering.fb.com/2025/10/06/d...

Looking for people to test the latest version of simd-sketch. It's now 2x as fast at sketching, and supports skipping over kmers containing N and other ambiguous bases (which is only ~35% slower). 'cargo install simd-sketch' is right there under your fingertips ;) github.com/RagnarGrootK...

FxHashSet::<u32>::contains throughput is wild! - Up to 4x slowdown for negative queries due to probing. - Positive queries are fast for small tables, but slow in RAM because they need 2 cache misses. Lots of variance depending on the load factor, ie whether n is close to 87.5% of a power of 2.

Pleased to see this pre-printed, highlighting the completeness/accuracy of @nanoporetech.com long-read genome assembly for clinical Enterobacterales: www.biorxiv.org/content/10.1... Thanks to colleagues @modmedmicro.bsky.social, @ukhsa.bsky.social, @genewiz.bsky.social and @oxfordbrc.bsky.social!

Terrific new feature presented by @theo.io on @pathoplexus.org called SeqSets for generating DOIs for sequence subsets used in publications, that can then be tracked for impact via CrossRef that will allow data generators to track impact! #IMMEMXiV

Front page of Hacker News 🫨 bsky.app/profile/bede...

Blogged about how zstd --long fills the gap between fast and slow-but-high-ratio genome compression methods log.bede.im/2025/09/12/z...

First rate BBC journalism that somehow didn't make the cut for #r4today www.bbc.co.uk/news/article... @andyverity.bsky.social

Sometimes you meet absolutely incredible bioinfo-magicians. It was a huge privilege when @shenwei356.bsky.social joined our group for a year on an @embl.org sabbatical. While here, he developed a new way of aligning to millions of bacteria, called LexicMap 1/n www.nature.com/articles/s41...

Zstandard's --long range mode works wonders for assemblies, but needs uninterrupted single line sequences. *AllTheBacteria 661k, multiline fasta* gzip (pigz): 751GB zstandard --long: 641GB (30% original size) *Single line fasta* gzip (pigz): 700GB zstandard --long: 232GB (10% original size)

"The Agriculture Department will not require poultry companies to limit salmonella bacteria in their products, halting a Biden Administration effort to prevent food poisoning from contaminated meat." www.seattletimes.com/business/usd...

For anyone who has used pling for comparing plasmids using rearrangement distances ("how many structural events apart are these plasmids"), here's how to tweak parameters, and integrate it with typing info, and the host phylogeny www.biorxiv.org/content/10.1... github.com/iqbal-lab-or...

Preprint out for myloasm, our new nanopore / HiFi metagenome assembler! Nanopore's getting accurate, but 1. Can this lead to better metagenome assemblies? 2. How, algorithmically, to leverage them? with co-author Max Marin @mgmarin.bsky.social, supervised by Heng Li @lh3lh3.bsky.social 1 / N

Claude's kryptonite: incredibly mundane biology questions