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