Build a clade collection
In this tutorial, we will create a clade project including all the genomes available for a species in RefSeq as well as any additional genomes you may have using MiGA alone. If you want to explore a more manual approach using
bash, see the Build a clade collection using BASH example. We will use Escherichia coli as the target species, but you can use any species (or any taxon) you want.
miga new -P E_coli -t clade
There are different stages of completeness defined in the NCBI Genome database, and you may want to include only some cases depending on your analysis. The stages (from higher to lower quality) are:
- 1.Complete: Genomes including all replicons in the organism(s) sequenced.
- 2.Chromosome: Genomes with complete chromosome (but missing other replicons).
- 3.Scaffold: Draft genomes with scaffold status.
- 4.Contig: Draft genomes with contig status.
In this example, we'll skip the draft genomes. However, if you want all of them, simply use the
Re-running and updating: If the following code fails at any point, for example due to a network interruption, you can simply re-run it, and it will take it from where it failed.
miga ncbi_get -T "Escherichia coli" -P . --complete --chromosome -v
Note that you can change the value of
"Escherichia coli"to any other species name, or even taxa of any rank such as genus or family.
miga ncbi_get -T "Escherichia coli" -P . --complete --chromosome --api-key ABCD123 -v
Or you can set it globally as an environmental variable before running
If you have any unreleased genomes, you can simply add them to the same project to be processed together with those publicly available. You can initialize datasets at different points, see input data. For the purposes of this tutorial, we'll assume that you have raw coupled reads from two sequencing lanes (1 and 2) in Gzipped FastQ files:
# Unzip and/or concatenate input files
# this is not necessary if your files are ready to use:
gzip -d -c ~/some/file/d1_ACTG_L_R1.fastq.gz > /tmp/sister.1.fastq
gzip -d -c ~/some/file/d1_ACTG_L_R2.fastq.gz > /tmp/sister.2.fastq
# Register the dataset
# change the dataset name MyDataset to whichever name you want:
miga add -P . -D MyDataset -t genome -i trimmed_reads_paired /tmp/sister..fastq
miga daemon start -P .