COVIDalignmentPhylogeny
Takes the consensus fasta SARS-CoV-2 genome assembly outputs of IlluminaPreprocessAssembly.wdl and NanoporeGuppyAssembly.wdl and a reference genome which should be stored in your Terra workspace data and does the following:
- concatenates the input fasta genome and the references into a multifasta
- aligns the sequences using mafft
- generates a pairwise distance matrix with snp-dists
- generates a phylogeny with iqtree
- exports the outputs to your chosen google
This is a companion discussion topic for the original entry at github.com/CDPHE/COVIDalignmentPhylogeny