2020年10月28日 软件 开放存取

罗布兰夫; 理查德·曼斯菲尔德

MARC21 XML导出

<?xml version='1.0' encoding='UTF-8'?>
<record xmlns="http://www.loc.gov/MARC21/slim">
  <leader>00000nmm##2200000uu#4500</leader>
  <controlfield tag="005">20201029002700.0</controlfield>
  <controlfield tag="001">4145970</controlfield>
  <datafield tag="700" ind1=" " ind2=" ">
    <subfield code="a">Richard Mansfield</subfield>
  </datafield>
  <datafield tag="856" ind1="4" ind2=" ">
    <subfield code="s">9831507</subfield>
    <subfield code="z">md5:54af4da08b0fdc8ba5d7c635e1c5cafc</subfield>
    <subfield code="u">//americinnmankato.com/record/4145970/files/roblanf/sarscov2phylo-18-10-20.zip</subfield>
  </datafield>
  <datafield tag="542" ind1=" " ind2=" ">
    <subfield code="l">open</subfield>
  </datafield>
  <datafield tag="260" ind1=" " ind2=" ">
    <subfield code="c">2020-10-28</subfield>
  </datafield>
  <datafield tag="909" ind1="C" ind2="O">
    <subfield code="p">software</subfield>
    <subfield code="p">user-covid-19</subfield>
    <subfield code="p">user-zenodo</subfield>
    <subfield code="o">oai:zenodo.org:4145970</subfield>
  </datafield>
  <datafield tag="100" ind1=" " ind2=" ">
    <subfield code="a">roblanf</subfield>
  </datafield>
  <datafield tag="245" ind1=" " ind2=" ">
    <subfield code="a">roblanf / sarscov2phylo：18-10-20</subfield>
  </datafield>
  <datafield tag="980" ind1=" " ind2=" ">
    <subfield code="a">user-covid-19</subfield>
  </datafield>
  <datafield tag="980" ind1=" " ind2=" ">
    <subfield code="a">user-zenodo</subfield>
  </datafield>
  <datafield tag="540" ind1=" " ind2=" ">
    <subfield code="a">Other (Open)</subfield>
  </datafield>
  <datafield tag="650" ind1="1" ind2="7">
    <subfield code="a">cc-by</subfield>
    <subfield code="2">opendefinition.org</subfield>
  </datafield>
  <datafield tag="520" ind1=" " ind2=" ">
    <subfield code="a">Citation and reuse
&lt;p&gt;Please cite this release as:&lt;/p&gt;
&lt;blockquote&gt;&lt;p&gt;Lanfear, Rob (2020). A global phylogeny of SARS-CoV-2 sequences from GISAID.  Zenodo 土井： 10.5281 / zenodo.3958883&lt;/p&gt;
&lt;/blockquote&gt;
&lt;p&gt;You can visit that 土井 here: &lt;a href="//americinnmankato.com/badge/latestdoi/260124648"&gt;&lt;/a&gt;&lt;/p&gt;
&lt;p&gt;If you publish papers that use this tree you must still follow the GISAID data sharing and attribution rules.&lt;/p&gt;
Details
&lt;p&gt;The trees in this release were generated with the following command line:&lt;/p&gt;
&lt;p&gt;&lt;code&gt;bash global_tree_gisaid_start_tree.sh -i [gisaid.fasta] -p [previous_iteration] -t 250&lt;/code&gt;&lt;/p&gt;
&lt;ul&gt;
&lt;li&gt;&lt;p&gt;[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&lt;/p&gt;
&lt;/li&gt;
&lt;li&gt;&lt;p&gt;[previous_iteration] is the filepath of the previous release, this is used to provide the &lt;code&gt;excluded_sequences.tsv&lt;/code&gt; and &lt;code&gt;ft_SH.tree&lt;/code&gt; files as the starting points of the current iteration.&lt;/p&gt;
&lt;/li&gt;
&lt;/ul&gt;
Filtering statistics
&lt;pre&gt;&lt;code&gt;sequences downloaded from GISAID
108202
//
alignment stats of global alignment
Alignment number:    1
Format:              aligned FASTA
Number of sequences: 106464
Alignment length:    29903
Total # residues:    3173059520
Smallest:            29030
Largest:             29903
Average length:      29804.1
Average identity:    100%
//
alignment stats of global alignment after masking sites
Alignment number:    1
Format:              aligned FASTA
Number of sequences: 106464
Alignment length:    29903
Total # residues:    3158473719
Smallest:            28961
Largest:             29675
Average length:      29667.1
Average identity:    100%
//
alignment stats after filtering out short/ambiguous sequences
Alignment number:    1
Format:              aligned FASTA
Number of sequences: 106273
Alignment length:    29903
Total # residues:    3152816993
Smallest:            28961
Largest:             29675
Average length:      29667.1
Average identity:    100%
//
alignment stats of global alignment after trimming sites that are &amp;gt;50% gaps
Alignment number:    1
Format:              aligned FASTA
Number of sequences: 106273
Alignment length:    29646
Total # residues:    3144363035
Smallest:            28437
Largest:             29646
Average length:      29587.6
Average identity:    100%
//
After filtering sequences with TreeShrink
Type:   Phylogram
#nodes: 186198
#leaves:    106272
#dichotomies:   75928
#leaf labels:   106272
#inner labels:  71407
Number of new sequences added this iteration
340 alignment_names_new.txt
&lt;/code&gt;&lt;/pre&gt;
Notable changes to the scripts in this release
&lt;ul&gt;
&lt;li&gt;Some more issues occurred with sequence names, so I have shifted to using exclusively just the EPI-IDs for sequence names in the alignments, trees, and scripts. This should avoid any future issues with sequence names. &lt;/li&gt;
&lt;/ul&gt;
&lt;p&gt;If you have a GISAID data stream, this means you can get the data as follows:&lt;/p&gt;
&lt;pre&gt;&lt;code&gt;curl -u "${GISAID_LOGIN}:${GISAID_PASSWORD}" ${URL} 2&amp;gt;/dev/null \
    | xz -d -T22 \
    | jq -r 'select(.covv_subm_date &amp;lt;= "'"$DATE"'" and .is_complete == true and .is_high_coverage == true) | ["&amp;gt;" + .covv_accession_id + "\n" + .sequence] | join("\n")' \
    | xz -c -T22 &amp;gt; gisaid-${DATE}.fasta.xz
&lt;/code&gt;&lt;/pre&gt;
Notable aspects of the trees
&lt;ul&gt;
&lt;li&gt;Names have changed to include only the EPI-IDs. This has the side-benefit of roughly halving the size of the tree file.&lt;/li&gt;
&lt;/ul&gt;</subfield>
  </datafield>
  <datafield tag="773" ind1=" " ind2=" ">
    <subfield code="n">url</subfield>
    <subfield code="i">isSupplementTo</subfield>
    <subfield code="a">//github.com/roblanf/sarscov2phylo/tree/18-10-20</subfield>
  </datafield>
  <datafield tag="773" ind1=" " ind2=" ">
    <subfield code="n">doi</subfield>
    <subfield code="i">isVersionOf</subfield>
    <subfield code="a">10.5281 / zenodo.3958883</subfield>
  </datafield>
  <datafield tag="024" ind1=" " ind2=" ">
    <subfield code="a">10.5281 / zenodo.4145970</subfield>
    <subfield code="2">doi</subfield>
  </datafield>
  <datafield tag="980" ind1=" " ind2=" ">
    <subfield code="a">software</subfield>
  </datafield>
</record>