This is a limited proof of concept to search for research data, not a production system.

Search the MIT Libraries

Title: roblanf/sarscov2phylo: 3-11-20

Type Software roblanf, Richard Mansfield (2020): roblanf/sarscov2phylo: 3-11-20. Zenodo. Software. https://zenodo.org/record/4279921

Authors: roblanf ; Richard Mansfield ;

Links

Summary

Citation and reuse

Please cite this release as:

Lanfear, Rob (2020). A global phylogeny of SARS-CoV-2 sequences from GISAID. Zenodo DOI: 10.5281/zenodo.3958883

You can visit that DOI here:

If you publish papers that use this tree you must still follow the GISAID data sharing and attribution rules.

Details

The trees in this release were generated with the following command line:

bash global_tree_gisaid_start_tree.sh -i [gisaid.fasta] -p [previous_iteration] -t 250

[gisaid.fasta] is the fasta file of high coverage and complete raw sequences from GISAID up to and including the date in the title of the release, determined by the 'submission date' filter on a GISAID data feed

[previous_iteration] is the filepath of the previous release, this is used to provide the excluded_sequences.tsv and ft_SH.tree files as the starting points of the current iteration.

Filtering statistics sequences downloaded from GISAID 131215 // alignment stats of global alignment Alignment number: 1 Format: aligned FASTA Number of sequences: 128821 Alignment length: 29903 Total # residues: 3839063069 Smallest: 29105 Largest: 29903 Average length: 29801.5 Average identity: 100% // alignment stats of global alignment after masking sites Alignment number: 1 Format: aligned FASTA Number of sequences: 128821 Alignment length: 29903 Total # residues: 3821789644 Smallest: 29036 Largest: 29675 Average length: 29667.4 Average identity: 100% // alignment stats after filtering out short/ambiguous sequences Alignment number: 1 Format: aligned FASTA Number of sequences: 128777 Alignment length: 29903 Total # residues: 3820487101 Smallest: 29036 Largest: 29675 Average length: 29667.5 Average identity: 100% // alignment stats of global alignment after trimming sites that are >50% gaps Alignment number: 1 Format: aligned FASTA Number of sequences: 128777 Alignment length: 29646 Total # residues: 3810742331 Smallest: 28498 Largest: 29646 Average length: 29591.8 Average identity: 100% // After filtering sequences with TreeShrink Type: Phylogram #nodes: 224247 #leaves: 128714 #dichotomies: 90656 #leaf labels: 128714 #inner labels: 85123 Number of new sequences added this iteration 6018 alignment_names_new.txt Notable changes to the scripts in this release None Notable aspects of the trees There is some notable lack of clustering in some pangolin lineages, but I have not yet had a chance to figure out whether this is an issue with my recently added QC code (quite possible), the tree, or something else.

More information

  • DOI: 10.5281/zenodo.4279921

Dates

  • Publication date: 2020
  • Issued: November 18, 2020

Rights

  • info:eu-repo/semantics/openAccess Open Access

Much of the data past this point we don't have good examples of yet. Please share in #rdi slack if you have good examples for anything that appears below. Thanks!

Format

electronic resource

Relateditems

DescriptionItem typeRelationshipUri
IsSupplementTohttps://github.com/roblanf/sarscov2phylo/tree/3-11-20
IsVersionOfhttps://doi.org/10.5281/zenodo.3958883
IsPartOfhttps://zenodo.org/communities/covid-19
IsPartOfhttps://zenodo.org/communities/zenodo