Skip to content

Instantly share code, notes, and snippets.

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
Lineage report
--------------
Lineage Estimated Fraction Family Spoligotype Rd
lineage4 1.000 Euro-American LAM;T;S;X;H None
lineage4.3 1.000 Euro-American (LAM) mainly-LAM None
lineage4.3.4 1.000 Euro-American (LAM) LAM RD174
lineage4.3.4.2 1.000 Euro-American (LAM) LAM1;LAM4;LAM11 RD174
@awitney
awitney / gist:60b49d0a2cb8ab74d3376f955ee03721
Created January 11, 2019 21:26
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
@awitney
awitney / gist:6827524
Created October 4, 2013 15:09
fetching genes from ensembl
#!/usr/bin/env perl
use strict;
use warnings;
use Bio::EnsEMBL::LookUp;
use Bio::SeqIO;
use Bio::Seq;
my $string = $ARGV[0] || 'gyrA';