70e001e0f7a71f53c7a70b45feabac07d9d3751b
angie
  Mon Dec 14 16:35:16 2020 -0800
Update sarscov2phylo process to use tree and multi-seq ali's downloaded from GISAID instead of tree from github and aligning nextfasta locally.

diff --git src/hg/utils/otto/sarscov2phylo/updateSarsCov2Phylo.sh src/hg/utils/otto/sarscov2phylo/updateSarsCov2Phylo.sh
index 76ce0cc..0c66d53 100755
--- src/hg/utils/otto/sarscov2phylo/updateSarsCov2Phylo.sh
+++ src/hg/utils/otto/sarscov2phylo/updateSarsCov2Phylo.sh
@@ -1,40 +1,38 @@
 #!/bin/bash
 set -beEu -x -o pipefail
 
 #	Do not modify this script, modify the source tree copy:
 #	kent/src/hg/utils/otto/sarscov2phylo/updateSarsCov2Phylo.sh
 
 usage() {
-    echo "usage: $0 releaseLabel metadata_date.tsv.gz sequences_date.fasta.gz epiToPublicAndDate.date"
+    echo "usage: $0 releaseLabel metadata_date.tsv.gz epiToPublicAndDate.date"
 }
 
-if [ $# != 4 ]; then
+if [ $# != 3 ]; then
   usage
   exit 1
 fi
 
 releaseLabel=$1
 nextmeta=$2
-nextfasta=$3
-epiToPublic=$4
+epiToPublic=$3
 
 scriptDir=$(dirname "${BASH_SOURCE[0]}")
 ottoDir=/hive/data/outside/otto/sarscov2phylo
 problematicSitesVcf=/hive/data/genomes/wuhCor1/bed/problematicSites/20-08-26/problematic_sites_sarsCov2.vcf
 genbankFa=$ottoDir/ncbi.latest/genbank.fa.xz
 cogUkFa=$ottoDir/cogUk.latest/cog_all.fasta.xz
 cncbFa=$ottoDir/cncb.latest/cncb.nonGenBank.fasta
 
 mkdir -p $ottoDir/$releaseLabel
 cd $ottoDir/$releaseLabel
 
-$scriptDir/getRelease.sh $releaseLabel $nextmeta $nextfasta >& getRelease.log
+$scriptDir/processGisaid.sh $releaseLabel $treeDir $msaFile $nextmeta $problematicSitesVcf \
+    >& processGisaid.log
 
-$scriptDir/processRelease.sh $releaseLabel $problematicSitesVcf >& processRelease.log
-
-$scriptDir/mapPublic.sh $releaseLabel $problematicSitesVcf $epiToPublic >& mapPublic.log
+$scriptDir/mapPublic.sh $releaseLabel $nextmeta $problematicSitesVcf $epiToPublic >& mapPublic.log
 
 $scriptDir/extractUnmappedPublic.sh $epiToPublic $genbankFa $cogUkFa $cncbFa \
     >& extractUnmappedPublic.log
 
 $scriptDir/addUnmappedPublic.sh $releaseLabel >& addUnmappedPublic.log