Package information

• pandas 1.0.5
• biopython 1.77

Package information

• ggplot2
• hrbrthemes
• dplyr
• tidyr
• viridis

All datasets for the analysis are avaliable at https://zenodo.org/record/4007666#.X0qEK3UzacC

Data_filt.py : The python script performs the data cleaning and writes the cleaned data to a new file. This code gets rid of NAs and unknowns from the clinical dataset. The dataset originally contained 4592 samples which after running the python script were 3620.

• Input: The script takes input the GISAID clinical data containing all samples, along with the complete GISAID genomic data for 4592 samples.
• Output: Two files containing filtered samples.

Data was filtered and classified into Best Reaction, Mild Reaction and Worst reaction. The age was made a numerical variable and then a plot was made against status. This was to see the immunity distribution across ages

The phylogenetic analysis was performed on the https://usegalaxy.org/ server using the filtered genomic datasets for 3620 genomes.
MAFFT was used for multiple sequence alignment with default paramteres, and FASTTREE was used for the phylogenetic tree reconstruction.

To generate the annotations for the tree ORCaD app was used, along with OR calculations
For OR calculations see ORCaD

Please contact [email protected] in case of any issues

The project was licensed under the Apache 2.0 license