Lee Katz

I am the senior bioinformatician in the Enteric Diseases Laboratory Branch at CDC. In my free time, I try to be a good husband and parent, and in my imaginary free time, I podcast with my co-hosts on the Micro Binfie podcast.

less than 1 minute read

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.com/lskatz/lskScripts/blob/master/qsub/launch_kraken.sh

The main points of the script are shown below.

run $KRAKENDIR/kraken --fastq-input --paired --db=$KRAKEN\_DEFAULT\_DB --preload --gzip-compressed --quick --threads $NSLOTS --output $KRAKENOUT $READS  
  
run kraken-translate --db $KRAKEN\_DEFAULT\_DB $KRAKENOUT | cut -f 2- | sort | uniq -c |\\  
  perl -lane '  
              s/^ +//;   # remove leading spaces  
              s/ +/\\t/;  # change first set of spaces from uniq -c to a tab  
              s/;/\\t/g;  # change the semicolon-delimited taxonomy to tab-delimited  
              print;  
             ' |\\  
  sort -k1,1nr > $KRAKENTAXONOMY  
  
  
run $KRONADIR/ktImportText -o $HTML $KRAKENTAXONOMY

Updated:

You may also enjoy

Pixi basics Permalink

less than 1 minute read

I wrote a post on the microbifie blog on how to get started with Pixi.

Call for action: bioinformatics conferences

less than 1 minute read

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

10 minute read

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...