Skip to content

Instantly share code, notes, and snippets.

@awitney
Created January 11, 2019 21:26
Show Gist options
  • Select an option

  • Save awitney/60b49d0a2cb8ab74d3376f955ee03721 to your computer and use it in GitHub Desktop.

Select an option

Save awitney/60b49d0a2cb8ab74d3376f955ee03721 to your computer and use it in GitHub Desktop.
mycobacterium build kraken db
rsync -avL rsync://ftp.ncbi.nlm.nih.gov/genomes/refseq/bacteria/Mycobacterium_* refseq/Mycobacterium/
rsync -avL rsync://ftp.ncbi.nlm.nih.gov/genomes/refseq/bacteria/Mycolicibacterium_* refseq/Mycolicibacterium/
mkdir fna
find refseq/Mycobacterium/ -name "*_genomic.fna.gz" -exec cp {} fna \;
find refseq/Mycolicibacterium/ -name "*_genomic.fna.gz" -exec cp {} fna \;
rm -f fna/*_cds_from_genomic.fna.gz
rm -f fna/*_cds_from_genomic.fna.gz
rm -f fna/*_rna_from_genomic.fna.gz
gunzip fna/*.gz
for i in fna/*.fna ; do echo $i; kraken2-build --add-to-library $i --db mycobacterium; done
kraken2-build --download-taxonomy --db mycobacterium
kraken2-build --build --db mycobacterium
rm -r fna
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment