Pixi basics Permalink
I wrote a post on the microbifie blog on how to get started with Pixi.
Bioinformatics blog
2024
Call for action: bioinformatics conferences
Every fiscal year, we have to predict where we are going to travel over the whole next year. However, conferences mainly post on their own website, making no...
Hash collision experiment
I have talked a bit about hashes with MLST but one crucial issue with hashes is whether they have collisions. In other words, if we have two sequences, there...
What does it take to switch over to BAM from FASTQ? Permalink
I wrote a post on the microbifie blog.
Hash MLST analysis
With EToKi and perhaps future MLST callers, we are possibly looking at a new class of MLST databases 1. I’m not sure which caller was actuall...
2023
List of bioinformatics conferences
Every fiscal year, we have to predict where we are going to travel over the whole next year. However, conferences mainly post on their own website, making no...
Converting a set of VCFs to distances
I think that converting a bunch of VCF files to an alignment is still an artisinal thing. I made this sort of task in the Lyve-SET pipeline back in 2013 and ...
Open source alternatives
BioNumerics is sunsetting in December 2024. This is an issue especially in the PulseNet community where the network relies on this software to do many things...
MLST with colorid
I am pleasantly surprised by how well ColorID works with MLST. For some background, please see my previous post on cgMLST.
We can make this beautiful and searchable
You might have seen this message on GitHub for a tab separated file.
Benchmark datasets
Over the last several years, we have been generating benchmark datasets.
The memories
I was cleaning out my office and found a flag that was given to me a few years ago. The legendary Peter Gerner-Smidt had become a citizen and had given some ...
How to call whole genome MLST on the command line
I have been asked in many different contexts how to call alleles on the command line. In August 2022, I gave an internal workshop on how to call alleles on t...
How to be a good citizen on Linux
This is a short tutorial on how to be a good citizen in a shared Linux environment. At my own institution, there are many users per computer and we all have ...
2022
mstdn.science
I am using my site to verify me on Mastodon. Nabil and Duncan spun up a Mastodon server at https://mstdn.science a few weeks ago and it has been wildly succe...
GenBank CV
I was inspired by Andrey Kislyuk’s post about 15 years ago on his personal website where he formatted his personal website to look like a GenBank file. I rec...
2021
NFT of BRCA1
I’ve been sort of amazed ever since I heard about NFTs from Planet Money and now the New York Times. Basically you associate a place on the blockchain with ...
2020
Downloading the breadth of SARS-CoV-2
I was trying to figure out how to download the breadth of all of SARS-CoV-2 genomes and so I started out with the two major repositories: NCBI and GISAID.
2019
Minimizer perl module
I looked all over cpan but did not find a module to make Minimizers. Therefore I went ahead and made it in Bio::Minimizer.
Cluster detection using Mash
I have been racking my brain on what the command line version of a rapid cluster detection should be. As in, how can I rapidly figure out if two random geno...
Fixing the R package manager
R’s package manager is broken. I present to you my solution for it. $HOME is your home directory. 1) mkdir -pv $HOME/R/tmp 2) Edit $HOME/.Renviron and ad...
Mash perl module
Hey y’all, I made a new perl module to read and write Mash files. I made it a complete package, adding it to CPAN with documentation and adding unit testing...
2017
Sampling the taxonomy database
I was a little frustrated that every time I wanted to try out my new Bio::DB::Taxonomy-based script, it would take a few minutes to run…. and then I would fi...
Custom Kraken database
Like many labs around the world, we use Kraken for contamination detection. This isn’t its intended purpose however because it is supposed to be a taxonomic...
2016
Tree To Reads
It looks like Tree To Reads is online in a preliminary draft! I encourage everyone to try it out! http://biorxiv.org/content/early/2016/01/22/037655 Using...
Perl writing style
I’m taking a moment to reflect on a post that Torsten Seemann wrote a couple of years ago called “Minimum standards for bioinformatics command line tools.” ...
2015
Krona snippet
I couldn’t find any kind of conversion Kraken and Krona recently and so I wrote up a little pipeline. Full script is available at lskScripts: https://github...
Edirect snippets
I have been really excited by the new-ish edirect utilities. I thought I’d put out my snippets and encourage anyone to show theirs too. I haven’t found a l...
The confusing vocabulary of public health and bioinformatics
In the very beginning of our Listeria project, we basically had an idea that we could use evolution to help in epidemiological investigations. The idea was ...
Blog posts
A few blog posts that I’ve really liked for my niche of public health and bioinformatics. I think one purpose of this blog can be to start gathering these p...
FASTQ must die! Long live SAM/BAM! (seconded!)
Here is a short rant but it becomes proactive. For those of you who don’t know me, I dislike having two fastq files per genome if they are paired end (PE). ...