# install.packages("BiocManager")
# BiocManager::install("msa")
#
https://www.ncbi.nlm.nih.gov/nuccore/NC_045512.2?report=fasta # SARS-CoV2
#
https://www.ncbi.nlm.nih.gov/nuccore/AF033819.3?report=fasta # HIV-1
library("msa"); library("seqinr"); library(ape)
SARSCoV <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/SARSCoV.fasta")
MERS <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/MERS.fasta")
SARSCoV2 <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/SARSCoV2.fasta")
HIV1 <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/HIV1.fasta")
HIV2 <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/HIV2.fasta")
Variola <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/Variola virus.fasta")
H1N1 <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/H1N1.fasta")
Zika <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/Zika.fasta")
Ebola <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/Ebola.fasta")
HAV <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/HAV.fasta")
HBV <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/HBV.fasta")
HCV <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/HCV.fasta")
HDV <- readAAStringSet(filepath = "D:/Work/Projects/SmallSciTheses/2019ncov/HDV.fasta")
da <- msa(inputSeqs = c(SARSCoV, MERS, SARSCoV2, HIV1, HIV2, Variola, H1N1, Zika, Ebola, HAV, HBV, HCV, HDV))
save(da, file = "D:/Work/Projects/SmallSciTheses/2019ncov/da")
da2 <- msaConvert(da, type="seqinr::alignment")
da2$nam <- c("HIV-1", "HIV-2", "SARS-CoV", "SARS-CoV2", "Zika virus", "HDV", "MERS", "Variola virus", "HCV", "Ebolavirus", "HAV", "H1N1", "HBV")
da.dist <- dist.alignment(da2, "s")
dt <- nj(da.dist)
plot(dt, main="Phylogenetic Tree of some viruses", cex = 1.5)