Skip to content

Instantly share code, notes, and snippets.

@awitney
Created June 25, 2019 11:54
Show Gist options
  • Select an option

  • Save awitney/7ac35d3fcdecbff5ffc254bd6e526475 to your computer and use it in GitHub Desktop.

Select an option

Save awitney/7ac35d3fcdecbff5ffc254bd6e526475 to your computer and use it in GitHub Desktop.
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