roblanf/sarscov2phylo: 16-08-2020
roblanf · Zenodo (CERN European Organization for Nuclear Research) · 2020
Citation and reuse Please cite this release as: Lanfear, Rob (2020). A global phylogeny of SARS-CoV-2 from GISAID data, including sequences deposited up to 16-August-2020. Zenodo DOI: 10.5281/zenodo.3958883 You can visit that DOI here: Please note - you cannot publish papers that use this tree without following the GISAID data sharing and attribution rules. These rules are important - they protect the data uploaders, and create trust in a global system of data sharing with potentially vast public health benefits. By building and maintaining trust we ensure that people keep sharing their data, and that the public health benefits keep flowing. I do not want the existence of this tree to be some kind of attribution laundering service (e.g. where people feel free to use the tree without following the GISAID data sharing rules), so please don't use it in that way. For example, if you are going to interpret other people's data from GISAID and publish the results, including by using this tree, you should get in touch with the people that submitted the data. The code in this repo is covered by the GNU license, and you can use that however you like. Details The trees in this release were generated with the following command line: bash global_tree_gisaid_start_tree.sh -i gisaid_hcov-19_2020_08_18_22.fasta -o global.fa -s ft_SH_14-08-20.tree -t 54 The raw sequence file contains all available SARS-CoV-2 genomes in GISAID available on the 16th of August 2020, determined by the 'submission date' filter on GISAID. The ZIP file contains the code necessary to reproduce the trees themselves, and the README in the zip file also describes the methods used in detail. I also include the tree itself here so that it can be easily downloaded without downloading the entire repo. The file ' ft_SH_14-8-20.tree' is the 'ft_SH.tree' file from the 14-8-20 release. The lnL of the final tree is: -526648.019 Filtering statistics sequences downloaded from GISAID 53910 // alignment stats of global alignment Alignment number: 1 Format: aligned FASTA Number of sequences: 52730 Alignment length: 29903 Total # residues: 1573059735 Smallest: 29030 Largest: 29903 Average length: 29832.3 Average identity: 100% // alignment stats of global alignment after masking sites Alignment number: 1 Format: aligned FASTA Number of sequences: 52730 Alignment length: 29903 Total # residues: 1571796972 Smallest: 29030 Largest: 29903 Average length: 29808.4 Average identity: 100% // alignment stats after filtering out short/ambiguous sequences Alignment number: 1 Format: aligned FASTA Number of sequences: 52463 Alignment length: 29903 Total # residues: 1563852377 Smallest: 29131 Largest: 29903 Average length: 29808.7 Average identity: 100% // alignment stats of global alignment after trimming sites that are >50% gaps Alignment number: 1 Format: aligned FASTA Number of sequences: 52463 Alignment length: 29782 Total # residues: 1558617245 Smallest: 28379 Largest: 29782 Average length: 29708.9 Average identity: 100% // After filtering sequences with TreeShrink Type: Phylogram #nodes: 94825 #leaves: 52436 #dichotomies: 40650 #leaf labels: 52436 #inner labels: 38514 Number of new sequences added this iteration 102 alignment_names_new.txt Notable changes to the scripts in this release None Notable aspects of the trees None